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The brick masonry work over beams and lintels in 
buildings has been considered only as an overburden in 
conventional method of design. However, with the increased 
understanding of reinforced brick masonry, it has been 
established that reinforced brick masonry on reinforced 
concrete beams gives additional strength to the latter. 
Earlier investigators have established that composite 
action between the beam and the masonry over it is possible 
i.e, additional strength is available only if height of 
brick masonry above the beam is more than half the span 
and load is applied at the top of brick work, Therefore, 
no change in design procedure has been recommended for 
situations where height to span ratio is below 0,5 and 
also in situations where the load is applied at the junction 
of the beam and masonry wall. 
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Tlie present work addresses itself to the study of 
the extent of interaction between the concrete beam and 
the masonry above it by considering tv/o types of loading 
( i) system loaded at the top of brick work (2) system 
loaded at the ijunction of concrete beam and masonry wall 
which is of common occurrence in multis ooreyed framed 
construction. Tlirough out the e:cperimental investigation, 
the supporting beam was a precast concrete member, 8.0 cm 
thick, in 1:2:4 nominal mix except in tvTO samples where 
thickness vra.3 increased to 16 cm . Single legged Z shaped 
shear cum tensile connectors were p^rovided through out the 
span at a uniforcni spacing to make the concrete beam and 
the brick panel act as a composite structure. Parametric 
study has been made to study the effect of following 
parameters on the ultimate load-carrying capacity of 
composite structure. 

(1) Variation of cement sand ratio in mortar for 
brick T^rork. 

(2) "ariation of height to span ratio. 

(3) Location and size of symmetric and non-symmetric 
openings in the brick ma,sonry. 

In order to conduct the above study, 30 full 
scale load tests have been conducted in the laboratory 
at I.I.T. Kanpur using local materials. Mechanical 
properties of brick v/ork, concrete and reinforcing steel 
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has 1)0311 established rad reported in the nresent work. 

Some of the Woll^boam inbocaction tests have been conducted, 
loaded from top of brick masonry while others loaded at the 
junction of concrete beaia and brick masonry. The experimental 
’results obtained have been verified by analytical study, 

For the purpose of analytical study, the composite structure 
has been idealised as a plane stress problem and solved by 
finite element method. A linear strain triangular element 
has been used for concrete and brick work while linear 
strain b-ar element has been used for reinforcement. The 
sub'-regioning of linear strain triangular element and 
linear strain bar elements have been carried for the purposes 
of accuracy of results and economy in computations. The 
ultimate load carrying capacity of the composite structure 
has been obtained by analysing cracking of concrete and/or 
brick work and yielding of concrete, brick work and steel. 

The uniaxial stress strain laviT for concrete is represented 
by an idealised bilinear curve. Brick vkirk is assumed to 
be elastic and orthotrojDic upto yield point beyond which 
its behaviour is assumed to be perfectly plastic. The 
steel is idealised as an elasto-plastic material. Maximum 
stress theory is used to check for cracking. Yonmises 
yield criterion is used for yielding of concrete and Hill*s 
anisotropic yield criterion is used for brick vrork. 
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A computer progr amm e has been developed, making 
use of criteria described above to trace the load deflection 
response and crack propagation from zero load to ultimate 
load. The salient conclusions from the present study are: 

(1) Single legged Z shaped shear cum tensile connectors 
have been found adequate to make the structure 
behave as a composite structure. 

(2) The entire structure behaves as an under- reinforced 
beam. 

(3) Cement sand mortar of 1:6 ratio has been found 
adequate incase the load is applied at the top of 
brick work, while a rich mortar of 1:3 cement sand 
ratio is necessary when the load is applied at the 
junction of concrete beam and the brick panel 
above it. 

(4) Openings at the centre of wall do not affect the 
load carrying capacity of the structure but an 
opening near the support affects the load carrying 
capacity adversely for the case of the load applied 
at the top of the system. However, if the load is 
applied at the junction of concrete beam and the 
wall, both types of openings affect the load carrying 
capacity. 
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finite clement model used in tlie present work 
predicts ceflection, ultimate load capacity and 
crack propagation quite accurately for tke case 
of load applied at the top of i^rick work. On the 
other hand this mor'ol predicts accurately ultimate 
load hut fa].ls to predict the deflection and crack 
propagation accurately in case when the load is 
applied at the ounction of concrete beam and brick 


masonry wall 



INTRODUCTION 


Berlins and lintels are important structural components 
of buildings. They contribute appreciably towards the cost 
of building. This is more so in case of framed construction. 
The lintels above ail the doors and windows, the verandah 
beams and the grade beams support the masonry upto the ceiling 
level. In framed structures, this masonry work serves as 
partition walls. The brick masonry work over beams and 
lintels in buildings has been considered only as an over- 
burden in the conventional method of design. However, with 
the increased understanding it has been established that 
brick masonry on reinforced concrete beams behaves differently 

f 

and gives additional strength to the composite system, 

lilarlier investigators observed that composite action 
between beams and masonry is possible i,e, additional strength 
is available only if the height of brick masonry above the 
beam is more than half of the span and load is applied at the 
top of brickwork. Therefore, no change in' design procedure 
is recommended for situations where height to span ratio is- 
le ss than 0.5. 

The present work addresses itself to the study of 
interaction between the concrete beam and the masonry above 
it, using single legged Z-shaped vertical connectors of mild 



steel b "srs embedded a,lon.g the length of the composite 
construction. Two types of loads have been considered in 
the preseirc study both being inplane loads applied at (l) the 
top of brickwork henceforth called the compressive loading and 
(2) the junction of concrete beam and masonry wall henceforth 
called the teasile loading. The ultimate load carrying 
capacity of the composite structure is obtained experimentally 
and verified analytically. The crack propagation with 
increased load is also studied. 

Existing work on brick-masonry and brick-masonry 
supported on R.C. beams has been reviewed and presented in 
Chapter I section 1.1. The work on analytical tools to study 
the behaviour of such composite systems is also briefly 
reviewed in this chapter in section 1.2. The scope of 
present investigation is described in section 1 .3 of this 
very chapter. 

In any investigation of this kind the importance 
of material properties can itot be over emphasized, while 
the material properties of reinforced concrete are by and 
large well-known and can be easily “'established for a given 
situation, the same is not true for brick-vi/ork. The material 
properties of brick-masonry depend upon the quality of bricks, 
the mortar used and, the last but not the least, the work- 
manship. Quantifying the mechanical material properties for 
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brickwork and. establishing suitable constitutive laws for 
brickwork, concrete and reinforcing steel is the subject of 
discussion of Cha.pter II. 

d’cperimental studies h.ave been conducted on about 
30 specimens to stucfy the effect of following parameters on 
the ultimate load carrying capacity of the composite structure, 

(1) Variation of cement sand ratio in mortar for brick- 
work, 

(2) Variation of height to span ratio. 

(3) Location and size of symmetric and non-symraetric 
openings in the brick masonry, 

Full scale load tests have been conducted in the 
Structural Engineering laboratory at III Kanpur* Ihe 
experimental details, loading arrangements and results of 
experimental studies form the subject matter of Chapter III, 

Ihe experimental results obtained in the present 
work have been verified by analytical study. For the 
purposes of analytical study, the composite structure has 
been idealised as a plane stress problem and solved by 
finite element method, A linear strain triangular element 
has been used for concrete and brickwork while linear strain 
bar element has been used for reinforcement. Ihe sub- 
regioning of linear strain triangular elements and linear 
strain bar elements have been carried out for purposes of 



accuracy in identifying the, cracked region. The ultimate 
load carrying capacity of the composite structure has been 
obtained by analysing cracking of concrete and/or brickwork 
and yielding of concrete, brickwork and steel. The uniaxial 
stress-strain law for concrete is represented by an idealised 
bilinear curve. Brickwork is assumed to be elastic and 
orthotropic upto yield point beyond which its behaviour is 
assumed to be perfectly plastic. The steel is idealised as 
an elastoplastic material. Maximum stress theory is used to 
check for cracking. Von Mises yield criterion is used for 
yielding of concrete and Hill's anisotropic yield criterion 
is used for brickwork. These are discussed in Chapter IV, 

1 computer programme in Fortran IV has been developed making 
use of the model developed and criteria described in order to 
trace the load deflection response and crack propagation from 
initial load to failure load. The analytical results of 
various experimental specimens have been obtained using the 
developed programme which are presented in this Chapter, 

There is a good agreement between the experimental and analytical 
results obtained, Bifferences, if any, are discussed herein. 

The conclusions drawn on the basis of experimental 
and analytical work are s’ummarised in Chapter V, A design 
methodology for such composite systems based on the experience 
of present work is described herein. Some suggestions for 
future developmental work in this regard are also indicated 
in this ChaJ)ter, 



GHiiPTER I 


REVIEW OE EXIST IM& LITER .HURE 


1.1 jJJCEERLvIEJT'ED STUDIES 

1.1.1 Studies on Brickv/ork 

JBrick iuasonry has been in use for times immemorial. 
However, use of reinforced brick masonry (R.B.M. ) dates 
back to about one .and h;ilf century. Sir Mark Isainbard Brunei 
a Chief Engineer of ITewyork was probably the first man to 
introduce iI.B.m, in the year 1813. Ir. 1836, he tested some 
structures to determine the strength contribution of reinforce 
ment in brick masonry. In 1837, Colonel Pasley conducted a 
series of tests on R.B.M.beamg* In 1851 , several' tests on R.B.M 
beams were done in London using Portland cement. Corsion 
in 1872 computed allowable stresses for use in masonry 
lintels and discussed relation of tensile strength of 
masonry to that of mortar, Hugofillipi in 1853 and Mench 
in 1919 made some tests on R.B.M. beams. Government of 
India in 1923 published a technical paper Mo. 38 with various 
test results of R.B.M, beams, columns, arches etc. This 
work done by Brener is a hall mark towards the development 
of R.B.M, in India. Pearson, Stang and Me J3urney (1) in 
1932 reported the test results on R.B.M, beams, Withey 
(2,3) also carried out some tests on R.B.M, beams and 
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columns. During the period 1922 to 1950 research was 
conducted on both plain and reinforced briorkv/ork at National 
Bureau of standards and practically all Engineering Institutions 
of U.S.i. It has been established some 45 years ago that 
height/thickness ratio of compression specimens of brickwork 
is important. In order to obtain data on the effect of 
varying the strength of brick and mortar, very large' mumber 
of tests have been conducted in Building Research Station 
England by Thomas, E.G. (4) on brick piers, 9 inches squares 
and 36 inches high. He showed that strength of brickwork is 
a function of both i.e. the strength of bricks and 
of mortar used. Results obtained by him further indicrj,te 
that the effect of increasing mortar strength, in teims of 
the resulting brickwork, is not proportional. In 1973 
Purshothaman (5) carried out some tests at I.I.T. Kanpur 
on brick masonry and also determined the properties of 
brickwork, analytically using finite eleaient method. There 
is vast difference between the experimental and analytical 
results obtained by him. The Structural Clay Products 
Institute (SCPI) (6), has recoimuended that brickwork prism 
tests be conducted on specimen not less than 12 inches in. 
height and shall have a height /thickness ratio of not less 
than 2 and not more than 5. The Institute recommends that 
modulus of elasticity,E ,for brickwork be taken as 1000 f*m 
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where f'm is the 28 day compressive strength of bride 
prisms of h/t ratio of 5. '2he shear modulus is predicted 
to be 400 f 'm. 

1,1,2 Studies on Y/alls Supported on Beams 

T/ood R.H, (7) was the first to conduct a series of 
experiments on walls supported on beams. Tests on 9'' 
bride walls without supporting beams showed that even these 
unsupported walls can resist large vei’tical loads. When 
reinforced concrete beams were used with brick walls 5 
tension concentration in supporting beams was noticed. 

Walls v^ith door and window openings were also tested. Wood 
found that deep girder action did not apply to loading at 
the junction of concrete beam and masonry wall, unless 
tensile connectors could be placed between wall and beam. 

He suggested that some light reinforcement may be necessary 
in the walls to take care of shrinkage and also continuity 
in case of continuous walls. [Jltimate load could not be 
reached in these tests, 

Rosenhaupt and his associates (8,9,10,11,12) studied 
the behaviour of walls on continuous beams, effect of various 
types of openings with and v/ith out stiffening concrete 
frames surrounding the openings and the effect of prestressing 
the brick\i/ork, besides the effect of foundation seutlement. 
They found that reinforcement in supporting beams had secondary 
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effect, except when flexure controlled the failure as in 
the case of shallow walls or strong brickwork* On the other 
hand shear, bond and vertical compressive stresses in deep 
walls controlled failure. Inclusion of horizontal and 
vertical concrete ties in the beam wall system had significant 
effect on the internal stress distribution. The openings 
definitely weakened the structure but they showed that 
strength could be compensated by prestressing or by the 
introduction of vertical concrete ties. T/hile testing 
continuous composite walls, it was found that positive 
moments existed even over the middle supports. The failure 
modes peculiar to such walls and beams are the crushing of 
masonry over the supports and effects of shear in masonry 
panel. 

-Burhouse (13) used concrete beams and encased joists 
as supporting beams. He introduced a building paper joint 
at the mid span as previous investigators had found the 
separation of supporting beam from the masonry wall. 

Crushing of brickwork was found to be the predominant mode 
of failure and strong brickwork is advocated if the beneficial 
arching action in walls on beams is to be fully utilised. 

In shallow beams, progressive slipping of w^ais on beams 
was noticed, showing the need for an adequate connection 
between the v/alls and the beams. 



Purshothainan (5) coaducted a series of experiments 
at ITT ICanpur. He studied the effect of mortar ratio in 
briclov/orlc, effect of height to span ratio and effect of 
openings. Eo concludes that the beneficial interaction is 
available oiiLy if the load is applied at top of bricla’/orlc 
With SAmraietric door and v/indov/ openings. It is further 
concluded that the interaction exists only if rich mortar 
of 1:3 cement sand ratio is used. No interaction exists if 
the opening is eccentric or if the load is applied at the 
junction of beam and the wall. It is further recommended 
that for design purposes, the benefit of interaction should 
be used only if the height to span ratio is more than 0,5. 

Prasad Rao and Mallick (I 4 ) conducted tests on 
walls-on-beams giving particular attention to their ultiE.ate 
strength, 

ihuflesh et al (15) also tested the wa-ll on beams. They 
pp»vided tensile connectors in the form of two legged 
stirrups, close near the supports and wide apart at the 
midspan. It is concluded that by this arrangement, the 
load carrying capacity of the composite structure is 
increased and claimed a saving of 20 percent in reinforcement. 

Smith .6, Stafford (16) conducted a series of experiments 
to study the collapse of masonry walls on steel beams and 
reinforced concrete beams. Six tests were conducted on 
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wall beams 6 ft span and 4 ft Mgh, They concluded that 
the behaviour of two types of wall beam structure differ 
to the extent that modifications to the design method are 
necessary. They also proposed the recommendation for changes 
in design method, 

1.1.3 Review of I.S, Code of Pr^ictice 

ictually spealcing there are no codes of practice 

which deal specifically v/ith the problem of wall on beams. 

National Building Code of India 1970 has made an indirect 

use of the arch action for the design of capping beams over 

under reamed piles. It has suggested the use of niaximurn 

T7l2 

bending moment of where Yl is the uniformly distributed 

load per cm run and u is the effective span in cm provided 

the beams are supported during construction. If it is not 

supported during construction, an increased value of “7— for 

fu 

the maximum bending moment has been recommended. The minimum 
depth of capping beams recommended is 15 cm. 

1.2 MllYTIOiiL studies 

Earlier investigators have carried out analytical 
investigations on the interaction betvreen brick walls and 
their supporting beams by classical theory of elasticity, 
latice analogy, beam truss and virendel girder aiifalogy. The 
first simplification for the state of stress in a wall on 
beam was given by Wood (7). Based on the observation that 
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most of the compressive load was transferred to the supports 

through arch action, he suggested that the beams be designed 

WL , YiQj j -u • 

for a reduced bending moment of 7:77- to depending upon the 

50 100 

position of openings. An equivalent moment: arm of 2/3 of 
depth subject to a maximum of 0.7 times the span was also 
suggested to determine the area of reinforcement. He further 
suggested thcvt the limiting moment arm method be used in 
walls without openings. For smaller depths, arch action 
may be absent and hence no change in design moments was 
suggested. Rosenhaupt (9) also su';gested the equivalent 
moment arm method using 0,6 times the height of wall. Rosenhaup 
suggested the uses of truss analogy to deal with openings 
in wall. 

Prasad Rao and Mallick (14) have discussed several 
simple methods of computing the ultimate strength of walls 
on be£Ujis without opening. I'he simplest is to use an 
equivalent lever arm of 0.75 to 0.85 times of the depth and 
multipljr the sane by the area of reinforcement and yield 
stress of steel to obtain the ultimate moment. I'his approach 
is valid for tension failure. For compression failure 
Whitney’s theory could be used. 

Ghandrashekharan, K. and Abraham Jacob, K. (17) 
carried out two dimensional photo elastic analysis to study 
the composite action of walls supported on beams. !rhe3'' 



constructed the photoelastic composite model of t?/o different 
materials which p-awe the required ratio of elastic constants 
of wall and beam at an elevated temperature of 115°C!. The 
tvro materials used by them were 6olumbia resin OR 39 and 
Araldite CY 230, OR-39 was used to represent the beaiJ. and 

OY -230 represented the wall portion. The ratio of modulus 
of elasticity obtained at room temperature for OR 39- 
^iraldite composite was 1.32 while at elevated temperature 
(l15°c) this ratio was 23.5. They have shown that the inter 
face stresses can be obtained directly using the photoelas tic 
data along with continuity conditions at the interface. The 
stresses at the interior of wall portion were obtained by 
using ntimerical technique which requires only the boundary 
stresses to be known. They considered walls with and without 
openings • Their results compared well with those of 
Rosenhaupt ( 9 ) • 

1.2.1 Finite Element Study of Reinforced Concrete Structures 
and Brick Wall on Concrete Beams 

In recent years, the finite element technique has 
successfully been applied to perform the non-linear analysis 
of reinforced concrete structures. However, its use for 
brick masonry structures has been very limited. 

Scordelis (I 8 ) in his state of art report has 
reviewed the present status of research and application of 



finite element method to the analysis of reinforced concrete 
structures. The application of finite element technique to 
reinforced concrete was first reported by Nago and Scordelis( 19) . 
They preformed linear elastic analysis on simple beams with 
predefined crack patterns to determine the principal stresses 
in concrete, stresses in steel reinforcement and bond stresses. 
They also took into account bond slip by use of finite spring 
elements designed as bond links between steel and concrete 
spaced along the bar length, l^go, Scordelis and Franklin (20) 
used the same approach to study the shear in beams with diagonal 
tension cracks considering the effect of stirrups, dowel shear, 
aggregate interlock and horizontal splitting along reinf orcement 
near the support. iTilson (2l) extended the technique further 
by adding non-linear material and non-linear bond slip 
relationships. He modelled the cracking of concrete by 
separation of nodes. Thus propagation of cracks in the 
structure was constrained and cracks developed only along 
the inter element boundries. Moreover, nodal separation 
amounted to a continuous change in the structural topology 
and nodal connectivity of finite element mesh. 

Rashid (22) and Franklin (25) incorporated the 
cracking of concrete and non-linear material properties in 
the analjrtical model by altering the elasticity matrices 
of individual elements. Incremental loading, with iterations 



within each increment, were used to account for cracking 
of the elemeats and non-linear properties of materials* 
frai-iklin made use of special frame type elements, quadri- 
lateral plane stress elements, azial bar elements. He 
analysed the fr.imes with and without infils. Zienkiewicz 
et al (24) carried out the analysis of rock as 'no tension' 
material by finite element method. A linear elastic analysis 
was conducted to determine the stresses in the elements. 

Since the rock was assumed as a no tension material, tensile 
stresses from the elements were released and equivalent 
nodal forces due to stress release v/ere computed. The 
structure was reanalysed for these nodal forces and the 
process was repeated until convergence was obtained. 

Falliappan and Hath (25) assumed the concrete to initially 
exhibit a limited tensile strength which was reduced to 
zero once cracking occurred. An incremental loading procedure 
was employed and redistribution of stress due to cracking 
of concrete was achieved by iterative 'stress transfer 
approach', Zienkiewicz et al (26) presented on 'initial stress' 
finite element approach for the solution of elastoplastic 
problems. The pseudo load vector due to change in material 
properties was computed and material non-linearity/ was' 
included by applying and iterating these pseudo loads, 
Valliappan and Doolan (2?) applied in 'initial stress' finite 
element method to study the behaviour of reinforced concrete 



structures which included tensile cracking and yielding 
of conci'ete and steel. 'They made use of constant strain 
triangular elements for concrete and bar elements for 
reinforcement , 

Mufti et al (28) studied the non-linear behaviour 
of structural concrete. Plane stress triangular ;ancl rectangular 
elements -were used. The superiority of rectangular elements 
over triangulan elements for nonlinear analysis of reinforced 
concrete structures was shown. Concrete was assumed to be 
linear elastic in tension and nonlinear stress strain relation- 
ship was used in compression. Steel elements were connected 
to concrete elements through bond linkage elements. 

Purshothaman (5) developed a finite element model 
to predict the behaviour of brick masonry walls over reinfor- 
ced concrete beams. In this model, brickwork was idealised 
as a linearly elastic orthotropic material and concrete and 
steel as elastic perfectly plastic material. Von Miscs 
criterion was used for all the three materials. Bond slip 
was neglected, do slip between brick wall and concrete 
beaiTi was assumed. Elastic, post cracking and failure stages 
were incorporated in the model. Constant strain triangular 
elements for brickwork and concrete and beam elements for 
steel were used. Incremental iterative method modifying 
the stiffness of critical elements, was used to trace the 
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progressive craclcing, local compression fc-alure and yielding 
of steel, Suida’i and Schnorbrich (29) performed the three 
dimensional finite element analysis of reinforced concrete 
beams, fhree dimensional 20 noded isoparametric elements 
were used in the analysis, dlastoplastic beha\riour of 
concrete and reinforcing steel and cracking of concrete 
were accomodated in the finite element model. The results 
obtained compared well with those obtained experimentally. 

Colville and ibbasi (50) presented a finite element 
approach for nonlinear analysis of reinforced concrete plane 
stress problems where in the reinforcing steel was included 
in the stiffness formulation of element . Constant strain 
triangular elements and rectangular elements with linear edge 
displacements were used in the analysts. ITonhomogeneity and 
anisotropy due to reinforcement and effect of tension cracking 
were considered in the model. The initial stress approach was 
used to predict the extent and location of tension cracks. 
Houde (31) has proposed nonlinear bond slip, aggregate inter- 
lock and dowel action relationships, and has used these to 
study the propagation of cracks in simple beams. Mirza and 
mufti (32) have employed these relationships to analyse a 
reinforced concrete bracket and a beai'ii column ;]0iiit, Mam, 
Chung-Hyum and Salmon (33) studied the nonlinear behaviour 
of reinforced concrete beams under short term loading. They 
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used combination of incomp. stable isoparameteric quadri- 
lateral element and a linear bas element. Effects of 
progressive cr:i,clcing and yielding of both steel and concrete 
were included in the analysis. Incremental iterative method 
with constant and variable stiffness was used to study in 
nonlinear behaviour. In correctness of constant stifTness 
method in the evaluation of nonlinear behaviour of reinforced 
concrete due to cracking was shown, 

failure of brittle materials due to crack propagation 
was studied by Hag id and Hashimi (34) using finite element 
method, in incremental finite element method was used to 
detect the origin of initial crack and to trace the formation 
and pi^opagation of subsequent cracks upto failure. Two 
types of elements, simple constant strain triangular 
elements and eight noded isoparametric quadrilateral elements, 
were used. Cracking was accounted for by separating the 
nodes when an element indic.ated the tensile failure. To 
avoid termination of analysis and renumbering of nodes 
after cracking and to allow for automatic continuation of 
analysis after cracking, a multiple numbering was adopted in 
the solution procedure, failure loads obtained by using this 
method were found to be higher than those obtained experimentally 

Oedolin and Poli (35) conducted finite element studies 


of shear critical reinforced concrete beams, A nonlinear 
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model was developed to trace the history of strains, stresses 
and crack propagation in reinforced concrete beams subjected 
xo plane state of stress for a monotonic increase of external 
load. Concrete was assumed to be an ortho tropic material in 
the direction of existing principal stresses. The stress 
transfer between steel and concrete was achieved through 
experimentally determined nonlinear bond slip relationship. 
Cracks in concrete were assumed capable of transmitting some 
shear parallel to crack it self. Analysis was performed 
using constant strain triangular element and an incremental 
iterative procedure based on tangent stiffness approach. 

Cedolin et al (36,37) presented a finite element analysis 
of deep beams and prestressed concrete beams. The analysis 
included nonlinear constitutive relation of concrete, presence 
of main and web reinforcement, their relative movement with 
respect to concrete and crack propagation in concrete along 
with other related phenomena such as aggregate interlocking 
and dowel action. Constant strain triangular elements were 
used to idealize concrete. One dimensional linlcage elements, 
parallel to steel bars, with nonlinear bond slip characteristic, 
were used to oonnect the steel and concrete elements. An 
incremental iterabive method was used to predict the history 
of stresses, deflections and crack propagation upto failure, 

Paneerselvam (38) presented a nonlinear finite 
element analysis of reinforced concrete framed structures. 



A reini'oroed concrete element using 4- noded isoparameteric 
quadrilateral with incompatable modes wad developed. This 
element included reinforcement in any orienbation* Sub- 
regioning of elements was carried out for the computation 
of pseudo loads due to nonlinearities. The uniaxial curve 
for concrete was idealised by parabola-rectangle. Octahedral 
shear stress and shear strain criteria were used for biaxial 
yielding and crusliing of concrete, Nonlinearities due to 
nonlindar stress strain law for concrete, cracking, yielding 
and crushing of concrete and yielding of steel were included 
in the analysis, 'Initial stress' approach for incremental 
iterative procedure was used, A,B, Agarwal (39) studied 
the nonlinear analysis of reinforced concrete planar structures 
subjected to monotonic reversed cyclic and dynamic loads. 

He idealised ube concrete and steel as elastoplas tic material. 
The finite element model used accounted for material non- 
linearities developed due to cracking, yielding and crushing 
of concrete and plasticity of steel reinforcement, Rectangular 
plane stress element, with three degrees of freedom at each 
node was used. Steel reinf orcemenr was assumed as a un- 
axially stressed material, and smeared or uniformly distributed 
over the concrete in an element. The results obtained using 
this model compared well with experimental results. 



1.3 OBJWOT iUTO SCOPE OF PRESEIW IRVESl'IG-.l'i'IOE 


PEe preceodirLF li'beratare survey reveals tliat 
exper iiaeru: 'ition on masonry walls supported on concrete 
beams shows interaction between them only when the height 
of masonry vrall suppoj’ted on R.O. be<xfn is more than half the 
span and the load is applied at the top of brickwork. I'his 
interaction is basically due to arch action which helps to 
reduce the dimensions of concrete beam and the reinforcement 
provided therein. There is no accepted guideline for this 
reduction as expemnental results of various investigators 
show different reduction coefficients. The reason for this 
variation is that the extent of arch action depends upon 
the properties of brickwork used for masonry wall, ill the 
investiga/bors have agreed that little or no interaction is 
there in c ise of height to span ratio below 0.5. Eurthorinore , 
when the load is applied at the junction of concrete beam and 
the masonry wallj it is reported that no interaction exists 
for any height to span ratio. To achieve interaction in 
such situations Vfood (7) has suggested the use of some suitable 
shear connectors. 

The present work has been motivated to study the 
interaction between brick walls on R.G. beams using single 
legged Z shaped vertical connectors of mild steel bars 
embedded along the length of composite construction. In the 



t- I 


event of such an interaction, the thickness of R.O, beam 
can be reduced to bare ninimum (sufficient only to provide 
cover uo tensile reinforcement). Moreover, composite iciion 
shall tend to reduce the requirement of tensile reinforcement 
thereby reducing the cost of construction as compared to 
conventional practice. Therefore full load tests on 30 
specimens of brick masonry walls supported on thin reinforced 
concrete beams having a constant span of 3*25 metres have 
been conducted upto failure and their load response charac- 
teristics, crack propagation and ultimate load carrying 
capacity have been recorded. The studies have been made by 
loading the specimen either at (1) the top of brick masonry 
or (2) the ymction of concrete beam and masonry wall. In any 
effort of this kind, determination of mechanical properties of 
brickwork using local raateri.al is inevitable. The same has 
been c irried out in the present work in order to achieve 
rational results. To study the variation of interaction of 
such composite construction parametric study has been conducted 
by varying the following parameters, 

(1) Gem-nb sand ratio in the mortar used for brick work, 

(2) Heifjht of brick masonry , 

(3) Location and size of symmetric and nonsymmetric 

openings in brick masonry. 

The experimental results obtained have been verified 


by analytical study, kor the purpose of analytical study, the 
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composite structure has been idealised as a plane stress 
problem and solved by finite element method, A linear strain 
triaiigular element has been used for concrete and brickwork 
while linear strain bar element has been used for reinforcement. 
The sub-re gioning of linear strain triangular element and 
linear strain bar elements have been carried for the purposes 
of accuracy of results and econoxay in computations. The 
ultimate load carrying capacity of the composite structure 
has been obtained by analysing cracking of concrete and/or 
brickwork and yielding of concrete, brickv 7 ork and steel. The 
uniaxial stress strain law for concrete is represented by 
an idealised bilinear curve. Brickwork is assumed to be 
elastic and orthotropic upto yield point be37'ond which its 
behaviour is assumed to be perfectly plastic. The steel is 
idealised as an elasto -plastic material. Maximum stress 
theory is used to check for cracking. Von Mises yield 
criterion is use' for yielding of concrete and Hill's 
anisotropic yield criterion is used for brickwork. 

A computer programme in fortran IV has been developed, 
making use of the model developed and criteria described above 
to trace the load deflection response and crack propagation 
from initial load to failure load. 



CHaL-'TER II 


MATBRIilL PROPBRTIES 


2. 1 INTROBUGTION 

Brick walls supported on R.C, beams consists of 
three main materials i«e, concrete, brickwork and steel 
reinforcement. Properties of steel reinforcement are 
generally well defined as it is comparatively close to ideal 
material. It is the hetrogeneity of concrete and brickwork 
which makes the determination of the constitutive relations 
for these materials a difficult task. The structural 
response of reinforced concrete and brickwork is a function 
of the properties of component materials apart from other 
parameters. The realistic determination of the response 
of reinforced concrete and brickwork structures requires 
the knowledge of ':he inelastic behaviour of the component 
materials as well the ability to incorporate these into 
the rational analysis of real life structures. The properties 
of reinforced concrete are also by now well documented. 
However it is the brickwork whose properties are not so 
readily available, primarily because they vary widely from 
place to place. This chapter deals with the determination 
of mechanical properties of brickwork. The material 



H 


constitutive relations used to carry out nonlinear finite 
element analysis of composite system are also described herein, 

2,2 EXPERIivlHJiTfll DETERMM’IOi^ CO? li/IBCHAlIG^iL PROFERTIBS OF 
BRIGRWORiC 

Mechanical properties of brickwork, for example, 
compressive strength, modulus of elasticity and resistance 
to diagonal tension vary greatly, depending upon the properties 
of constituent material, the arrangement of brick layers, 
thickness of mortar beds, workmanship and direction of loading 
etcetra. The properties of bricks change from country to 
country. Furthermore, in a vast country like India, the 
properties of bricks vary even from place to place. Therefore 
mechanical properties of brickwork have been determined 
experimentally in the present work and are reported in the 
following section. 

Before embarking on the determination of mechanical 
properties of brickwork, it is imperative to determine the 
physical and mechanical properties of bricks and mortar 
used. This is the subject matter of study which follows, 

2.2,1 Bricks 

Bricks available at Kanpur, India are hand moulded, 
sun dried and burnt in country kiln with coal or fire wood, 

A consignment of 30,000 bricks was commercially procured 
for the present investigation. Out of each stack of five 





thousanc bricks, fiTe bricks were picked up in a random 
manner, Tlius a iotal of thirty bricks were tested to get 
the physical and mechanical properties of bricks through 
the following tests. 

2.2. 1.1 Dimensional tolerences 

Length, breadth and depth of each of thirty bricks 
was measured at four cross-sections to get the idea of 
average dimensions of the bricks. The results are given in 
Table 2,1, Simple statistical analysis showed that the 
length, breadth and depth of bricks had coefficient of 
variation of 0,88 percent, 1.45 percent and 1.43 percent 
respectively. The coefficient of variation in area and 
volume turns out to be 1,74 percent and 2,12 Percent 
respectively as given in Table 2.2. 

2.2. 1.2 Density and water absorbtion test 

All the specimen bricks were weighted to get the 
dry density of bricks. These were subsequently soaked in 
water for 24 hours, wiped and weighed again to get the 
water absorbtion. The results are shown in Table 2,2, 
Analysis of these results gives the coefficient of variation 
of dry density and percentage water absorbtion as 3.89 
percent and 30.1 percent respectively. 
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2 . 2 . 1 .J Oompressive strent^th 

After contacting the water absorbtion test, frog 
of all the thirty bricks was filled with 1:5 cement sand 
mortar and both surfaces top and bottom lerelled with 
this very mortar. The bricks, ciared for 24 . hours as per 
I.S. code, were finally sub 3 ected to compression test on 
a universal testing machine. The results are shown in 
Table 2.5, The average compressive strength of bricks 

p 

is 270,55 kg/cm and the coefficient of variation is 11,74- 
percent which shows that strength of bricks even in same 
consignment vary greatly. 

2,2.2 Cement Sand Mortar 

7 cm cubes of cement sand mortar were prepared for 
each of the five cement sand ratio i.e. 1 : 5 , 1 :4> 1^5, 1:6 
and 1 :o. The mortar for the preparation of these cubes 
was taken either .rom the mortar prepared for making the 
brick prisms or that used for the masonry work over the 
concrete beams. After curing for about 28 days, these 
specimens were tested in compression testing machine. The 
results are shown in Tables 2.4-, 2.5, 2,6, 2,7 and 2.8 
respectively. The mean compressive strength for five 

2 

mortars are 228.57, 185.77, 66.44, 54.92 and 25.55 kg/cm 


respectively. 
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2.2.3 Oompressire Strength and Modulus of Elasticity 
of Brickwork 

Since brickwork is an orthotropic material, its 
compressive strength and modulus of elasticity is required 
along two mutually perpendicular directions. Hence two 
types of masonry prisms have been tested (l) cast with 
horizontal mortar layers Big. (2,1a) and (2) cast with 
vertical mortar layers Eig. (2.1b). Henceforth, the brick 
masonry prisms with horizontal mortar layers are termed 
as type H and those with vertical mortar layers as type V. 

2. 2. 3.1 Breparation of specimen 

A total of sixty brick masonry prisms, six each 
of type H and Y for each of 1:8 and ranging from 1:6 to 
1 :3 cement sand mortar were prepared. All the specimens 
were cast by two skilled masons under supervision. The 
time taken in casting each prism was two mason hours. The 
dimension of prism for type H and type Y respectively 
turned out to be 35.5 cm x 35.5 cm x 143 cm high and 35 cm x 
3& cm X 435 cm high. The thickness of joints was kept 
constant at 1 cm through out the present work. 

2.2. 3. 2 Testing arrangement 

Compressive strength test was carried out on a 200 tonnes 
universal testing machine. A 1 2 mm thick rubber sheet was 
placed at bhe top and the bottom of every prism. A mild 
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steel plate of 40 cni x 40 cm x 8cm thick was placed at the 
top over the rubber sheet, This plate projected equally 
on all the four sides beyond column edges. At the top of 
steel plate, load was applied through a circular compression 
tool of 25cm diameter. For measurement of strains, a 
mechanical dial gauge of 0,01 mm least count was fixed 
between the two platfoims of universal testing machine, 

Ihe load was applied gradually in interval of 2,5 or 5.0 
tonnes depending upon the total load. In order to account 
for the compression of rubber sheet, a separate load defle- 
ction measurement was made by keeping the two rubber sheets 
one over the other covered by the same mild steel plate. 

The load was again applied by the same compression tool 
the same dial gauge was used to note the deflections. 
Finally, a curve was dravm between stress versus deflection 
for the combination of rubber sheet and mild steel plate. 

This deflection was subtracted from the deflections obtained 
for the assdmbly of brick prisms with rubber sheet and mild 
steel plate to obtain the net deflection of the brick 
masonry prism. From these deflect ions strains were calculated. 
To get the mean stress strain curve, stresses were inter- 
polated at unit interval of strains. It is assumed that 
strains go onincreasing without any increase of stress 
beyond the failure of prism due to splitting. Table 2,9 
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shows a specimen detailed computation of strains and 
stresses from the observed data for a brick prism of 1 :3 
cemenr so,xid murtar. 

2.3«3.3 lest results 

(i) Compressive strength: The maximum compressive 
strength obtained for thirty samples of each of type H and 
type V are given in Tables 2.10 and 2,11 respectively. 

(ii) Modulus of elasticity: The secant modulus of 
elasticity has been calculated for each brick prism separately 
at ten, twenty, forty, sixty and eighty percent of mean 
compressive strength. These values are given in Table 2,12 
for type H and Table 2,13 for type Y brick masonry prisms, 

(iii) Stress strain curve: Por an assumed value of 
strain, a set of three to six values of stresses were 
available (refer Table 2, lu and Table 2.11), The value of 
mean stress corresponding to a strain was obtained from 
this data for 1:6 mix using (l) arithmetic mean (2) root 
mean square. These stress strain curves are shown in 
Fig, 2,2, From this plot it is seen that the two averaging 
techniques do not give any appreciable difference. Therefore, 
only the arithmetic mean was considered for the plot of 
stress strain curves for the remaining specimen. These 
are shown in Fig, 2.3 for type H and in Fig. 2.4 for type Y 
brick masonry prisms cast in different mortar mix. For 
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finite eleraent analysis tli se curves have been idealised as 
linearly elastic upto yield limit and perfectly plastic 
af Tervvarcis . 'the idealised stress strain curves are shown 
in Eig. 2.5 for type H and in Eig. 2,6 for type Y brick 
masonry prisms, 

2.2.4 Resistance to Diagonal Tension of Reinforced Brick 

Masonry 

In order to carryout the analysis of reinforced brick 
masonry it is necessary to obtain the tensile strength of 
brickwork in bending and its resistance to diagonal tension. 
Since the actual failure in both cases takes place due to 
tensile stress being more than the tensile strength, the 
value determined in the present work, as discussed subse- 
quently, has been used for both, 

2. 2, 4.1 Preparation of specimen 

A total of 15 reinforced brick beams, three for 
each mortar defined earlier were prepared as per the 
arrangement shown in photograph affixed as Pig. 2.7. Eadh 
of these specimens was reinforced with three bars of sixteen 
millimetre diameter tor steel on tension side two bars 
oif twelve millimetre diameter tor steel on compression 
side. The dimensions of all beam specimens were 57.5 cm 
wide , 26 cm deep and 1'50 cm long. 
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2 . 2 , 4 *2 2e sting arrangement 

'tile sped, ens were tested as s-imply supported beams. 
Centre to centre distance between the supports vvas kept as 
120 cm. Two point loads were applied with the help of 
2 4000 '3b s hydraulic jack as shown in Fig. 2.7. Botn the 
load points v/ere kept at 30 cm from centre of supports. 
Deflections were measured at three points, at the centre and 
quarter spans, by the help of mechanical dial gauges, 

2, 2.4*3 Test results 

All the specimen failed near the supports in diagonal 
tension. Though the test specimen and loading was symmetrical, 
the diagonal tension crack developed on either the left or 
right support. The maximum shear stress taken as a measure 
of diagonal tension is calculated using the expression 

s 

q = — — where S is the maximum shear force, b is the width 
bjd 

of beam and jd is the lever arm. The lever arm is computed 
assuming the section to be uncracked i.e, brickwork is 
assumed to offer resistance to tensile forc^. The value of 

modulus of elasticity of steel reinforcement is taken as 

f) P " 

2.1 X 10 kg/cm and for brickwork the average of the 

corresponding value determined in previous section at 10 percent 

of average compressive stress is considered. Thus the value 

of modular ratio m varies vath mortar mix used. However, 

the value of modular ratio m was kept limited to a maximum 
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of 200, because with the increase in this value, the stress 
in compression steel exceeds the yield stress. 

The load at failure due to diagonal tension, the 
maximum deflection at mid point and the strength in diagonal 
tension for each specimen are shown in Table 2,14-, lb.® 
mean load versus mid point deflection curves are shown in. 
Pig, 2,8 which are more or less linear. 


2,3 GOdSTITUTITS RhLiTIOi^S POR BRICKWORK 
2,3*1 Brickwork in Plastic Range 

Brickwork in elastic range in compression as well 
as in tension is assumed as linear orthotropic material. 
Therefore ortho tropic constitutive law given by Darwin et.al. 
(40) has been used for the present work. 
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E , E are the modulus of elasticity in two directions and 

X y 

Vx’Uy"^^® poisson's ratio in two directions 
G- is the shear modulus 

Equivalent poisson’s ratio i) is defined as 
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Troin e iaat_ons ( 2 . 1 ) to (2*3) , we get , 
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There are four independent material constants in the 
above equations. Wiile modulus of elasticity 1 S„ and E 
can be obtained from the stress strain curves, it has not 
been possible to determine the value of the Poissons ratio,qj , 
experimentally. The value of 1} m the present work has been 
taken as 0.1435 and 0.1 as used by Jain A,K^(4.1 ) ,referring(42) .It i£ 
also difficult to determine experimentally the shear modulus, G, 
for brickwork. Grimm ( 42 ) suggested that shear modulus 
for brickwork is about half of the modulus of elasticity 
in compression. Darwon and Pecknold ( 40 ) while developing 
orthotropic nonlinear model for plain concrete derived the 
following expression for shear modulus in terms of modulus 
of elasticity and poisson's ratio 

S = ^ B.E -ZWtZJX ] (2.5) 

4( ^ ^ X y x y 

This expression heis been derived on the assumption 
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that it IS independent of the orientation of axis. The value 
of G- derived from this expression has been used in the present 
work. 

2.3,2 Yielding of Brickwork 

There is no evidence in the literature regarding 
strength envelope for brickwork in the state of combined stresses. 
Purshothaman (5) used Von Mises criterion to check for biaxial 
compression yielding of brickwork, though this criterion is 
valid for isotropic materials only. Therefore in the present 
work Hill's anisotropic yield criterion has been used as 
brickwork is orthotropic in nature before cracking. Hill's 
anisotropic yield criterion is written as 

P(^) = y 2 ^ “l 2 ^ “23^ ^ 22 " 53 ^ ■''“31 

12^'^55^23'^°^66^31 ^ ^ ^ " ^0 (2,6) 

where o- is the effective stress; are anisotropic 

0 J 

parameters, > °’22' '^33 normal stresses in the direction 

of anisotropic axes 1,2,3 and '^12’'^23^ ^31 stresses 

in planes 12, 23 and 31. 

How for plane stress problems, '^1 2^ xy 

and all other components are zero, therefore the yield 
criterion given by equation (2.6) reduces to 
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where “'11= '^12 “31 


“22 “23 “12 


The anisotropic parameters are determined from yield 
stresses in various directions obtained from independent 
tests. The initial anisotropic parameters ®22’ °''33 

a ^4 are obtained successively letting all stress components 
in yield criterion enual to zero except the one under 
consideration. Therefore from yield criterion and from an 
uniaxial test in X direction, we get 

( 2 . 8 ) 


0^- 2 


“12 “31 “11 


= 2(--^-) 


ox 


similarly, from an uniaxial t est in Y direction, we obtain 

(2.9) 

and from an uniaxial test in Z direction, we get 


®12 “23 ~ '^22 


= 2(g — ) 

oy 


Oq 2 

“23 ^ “31 ^ “33 ^ 

Similarly from a shear test, we obtain 


( 2 . 10 ) 
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where cr^^, cr^^, and are initial yield stresses 
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obtained from the uniaxial tests and is the initial 
effective stress adopted from one of the above four uniaxial 
test values. 


Ihe three unknown parameters *^23 

be obtained by solving eq^uations (2,8) to (2,10) and thus 
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tfo 2 00 2 /o 2 

= (:r-) - (T-) + fc) 


( 2 . 12 ) 

(2.13) 

(2.14) 
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Por brick element shown in Pig. 2.9, it is assumed that 

. If effective stress is adopted from the uniaxial 
test in X direction, the yield criterion given by equation (2.7) 
takes the following form on simplification (43) 
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In the finite element increment analysis the constitutive 
relation for material under going plastic deformation can be 
expressed as 

|dn' = [ Cgp] ^de] 


( 2 . 16 ) 
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where [^^p] is the elastopl'-istic matrix which is derived 
using flow rule of plasticity. The elastoplastic matrix 

A 

IS given by (26) 


H ' + J ■— [ C ] ) > 

[ bo j La J 


(2.17) 


where 0 is the elasticity matrix and H ' is the slope of 
uniaxial stress versus plastic strain curve at a particular 
value of The yielding of brickwork is not caused by 
actual plastic flow, Infact it is the cumulative effect of 
microcrack propagation which is responsible for onset of 
yielding in brickwork (44) • Moreover, the yielding of 
brickwork does not affect the overall behaviour of the 
structure to the extent the cracking of concrete, brickwork 
and yielding of steel does. Therefore, in the present v^ork, 
an 'unconstrained flow rule' suggested and used by Iin et al. 
( 44 ) has been used. This flow rule assumes that the flow of 
plastic strains is not constrained but stresses are fixed at 
the initial yield points on the yield surface. The stress 
strain relation for the unconstrained flow rule in incremental 
form IS expressed as follows. 


ido' = [ 0 ] |d -3 ^ 

where [o] is a null matrix. 


( 2 . 18 ) 
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2.3.3 Orackin^^ of Brickwork 

’’''I'ucla of literature is not available on failure 
criterion for cracking of brickwork. Pur shot haman ( 5 ) has 
used maximum stress theory to predict cracks in brick masonry 
and the same has been used in the 'present work, Iccording to 
maximum normal stress criterion, when one of the principal 
stress exceeds the tensile strength of brickwork, brickwork 
IS assumed to have cracked perpendicular to that principal 
stress. After cracking, normal stress at crack drops to zero 
and shear modulus also gets reduced. 


Pig, 2,10 shows a cracked brickwork element in 
global coordinate system x y, X' T' is a local coordinate 
system having the coordinate axis coinciding with the direction 
of principal stresses at the time of cracking. Material 
constitutive matrix in local coordinate system X’Y' is obtained 
by transforming the initial constitutive matrix from global 
coordinate system through the following transformation. 


[0-] = [Tgh [C] (2-19) 

where 0 is the initial constitutive matrix given by 
equation (2.4-) and T 2 transformation matrix given by 



Cos^P Sin^P 

Sin^p Gos^P 


-Cos P Sin p 
Oos P Sin p 
Cos^P-Sin^P 


( 2 , 20 ) 


2 SinpGosP -2CospSinp 
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Since bricl-cwork is ortho tropic in nature in global 
coordinate system XY it becomes anisotropic in any other 
coordinate system. Therefore material constiicutive matrix 
of brickwork in X'Y’ coordinate system given by equation 
(2,19) can be written as 



qi 

"I2 



[O' ] = 


^22 

^23 1 

(2.21) 
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where Cin- elements of matrix [C] « The stress 

i. J 

strain relations for brickwork in X'Y' coordinate system are 
written as 
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(2.22) 


Since element has cracked along Y' axis, stress 
perpendicular to crack would be zero for all values of strains. 
This requires that first row of [c] matrix should be zero 
arifj hence for symmetry of constitutive matrix first column 
should also be zero. Thus material constitutive matrix for 
cracked brickwqrk can be written as 
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0 

0 

0 

0 

^22 

0^3 

0 

0^3 

“b^33 


( 2 . 23 ) 


where cc^ is the shear retention factor. The shear 
retention factor has an upper and lower bound values of 
unity arfl zero. In case of brichworh the crachs are more 
uneven and due to the presence of shear cum tensile 
connectors in brichwork, value of ttig should be quite high. 
Therefore, in the present work, the upper bound value of 
unity has been used. The material constitutive matrix in 
transformed from local coordinate system to global coordinate 
system through following transformation. 




( 2 * 24 ) 


where is the constitutive matrix of the cracked 

element in global coordinate system and [ G '„] is the 

O JL 

constitutive matrix of cracked brickwork in local coordinate 
system. 


If both the principal stresses in an uncracked brick 
work element exceeds the tensile strength capacity of brick 
work it is assumed to have cracked in both principal stress 
directions. However, brickwork cracked in one direction 
can further crack along a second direction when the tensile 
stress perpendicular to that direction exceeds the tensile 
strength capacity. The brickwork cracked in two directions 
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IS assumed not to transfer any load in tension. Therefore, 
brickwork cracked in two directions is assumed to have zero 
stiffness. The constibutive matrix for such a situation 


is given by 


0 0 0 

0 0 0 ' 

I 

0 0 0 I 


(2.25) 


2.4 ll'IEGHMICil PROPERTIES OP CORGRETE 


2 . 4,1 Compressive and Tensile Strength of Concrete 

The experimental investigation carried out on 
concrete was the control tests conducted during the 
fabrication of concrete beams, 15 cm concrete cubes were 
cast to determine the compressive strength and 50 cm long, 
1 5 cm diameter cylinders were prepared to determine the 
tensile strength of concrete. Concrete mix used was 1:2:4 
by weight and a constant water cement ratio of 0,65 was 
maintained. The results of control tests are shown in 
Table 2,15 and 2.16. 


2 , 4.2 Stress Strain Curve for Concrete 

The uniaxial stress-strain curve for concrete is 
affected by numerous factors, such as shrinkage, creep 
and microcracking. In compression, its early deviation 
from a linear elastic path has mainly been attributed to 
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the micro cracking which develops at the aggregate mortar 
interface, Further disintigration and ultimate failure of 
concrete occurs due to propagation of these cracks through 
the mortar. For the quantitative description of the stress 
strain relationship of plain concrete, several emperical 
formulae are available in the literature, A good review 
of this area has been presented by Popovics (45). 

To study the behaviour of plane concrete under 
biaxial stress fields, Kupfer, Hilsdorf and Rusch ( 46 ) have 
conducted extensive experimental work. They have presented 
their experimental results in the form of a failure envelope. 
Their findings have also indicated that the strength of 
concrete under biaxial compression, only 16 percent 

larger than under uniaxial compression but the biaxial tensile 
strength of concrete is approximately equal to its uniaxial 
tensile strength, liu, Rilson and Slate (47) have also 
studied the behaviour of plane concrete in biaxial compression 
state . 

Kupfer and G-erstile ( 48 ) have used the experimental 
data of Kupfer et, al, ( 46 ), and have derived empirical 
expressions to describe the failure envelope of concrete, 
Romstad, Taylor and Herrmann (49) have also developed an 
elaborate multilinear biaxial constitutive material model 
for plain concrete. They have divided the principal stress 
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space in four damage zones. In each zone the damage level 
is assumed to be constant and the material properties are 
treated as being linear, isotropic and constant. Each 
isotropic state is described by an appropriate value of 
modulus of elasticity and the Poisson's ratio based on 
experimental evidence of Kupfer et al. (46). 

In early applications (19,27) plane concrete has 
been idealised as an isotropic and linear elastic or 
elastoplastic material. This is shown in Eig. (2,1l) 
Panneer selvam (38), adopted a nonlinear stress strain 
curve proposed by Rusch (50). This stress strain law 
shown in Eig. 2.12 is a second degree parabola-rectangle 
and is given by 

e2 

0 = f^ (e X 10^ _ X 10^) 0 < e < 0.002 

= K, f ' 0.002 < 0 < 0.0035 (2.26) 

3 c 

=0 6 > 0,0035 

The overall nonlinear behaviour of reinforced 
concrete and brickwork structures in the elastic stage 
is primarily due to tensile cracking. Hence the nonlinear 
portion (Eig. 2.12) can be idealized as a linear curve as 
shown in Eig, 2.11. This means that material non-linearity 
in the elastic range due to stress strain relationship 
can be neglected as compared to nonlinearity due to tensile 
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cracking, The computer programme developed can be used to 
accomodate any of the two stress strain laws. However, in 
the present work nonlinear stress strain law has been used. 

The stress strain relationship for concrete in 
tension is a curve corresponding to linearly elastic brittle 
material, 

2.5 00HSTITUTI7E RBIjiTIOHS FOR COHCRETB 

V. 

Concrete, in general, is in a state of biaxial stress 
condition and for a rational analysis, the behaviour and 
constitutive law under biaxial stress-state must be known, 

2,5.1 Concrete in Elastic Range 

Case 1 : Concrete in biaxial tension 

Under this condition, concrete is assumed to be an 
isotropic homogeneous material. Thus material constitutive 
matrix is given by ^ . 

il V 0 

[C] = |u 1 0 

1-lr 1 1 - 

I 0 0 — 

2 

where is the modulus of elasticity and V is the Poisson’s 
ratio for concrete. 

Case 2: Concrete in biaxial compression 

In case it is assumed that concrete follows a linear 
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stress strain relationship, the constitutive law given by 
equation (2,27) still holds good for the biaxial compression 
state of stress. If nonlinear stress strain relationship is 
used as in present v/ork then determination of modulus of 
elasticity for concrete in biaxial compression is difficult, 
since little is known about the interrelationship between 
biaxial modulii and associated Poisson’s effect, Kupfer and 
G-erstle (48) studied the nonlinear response of concrete under 
biaxial compression and proposed a biaxial constitutive law 
in terms of bulk and shear modulii, Darwin et al ( 40 ) 
considered the concrete in biaxial compression state to be 
orthotropic and proposed a biaxial stress strain law, 
Panneerselvam (38) determined the two elastic modulii corres- 
ponding to two principal strains and used the minimum of 
two to define the constitutive law in biaxial compression. 
This approach though simple, is not accurate. In the present 
work, the constitutive matrix in biaxial compression is - 
obtained as follows. 


A simple modulus of elasticity in biaxial compression 


is found out using the concept of equivalent strain given in 
( 51 ). An equivalent one dimensional effective strain 6 is 
defined as 

2 


e* = 




2 ^'^xy'" Vz'Yzx" 

( 2 . 28 ) 
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Now for plane stress condition and ‘5^ = + 8„), 

yz zx z '|_y ^ y 

Therefore equivalent strain for a plane stress case will be 


_ 


[('i - 

2(1+Ii) ^ 


y 


+ (3 


z 


+ 


(-^ -e 

^ Z X 


3 

+ 

2 


2 1/2 

I'xy 


(2.29) 


Single modulus of elasticity iiln. is found out at 

’ ^ 

equivalent strain 8 from the assumed uniaxial stress strain 
curve. Hence for elastic concrete in biaxial compression, 
the constitutive matrix [ C] is give-n by 


[C] 


E 


n 


1-i»2 


1 

V 


V 

1 


0 


0 


2,5.2 Yielding of Goncrete 


0 

0 

1-V 


(2.30) 


The failure criterion for concrete in a biaxial 
stress state was proposed by Kupfer et al. (4.6,48) based on 
their experimental investigation. The yield failure surface 
proposed by them is shown in Eig. 2,13. This yield-failure 
surface was modified by Wanchoo and l.iay (52) as shown by solid 
lines in the same figure. The modified surface assumes that 
compressive yielding is governed bjj’ Yon Mises criterion and 
associated flow rule and is slightly more conservative than 
Kuffer's surface in the biaxial compression zone. It is 
less conservative in tension compression zone but since it 
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simplifies the computational procedure it is used in 
the present work* 


The Yon Mises yield criterion in terms cf octahedral 
shear stress is given by 


p(a) 


3 

^2 "^oct 



= 0 


(2.31) 


where a is the yield stress of concre be in uniaxial case and 


oct 


1 [(a _(j )^ + (a -0 )^ + ( 0 - 0 )^ + 6(t ^ + t ^ + t 

^ ^ X Y ^ y ^ z xy yz ‘'zx^ -■ 


If the value of is substituted in the above 

equation, the yield criterion would be 


F(o') = — r(a-a)‘^ + (a-a)‘^+ (0-0}"^ + 6 (v'^ + 

^ ' Y^2^ x y ^ z X ^ xy yz 


'ZX 


) J 


1/2 


-d, =0 


For a plane stress case t 


yz 


"^zx ^z 


criterion becomes 


1/2 


P(o) = (0^ + <^ ) - 5o = 0 


( 2 . 32 ) 

0, so the above 

(2*33) 


In finite element incremental analysis the constitutive 
relation for a material undergoing plastic defcs-rmation can be 
expressed as follows (26) 


(_4 c' j = [Oep ] ) de ; 


( 2 . 34 ) 
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where [C^p] is the elastoplastic matrix which is derived 
using the flow rule of plasticity. The elastoplastic 
matrix is given by (26) 


[0,J.]= [G] 


1[C]: 


[ 0 ] 


®j_ i ® J, ® 

.6 o ) (_d o j 

LaaJ 15^1 


(2,35) 


where [c] is the elasticity matrix and H' is the slope of 
uniaxial stress versus plastic strain curve at a particular 
value of Cq. The elastoplastic matrix for Von Mises criterion 
and plane stress is obtained by simplifying the above equation 
as (51) 


^c 




' + 2P 


q/+2W) 

t t 

o’x + "^ 

(„ 

1+v 


(-hj,' + 2:i^P) 

^""x' x' "^^P) 


-a ' + pa ' 
_.X.. ^ ^ 


■) 


'^x +Vcr ' 

(_ 

' +V0y- ’ 

1 + 1^ xy 


(■ 


R 


2H' 


+ 


1+1> ^ xy 2(l+y) 


(l-i>) 


where 


P = 


2H' "a. 


9E, 


'xy 
^ + V 


(2.36) 


R = 




and 


!^ = [ R + 2( 1-/) P ] 

or ' , a' are the deviotric stresses i.e. 
X » y 
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“^x ~ ^x ” ^ 

CTx + + 0^2 

y ” °y " ' ‘3 

Most of the finite element investigators have used 
normality law of flow rule valid for yielding of concrete. 

The physical meaning of this flow rule, is that the increment 
of plastic strain has to be normal to yield surface. This 
rule has been shown to be applicable to most of the ductile 
metals. However, for concrete, the flow rule with respect to 
yield surface is not well established. The apparent plastic 
yielding of concrete is caused, not by actual plastic flow, 
but by cumulative effect of microcrack propagation (4-4) • 
Moreover, the yielding of concrete does not affect the overall 
behaviour of the structure to the extent the cracking of 
concrete, brickwork and yielding of steel does. Therefore, 
in the present work, an 'unconstrained flow rule ' suggested 
and used by lin et al (44) has been used. This rule assumes 
that the flow of plastic strains is not constrained but 
stresses are fixed at the initial yield points on the yield 
surface. The stress strain relationship for the unconstrained 
flow rule' in incremental form is expressed as follows 

idcy! = [0] [del 

where [ 0 J is a null matrix 


(2.57) 



2.5.3 Crusliins of Concrete 


A crush surface, analogous to yield surface but 
interms of strains, is postulated to define the complete 
collapse (or crush) for the yielded concrete (4-4? 53). fhe 
crush surface in terms of octahedral shear strain is expressed 
as 


5) = f2 


'■'oct 


cu 


0 


( 2 . 38 ) 


where is the ultiinate strain in uniaxial case and 


'oct' 


1 O O 0^99 

~[(e,-ey)h(ey-eph(.^-ep2 + | (Y,y + + yJ 1 


iNTow for plane stress condition Y„ = Y =0 and 

'yz 'zx z X y 


therefore , 


PcC^) 


If. 

3 


[ (8 -S )‘ 

U \ ^ y/ 



Z' 





— 6 

cu 


(2.39) 


The crush surface dLs shown in I'lg, 2,14 is Ihe 
principal strain plane. The boundary between crushed and 
noncrushed region defined by equation (2.39) is valid only for 
biaxial compression state . If one or both principal strains 
are positive i.e. tensile, no such boundary is defined. In 
such a condition, concrete will crack much before crushing 
strain is reached. Therefore, it is to_ sufficient to define 
bhe failure envelope only for the case of biaxial compression 
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Once concrete is crusliedjit is assumed to have lost all 
its stiffness. Thus, crashed concrete is assumed to have 
zero stiffness and therefore ,its material constitutive matrix 
is given hy 

:o 0 0*j 

[ o] = jo 0 0 ' (2.40) 

io 0 0 ’ 

2 . 5.4 Cracking of Concrete 


When one of the principal stresses exceeds the 
tensile strength of concrete » concrete is considered to have 
cracked in a direction perpendicular to that of principal 
stress. The normal stress at crack drops to zero and shear 
modulus gets reduced due to cracking. 


Referring to Rig. 2.15, the stress strain relation 
for cracked concrete in X'Y' coordinate is given by 



(2.41) 


where G- is the uncracked shear modulus, cc^ is a factor to 
account for aggregate interlock, dowel action etc. that may 
he present is always less than uni by y. X'T is the 
coordinate system having^^ coordinate axis coinciding 

with the principal stress directions at the tin^jOf 

QtM ^ 
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The above equation can be ■written as 

|cr'|= [O'] I e-J (2.42) 

oc^ the shear retention factor retains certain amount 
of shear stress in the cracked concrete. Since shear s fcrensth 
along the cracks is a function of crack width , the 
possible upper and lower bound values for are unity and 
zero , A value of zero will mean that the cracked element 
behaves as a b-undle of uniaxial fibres capable of sustaining 
tensile or compressive loads only parallel to the direction of 
crack. This is not true as cracks in concrete are not smooth, 
parallel and fnictionless slippage planes. Instead they are 
irregular rough planes at unequal distance apart. Based upon 
the numerical experimentation done by Jain A.K. (41) as value 
of 0,4 for cc has been used in the present 'work. 

The constitutive relations of equation ( 2 , 42 ) are 
transformed to global coordinate system using the following 
transformation. 

[«or] = P2J (2.43) 

where [O'] is the constitutive matrix of cracked element in 
local coordinate system X'Y' , [ 0 ^,^ ] is the constitutive 

matrix of cracked element in global coordinate system and 
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"Cos^^ Sin P Cos P 

^^2^ ~ Sin^iP Cos^iP -Sin p Gos P (2.44-) 

-2Sin;3CosP 2SinpCosp Gos^P-Sin^P j 

wliere P is the inclination of principal stress measured positive 
in anticlockwise direction from X axis in global coordinate 
system. 

If both principal stresses in an uncracked concrete 
element exceeds the tensile strength of concrete, concrete 
will crack in both principal stress directions. However , 
concrete cracked along one direction can further crack along 
a second direction when tensile stress in that direction 
exceeds the tensile strength of concrete. Like the first 
one jthis crack is also assumed to develops perpendicular 
to the direction of principal tensile stress. Ihe concrete 
cracked along two directions is assumed to be unable to 
transfer any load in tension. Therefore, concrete cracked 
along two directions is assumed to have zero stiffness. 

The constitutive matrix for such a case is given by 

"o 0 0~ 

t *^cr^ 0 0 0 (2.45) 

0 0 0 _ 

2.6 PROPERTIES OP STEEL REiMPORGBll/IEHT' 

Mild steel plain rounds have been as reinforcement in 
bending and as single legged Z shaped vertical connectors. 

¥ 

m 



Mild steel bars of all diameters used i.e. 6 to 16 mm were 
subjected to control tests. The yield stress, ultxmate stress 
and the p'jrcentago elongation were determined* Table 2*17 
gives the details of results obtained from control test* 

T’ig. 2,16 shows the stress strain curve obtained for the 
steel used. 

2.7 OOhSTlTUTI'TE EELATIOUS POR STEEL REIltFORGEMEOT 

The steel reinforcement has been idealised as an 
elastoplastic material in the present investigation with 
yield stress + and Modulus of Elasticity E„» The idealised 
stross-strain curve is shown in Eig. 2.17. It is further 
assumed that the reinforcing bars carry only axial stresses. 
When stresses in steel are in elastic range, the stress strain 
relation with respect to X'Y' coordinate system Fig. 2.18 
are 


Wx 


0 

0 j 

! 1 

1 1 

= 1 
i ! 

0 

0 

0 

i ®y't 

1 i 

It' 1 

V.^xyJ 

1 

0 

0 ^ 

1 i Y' 

■ 'v xy J 


Equation 2.44 can bo used to transform the above 
relation is global coordinate system. 

When reinforcing steel has yielded it is assimicd to 
have zero incremental s tiff ness. The constitutive relations 

in such cases are as follows. 

idcrj = [O] Uel 

where [O ] is a null matrix 


(2.47) 



EXPERHili^AL WORK 


SI. Wo, Length in cm 



^1 

L2 

^3 

h 

1= ZL;l /4 

1 

23.10 

23.30 

23.35 

23.10 

23.21 

2 

23.00 

22.90 

23.15 

22.95 

23.00 

3 

22.96 

22.70 

22.75 

23.00 

22.85 

4 

22.90 

23.10 

22.95 

23.05 

23.00 

5 

22.85 

23.50 

23.00 

22.85 

23.05 

6 

23.35 

23.30 

23.35 

23.35 

23.34 

7 

23.35 

22.75 

22.80 

23.40 

23.08 

8 

25.30 

23.10 

23.00 

23.35 

23.19 

9 

22.85 

23.10 

23.10 

22.95 

22.96 

10 

23.10 

23.50 

23.40 

23.15 

23.29 

1 1 

22.75 

22.40 

22.50 

22.70 

22.59 

12. 

23.30 

23.15 

23.50 

23.35 

23.33 

15 

22.95 

23.15 

23.20 

22.95 

23.05 

14 

23.21 

23.30 

23.35 

23.30 

23.29 

15 

22.60 

22.65 

22.65 

22.70 

22.65 

16 

23.20 

23.30 

23.30 

23.25 

23.26 

17 

22.85 

23.05 

23.05 

22.95 

22.98 

18 

23.10 

22.95 

23.00 

23.05 

23.03 

19 

23.30 

23.10 

23.05 

23.25 

23.16 

20 

23.20 

23.15 

23.15 

23.15 

23.16 

21 

23.25 

23.35 

23.40 

23.40 

23.35 

22 

23.15 

23.25 

23.20 

23.30 

23.23 

23 

22.95 

23.05 

23,00 

22.95 

22.99 

24 

22.90 

23.00 

22.95 

22.85 

22.93 

25 

22.95 

23.10 

23.20 

23.15 

22.85 

26 

23. 10 

23.15 

23.05 

23.10 

23.11 

27 

23.15 

23.05 

23.20 

23.10 

23.24 

28 

22.90 

23.10 

23.10 

22.95 

23.01 

29 

23.30 

23.20 

23.35 

23.40 

23.31 

30 

22.80 

22.75 

22.65 

22.70 

22.73 


Arithmetic mean = 


Sum 

*30 


= I 


Standard deviation - 


= _ 23 .^_ 
' n-1 


Coefficient of variation’ 
=0,88 percent 

= 0,203 

Table contd, .on page 56 


dUJ-t: I con'td. . . 


Breadth in cm 



^1 


% 


B=EB3.A 

1 

11.10 

11.15 

11.50 

11.50 

11.31 

2 

11.05 

11.10 

10.95 

11.00 

11.14 

3 

10,85 

10.75 

10.75 

10.85 

10.85 

4 

10.95 

11.00 

10.90 

11.05 

11.09 

5 

10.95 

11.15 

11.15 

11.15 

11.10 

6 

1 1.25 

11.20 

11.15 

11.25 

11.21 

7 

10.70 

11.10 

11.50 

10.80 

11.03 

8 

10.90 

11.25 

11.40 

10.95 

11.13 

9 

10.80 

10.95 

11.50 

10.90 

11.04 

10 

10,85 

10.85 

10.85 

10.95 

10.88 

1 1 

11.00 

11.15 

11.15 

11.15 

11.11 

12 

10.65 

10.95 

11.50 

10.70 

10.95 

13 

11.20 

11.20 

11.15 

11.15 

11 .18 

14 

11.35 

11.50 

11.55 

11.30 

11.43 

15 

10.75 

11.00 

10.80 

11.00 

10.89 

16 

10.65 

10.95 

10.95 

10.70 

10.81 

17 

11.05 

11.30 

11.30 

11.05 

1 1.18 

18 

10.95 

10.70 

10.75 

10,85 

10.81 

19 

10.90 

10.90 

10.95 

11.00 

10.94 

20 

10.95 

10.85 

10.85 

10.90 

10.89 

21 

11,10 

11.05 

11.15 

1 1 .00 

11.08 

22 

11.05 

10.85 

10.95 

11.00 

10.96 

23 

11.25 

11.30 

11.15 

11.20 

11.25 

24 

11.00 

11.15 

11.25 

11,15 

11.13 

25 

11.10 

10.95 

10.95 

1 1 .00 

11.00 

26 

11.15 

11.35 

11.35 

11.25 

11.28 

27 

10.75 

11.00 

10.95 

10.90 

10.90 

28 

10.85 

10.95 

11.20 

10,95 

10.99 

29 

10.65 

10.95 

11.00 

10.95 

10.89 

30 

10.80 

11.00 

10.95 

10.90 

10.91 


Arithaetic mean - 
Standard, deviation 


= B = 11.045 Coefficient variation 
30 =1»4-5 percent 

a =f E = 0.158 

^ Table contd, .on page 
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Table 2.1 contd... 


SI. No. 


Depth 

in cm 




^1 

D2 


“4 

D=2Di/4 

1 

6.65 

6.25 

6.10 

6.50 

6.38 

2 

6.60 

6.40 

6.35 

6.25 

6.40 

3 

6.50 

6.35 

6.65 

6.45 

6.49 

4 

6.40 

6.55 

6.45 

6,40 

6.45 

5 

6.10 

6.35 

6.45 

6.25 

6.29 

6 

6.55 

6.64 

6.45 

6.54 

6.54 

7 

6.45 

6.60 

6.50 

6.70 

6.56 

8 

6.50 

6.45 

6.35 

6.30 

6.40 

9 

6.75 

6.50 

6.55 

6.70 

6.63 

10 

6.35 

6.25 

6.50 

6.41 

6.38 

1 1 

6.30 

6.55 

6.30 

6.20 

6.34 

12 

6.65 

6.30 

6.10 

6.30 

6.34 

13 

6.65 

6.70 

6,65 

6.55 

6.64 

14 

6.60 

6.45 

6.50 

6.30 

6.46 

15 

6.60 

6.65 

6.60 

6.25 

6.53 

16 

6.65 

6.40 

6.45 

6 .60 

6.53 

17 

6.75 

6.70 

6.35 

6.60 

6.60 

18 

6.75 

6.35 

6 . 40 

6.50 

6.50 

19 

6.55 

6.50 

6.45 

6.60 

6.53 

20 

6.55 

6.60 

6.60 

6.65 

6.60 

21 

6.60 

6.45 

6.50 

6.50 

6.51 

22 

6.55 

6.45 

6.40 

6.60 

6*50 

23 

6.45 

6.60 

6.55 

6.50 

6.53 

24 

6.35 

6.65 

6.50 

6.40 

6,48 

25 

6.30 

6.55 

6.35 

6.45 

6.41 

26 

6.45 

6.55 

6.60 

6.50 

6.53 

27 

6.60 

6.45 

6.50 

6.35 

6 . 48 

28 

6.70 

6.55 

6.60 

6.70 

6.64 

29 

6.65 

6.45 

6.20 

6.50 

6.45 

30 

6.55 

6.60 

6.65 

6. 45 

6.55 


mean 


- Sam 
50 


Arithmetic 
Standard deviation 


= V 


C = 6.42 CO€ 

r = 0.092 

n~1 


_J.LiJLCJ.XU vcox- 

=1.43 percent 



PERCEOTAGil WilJTER JIBSORBTIOR OF BRICKS USED IN 
EXPERDIEUTiUD WORK 


SI. 

No. 

Area 

AadLxB 

Volume 

Y-lixBxD 

Dry 
weight 
in kg 

V/et 

we ight 
in kg 

Dry 

density 

in 

Percentage 

absorbtion 

1 

262.51 

1674.78 

2, 646 

2.863 

1.580 

8.20 

2 

256.22 

1639.81 

2.615 

2.860 

1.595 

9.37 

3 

247.92 

1609.02 

2.534 

2.859 

1 .575 

12.83 

4 

255.07 

1645.20 

2.545 

2.880 

1.547 

13.16 

5 

255.86 

1609.33 

2.615 

2,910 

1.625 

11.28 

6 

261.64 

1711 .13 

2.498 

2.821 

1.460 

12.93 

7 

254.37 

1669.98 

2.687 

3.082 

1.609 

14.70 

8 

258.10 

1651.84 

2.631 

2.968 

1.593 

12.81 

9 

253.48 

1680.57 

2.544 

2.887 

1.514 

13.48 

10 

253.40 

1616.69 

2.519 

2.908 

1.558 

15,44 

1 1 

250.97 

1591 .15 

2.531 

2.848 

1.591 

12.52 

12 

255.46 

1619.62 

2.554 

2.866 

1.577 

12.22 

13 

257.70 

171 1.13 

2.722 

2.842 

1 .591 

4.41 

1 4 

266.20 

1719.65 

2.6O7 

3.109 

1.516 

19.26 

15 

246.66 

1610.68 

2,682 

2.976 

1 . 665 

10.63 

16 

251.44 

1641.90 

2.772 

3.021 

1.688 

8.98 

17 

256.92 

1695.67 

2.651 

3.034 

1.563 

14.44 

18 

248.95 

1618.18 

2.633 

3.056 

1.627 

16.07 

19 

253.37 

1654.57 

2.447 

2.903 

1.479 

18. 64 

20 

252,21 

1664.59 

2.691 

2.773 

1.617 

3.05 

21 

258.72 

1684.27 

2.770 

3.015 

1.645 

8.84 

22 

254.60 

1654.90 

2.765 

3.120 

1.671 

12'. 84 

23 

257 . 41 

1 681 . 41 

2.575 

2.981 

1.531 

15.77 

24 

255.21 

1653.77 

2.642 

3,010 

1.598 

1 3.9'3 

25 

251.35 

1611.15 

2-. 775 

3.100 

1.722 

11.71 

26 

260.68 

1702.24 

2.615 

2.910 

1.536 

11.28 

27 

253.32 

1641 .51 

2.687 

2.990 

1.637 

11.28 

28 

252.88 

1679.12 

2.631 

2.972 

1.567 

12.96 

29 

253.85 

1637.33 

2. 605 

3.109 

1.591 

19.35 

30 

247.98 

1631 .74 

2.770 

3.025 

1.698 

9.20 

A.M 

254.830 

165'3.762 



1.592 

12.386 

S.D 

4.446 

35.038 



0.062 

3.728 

c^ 

1'.740 

2,119 




xn -1 r\o 



C;U 1 VLPRJJ]SSI 7 E 

STBEiIGTH OF BRICKS 


Sl.iTo. 

Area 

Failure load 
in Fg. 

Compressive 
strength _ 
in kg/cm^ 

1 

262.51 

77750 

296.179 

2 

256.22 

72000 

281.008 

3 

247.92 

70000 

282.370 

4 

255.07 

71050 

278.551 

5 

255.86 

74000 

289.221 

6 

• 

— 

(M 

64750 

246 . 480 

7 

254.37 

57250 

225.066 

8 

258.10 

77000 

298.334 

9 

253.48 

66500 

26 2.348 

10 

253.40 

54750 

216.060 

1 1 

250.97 

70250 

279.914 

12 

255 . 46 

77000 

301.417 

13 

257.70 

8O5OO 

312.379 

14 

266.20 

58500 

219.760 

15 

246.66 

56500 

229.116 

16 

251.44 

72500 

288,339 

17 

256.92 

72000 

280.265 

18 

248.95 

63250 

254.067 

19 

253.37 

83500 

329.558 

20 

252.21 

83000 

329.090 

21 

258.72 

64500 

249.309 

22 

254.60 

68000 

267.086 

23 

257.41 

71500 

277.770 

24 

255.21 

69500 

272.325 

25 

251.35 

59500 

236.722 

26 

260,68 

63250 

242.655 

27 

253.32 

70250 

277.317 

28 

252.88 

57750 

228,569 

29 

253.85 

62500 

246.208 

30 

247,98 

78500 

316.558 

A.M. 


a 

270.550 

S.D 



31.778 

Gv 



11.745 



TABLE 2 4 ; COl'IPRESSIOK TiJT ON 7 CM MORTM CUBES 
PROPORTION 1:3 BY lEI&HT 


SI. No. 

Jo. of Days 
of Curing 

P allure 
Load 
in Kg. 

Compressive 
Sfress in 
lcg/cm2 

1 

29 

12900 

25a 

2 

29 

12100 

242 

3 

29 

11500 

230 

4- 

28 

12200 

244 

5 

28 

11100 

222 

6 

28 

12700 

254 

7 

30 

10800 

216 

8 

30 

10200 

204 

9 

30 

11100 

222 

10 

28 

11300 

226 

1 1 

28 

13500 

270 

12 

28 

12000 

240 

13 

29 

10200 

204 

14 

29 

10400 

208 

15 

29 

11000 

220 

16 

30 

1 1500 

230 

17 

30 

11400 

228 

18 

30 

12100 

242 

19 

30 

12100 

242 

20 

30 

10800 

216 

21 

30 

11300 

226 

22 

27 

12100 

242 

23 

27 

10700 

214 

24 

27 

1 1000 

220 

25 

28 

11300 

226 

26 

28 

10200 

204 

27 

28 

10800 

216 


Mean 


228.37 



I 


TiLBLii 2,5: COWJCPRESSIOiNF TE&J OE 7 CM MORTAR CUBES 
PROPO-^TIOU 1 :4 BY YCEIG-HT 


SI . Uo . 

Uj. of days 
of curing 

Failure 
load in 
kg. 

Compressive 
strength, in 

p 

kg/ cm 

1 

29 

10100 

202.0 

2 

29 

9100 

182.0 

3 

29 

9000 

180.0 

4- 

30 

9500 

190.0 

5 

30 

9200 

184.0 

6 

30 

8700 

174.0 

7 

28 

8500 

170.0 

8 

28 

8800 

176.0 

9 

28 

9800 

196.0 


Me^ ii 


183.77 



62 


TiBiB 2.6 : GOlvIPBBSSIOi^ TEST ON 7 CM MORTiB CUBES 
BROBuRTION 1:5 BI \BIGHT 


SI. No, 

No . of days 
of caring 

E allure load 
in Kg. 

Compressive 
strength in 

kg/ cm^ 

1 

28 

5700 

74.0 

2 

28 

3200 

64*0 

3 

28 

3500 

70.0 

4 

30 

4000 

80.0 

5 

30 

3100 

62.0 

6 

30 

5400 

68.0 

7 

31 

3000 

60.0 

8 

31 

3100 

62.0 

9 

31 

2900 

58.0 


Mean 


66,444 



TABLJ 2.7 : 


COMPRESSION TEST ON 7 CM CUBES 
PROPORTION 1 :6 BY WEIOHT 


63 


SI . No . 

No . of days 
of curing 

P allure load 
kg. 

Compressive 
stress in 

kg/ cm2 

1 

29 

2800 

56.0 

2 

29 

2700 

54.0 

3 

29 

2400 

*48.0 

4 

29 

2900 

58.0 

5 

29 

2750 

55.0 

6 

29 

2800 

56.0 

7 

28 

2500 

50.0 

8 

28 

2200 

44.0 

9 

28 

2800 

56.0 

10 

30 

2900 

58.0 

1 1 

30 

3100 

62.0 

1 2 

50 

2800 

56.0 

13 

31 

2600 

52.0 

14 

31 

2800 

56.0 

15 

31 

3200 

64.0 

16 

29 

2750 

55.0 

17 

28 

2400 

48.0 

18 

28 

2600 

52.0 

19 

30 

2500 

50.0 

20 

30 

2700 

54.0 

21 

30 

2900 

58.0 

22 

30 

3200 

64. 0 

23 

30 

3000 

60.0 

24 

30 

2600 

52.0 


Mean 


54.916 
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TA3L2 2 8 : G014PEBSSI0F l^ST 02 7 CM MORTM CUBE 

PR0P0R1'I02 1;8 BY Y®IGHT 


.No . 

No. of days 
of curing 

failure load 
in kg. 

Compressive 
strength, 
in kg/cm2 

1 

29 

1250 

25.0 

2 

29 

HOO 

28.0 

3 

29 

1100 

22.0 

4 

29 

1000 

20.0 

5 

29 

1150 

23.0 

6 

29 

1300 

26,0 

7 

30 

1200 

24.0 

8 

30 

1 150 

23.0 

9 

30 

1050 

21.0 


Mean 


23.555 



TABLiL 2.9; CAlCUllTION OS' STRSSS il© STRAIU FOR ONE PRISM 
SPEC I EN OP TYPE H IN 1:3 MORTAR 


] 

« i 

• 

H 

CO 

Load in tonnes 

o 

f1 

H 

CQ 

CQ 

CD 

: 

CO 1 

j Deflection of 
i rubber sheet 
! corresponding to 

1 stress in col, ( 3 ,) 

1 in mm 

Reading of dial 
gauge for test 

1 assembly in nrni 

Deflection of j 

|test assembly | 

I in mm f 

Net deflection 
of prism in mm 
(6)-(4) 

t 

1 

; 

[ Strains x 10“'^ 
corresponding to 

stress in Col. (31 

1 

Assumed strains 

xlO-4 

Interpolated 
stress in kg/cm^ 
for strain in( 9 )' 

( 1 ) j 

1 ( 2 ) ! 

(3) 

(4) 

|(5) ! 

(6) 



(9) j ( 10 ) 


1 -1 



i , - -.| 

I : 




1 

0.00 

0.000 

0.00 

6.39 

0.00 

0.00 

0.00 

1 

1.289 

2 

1.00 

0.794 

0.10 

6.58 

0.19 

0.09 

0.63 

2 

2.629 

3 

5.00 

3.968 

0.50 

7.32 

0.93 

0.43 

3.00 

3 

3.968 

4 

10.00 

7.936 

0.97 

8,14 

1).75 

0.78 

5.45 

4 

5.588 

5 

15.00 

11.905 

1 .30 

8.72 

2.33 

1.03 

7.20 

5 

7.208 

6 

20.00 

15.873 

1.56 

9.13 

2.74 

1.18 

8.25 

6 

9.183 

7 

25.00 

19.841 

1 .75 

9.47 

3.08 

1.33 

9.30 

7 

11.452 

8 

27.50 

21.825 

1.82 

9.63 

3.24 

1.42 

9.94 

8 

14.929 

9 

30.00 

23.809 

1.92 

9.79 

3.40 

1.48 

10.35 

9 

18.706 

10 

32.50 

25.994 

2.00 

9.97 

3.58 

1.58 

11.05 

10 

22.116 

11 

35.00 

27.777 

2.07 10.30 

3.91 

1.84 

12.87 

11 

25.651 

12 

37.50 

29.762 

2.14 10.62 

4.23 

2.09 

14.60 

12 

26.829 

13 

40.00 

31.746 

2.20 11.03 

4.64 

2.44 

17.05 

14 

29.074 

14 

42.50 

33.730 

2.26 11.22 

4.83 

2.57 

18.00 

1 6 

30.810 

15 

45.00 

35.714 

2.33 

- 


- 

- 

18 

33.730 

16 

47.50 

37.698 

2,39 

- 

“ 

- 

— 

20 

33.730 

17 

50,00 

39,683 

2.45 


- 

- 

- 

25 

33.730 
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TABLji) 2.10 

: COICPRESSIVE 

STBEMH 

OF BRICK 

PRISM 

TYPE H 


SI. 

Eo, 

Mortar Gompressive Strength in kg/cm^ 


, Mean 




Sample 

No. 



value 



1 

2 

3 

4 

5 

6 


1 

1 j3 

35.714 

42.857 

41.667 

42.857 

37.698 


40.159 

2 

1 :4 

27.778 

32.738 

31.746 

30.357 

31.746 

34.524 

31 .482 

5 

1 :5 

25.794 

27.778 

29.762 

29.762 

27.778 

30.952 

28.638 

4 

1 :6 

26.825 

21.429 

26.587 

26.190 

26.984 


25.603 

5 

1 i8 

20.238 

27.222 

16.508 

22,222 



21 .548 

TiLBlIil 2.11 

: G0i,viPRSS3lVE 

STRENGTH 

OF BRICK PRISM 

TYPE V 


SI. 

Eo, 

Mortar < 

2 

Compressive Strength in k^/cm 


Mean 

Value 

MIX 



Sample 

No. 





1 

2 

3 

4 

5 

6 


1 

1 :3 

31.746 

33.730 

29.762 

42.460 

35.714 

39.683 

35.516 

2 

1 :4 

22.222 

29.762 

29.762 

30.555 

21 .826 

32.540 

27.778 

3 

1 °.5 

14.286 

14.683 

17.659 

16.071 

17.063 

18.650 

16.402 

4 

1 :6 

13.492 

15.159 

18.651 

10.714 

16.270 

* 

14.857 

5 

1 :8 

6.825 

8.730 

9.127 



# 

8.227 


^ Broken during liandling i.e. ■transport; ing from curing sife 
to the Universal Testing Machine* 



T^lBLE 2.12: MODULUS OF EL.iSl’IOITY OE BRICK PRISMS TYPE H 


Spec me r 
number 

Secent Modulus 
percentage of 

in kg/ cm^ X 1 0 at 
mean compressive stress 



10 

20 

40 

60 

80 

Mortar mix 

1 

i :3 

1.323 

1.457 

1 .924 

2.300 

1 .862 

2 

1 .418 

1.536 

2.736 

3,662 

4.372 

3 

1 .418 

1.383 

1.783 

2.267 

2.685 

4 

1 .418 

1.554 

2.248 

2.861 

3.488 

5 

1.383 

1 .367 

1 .907 

2.467 

2.701 

Mortar mix 






1 :4 






1 

0.908 

1 .007 

1.645 

1 .855 

1.543 

2 

1.091 

1.143 

1.756 

1,654 

1.482 

3 

1.342 

1.408 

1.855 

2.006 

1.492 

4 

1 .1 24 

1,176 

1,657 

1.757 

1.613 

5 

1.115 

1.134 

1.543 

1.578 

1.575 

6 

0.966 

0.990 

1 .478 

1 .514 

1.249 

Mortar mix 

1 :5 






1 

0.823 

1,042 

1.514 

1.342 

1.167 

2 

0.912 

0.934 

1.146 

1.267 

1.085 

3 

0.857 

0.914 

1.312 

1 .214 

0.986 

4 

0.842 

0.872 

1.356 

1.425 

1.384 

5 

0.861 

0.965 

1.524 

1 .617 

0.7 27 

6 

0.868 

0.811 

0.914 

0.695 



Strain measuring dial gauge was removed before reaching 
80 percent average stress. 

Table contd.,,0'n page 68 


Table 2,12 contd... 


Specimen Secant modulus in kg/cm^ x 10^ at percentage 

number of mean compressive stress 



10 

20 

40 

60 

80 

Mortatr mix 1 :6 

1 0.798 

0.806 

0.902 

0.938 

0.855 

2 

0.699 

0.759 

1 .210 

1.493 

0.755 

3 

0.521 

0.600 

0.702 

0.632 


4 

0.846 

0.852 

1 .027 

1 .229 

1.246 

5 

0.747 

0.772 

1 .037 

1 .1 12 

0.845 

Mortar mix 

1 :8 






1 

0.510 

0.523 

0.732 

0.625 

*Tf 

2 

0.5 16 

0.571 

0.766 

0.815 

0.727 

3 

0.726 

0.731 

0.841 

0.855 

0.820 

4 

0.732 

0.782 

0.829 

0.873 

0.805 


Strain measuring dial gauge was removed before reaching 
80 percent average stress. 



TABLE 2.15*. modulus QP ELASTICITY OP BEICK PRISMS TYPE V 


Specimen Secent modulus in kg/cm^ x lO'^ percentage 

number mean compressive stress 



10 

20 

40 

60 

80 

Mortar mix 

Til 





1 

1 .526 

1.373 

1.575 

1.580 

1.197 

2 

1.628 

2.377 

3.691 

3.674 

2.081 

5 

3.055 

3.060 

3.752 

3.065 


4- 

1.952 

3.153 

4.372 

4.694 

2.621 

5 

3.296 

3.8t7 

5.357 

3.763 

2.621 

6 

4.278 

3.693 

3.246 

3.829 

4.085 

Mortar mix 





1 H 






1 

1 .914 

1 .981 

1.869 

1 .406 

** 

2 

0.956 

1.236 

1 .674 

1 .435 

1.112 

3 

0.878 

0.980 

1 .295 

1 .314 

1.005 

4 

0.765 

1 .036 

1 .411 

1.173 

1.04-4 

5 

2.816 

2.724 

2.549 

1 .098 

** 

6 

1.627 

1 .806 

2.201 

2.153 

2.101 

Mortar mix 





1 :5 






1 

1.713 

1.752 

1 .801 

1 .811 

1.924 

2 

1.514 

1.601 

1 .623 

1.610 

1 .426 

3 

1.362 

1.371 

1 .412 

1.394 

1.012 

4 

1.028 

1 .029 

1 .124 

1.216 

1.314 

5 

1.276 

1 .314 

1 .213 

1.189 

0.784 


1.107 

1 .100 

1.217 

1 .312 



★^Strain measuring dial gauge was removed before reaching 
80 percent average stress. Table contd..on page 70 



Table 2,13 contd,.. 


Specimen Secant modulus in kg/cm^ x at percentage 

number of mean compressive stress 



10 

20 

40 

60 

80 

Morten mix iVS" 

1 

0.558 

0.662 

0.643 

0.675 

0.637 

2 

1.067 

1.277 

1.237 

0.977 

0*882 

3 

1.408 

1.374 

1.640 

1 .800 

1 .984 

4 

1.4S7 

0.575 

0.401 

0.339 


5 

1.073 

1.050 

1.317 

1.858 

1 *846 

Mortar mix 






1 :8 






1 

0.269 

0.268 

0.325 

0.293 

TTW 

2 

0.260 

0.260 

0.306 

0.263 

0.233 

3 

0.265 

0.265 

0.301 

0.263 

0.216 


Strain measuring dial gauge was removed before reaching 
80 percent average stress. 




TiiJLE 2 . 14 : MilXIMUM RESIST lICE TO DiAGOEillj TENSION IN 


REI FORCED brick BB4MS 


Specimen 

Number 

Maximum total 
load in kg. 

Maxiium deflection 
at mid point in 
mm 

Maximum shear 
stress as 
diagonal tension 

kg/cm^ 


Mortar Mix 1:3 

Modular ratio 

80 


1 

7292.516 

1 .81 


5.567 

2 

8625.850 

1 .70 


6.584 

5 

9306,122 

1.98 


7.104 

Mean 

8408,163 

I 

00 

• I 

^ 1 

1 

1 

j 


6.418 


Mortar mix 1 ;4 

Modular ratio 

100 


1 

6893.424 

2.10 


5.261 

2 

8163.265 

1 .70 


6.230 

5 

8163.265 

1 .78 


6.230 

Mean 

7739.985 

1.86 


5.907 


Mortar mix 1 :5 

Modular ratio 

160 


1 

4136,054 

0.910 


3.145 

2 

4952,381 

1.290 


3.766 

3 

3918.367 

1 .050 


2.980 

Mean 

4335.601 

1 .083 


3.297 


Mortar mix 1 :6 

Modular ratio 

200 


1 

2802.721 

1.340 


2.124 

2 

3401.360 

1.200 


2.577 

3 

2993.197 

1 .460 


2.268 

Me an 

3065.759 

1.333 


2,323 


Mortar mix 1 :8 

Modular ratio 

20 b 


1 

2476.190 

1.65 


1.876 

2 

2721 ,088 

1.60 


2.062 

2 

2938.775 

1.70 


2,227 

Mean 

2712.018 

1.65 


2.055 i 


i- 


T:IBLE 2 . 15 : COMPBESSIOU TEST OM 15OM CONCRETE CUBE 


SI. No. 

Number of 

3 - ays of 
curing 

E allure 
load 
in kg. 

Compressive 
strength 
in kg/cm^ 

1 

30 

46500 

206.667 

2 

30 

45500 

202.222 

3 

30 

48500 

215.555 

4 

28 

47000 

208.889 

5 

28 

47500 

21 1 ,11 1 

6 

28 

47000 

208,889 

7 

28 

47500 

2 t 1 .111 

8 

28 

45500 

202.222 

9 

28 

46500 

206,667 

10 

27 

46000 

204.440 

11 

27 

48000 

215.330 

12 

27 

47500 

21 1 .1 11 

13 

29 

47500 

21 1.1 11 

14 

29 

45000 

200.000 

15 

29 

45500 

202.222 

16 

29 

46 000 

204.444 

17 

29 

47000 

208.889 

18 

29 

45500 

202.222 

19 

28 

47000 

208,889 

20 

28 

46500 

206.667 

21 

28 

44500 

197.778 

22 

30 

47500 

21 1 .1 11 

23 

30 

45000 

200.000 

24 

30 

45500 

202.222 

25 

27 

44500 

197.778 

26 

27 

47000 

208.889 

27 

27 

47500 

211.111 

28 

29 

45000 

200.000 

29 

29 

44000 

195.556 

30 

29 

44500 

197.778 


Table contd 


on page 



Table 2,15 contd... 


Sl,iTo. 

Nombcr of 
da^^s of 
cui-'ing 

Failure 
load 
in kg. 

Compressive 
strengtlip 
in kg/cm'^ 

31 

29 

46000 

204 , 444 

32 

29 

45000 

200.000 

33 

29 

47000 

208.889 

34 

28 

44500 

197-778 

35 

28 

46000 

204,444 

36 

28 

47000 

208.889 

37 

30 

46500 

206.667 

38 

30 

47500 

211,111 

39 

30 

46000 

204,444 

40 

28 

47000 

208.889 

41 

28 

45500 

202.222 

42 

28 

47500 

21 1,111 

43 

28 

46000 

204.444 

44 

28 

46500 

206.667 

45 

28 

45000 

200.000 

46 

29 

47500 

21 1. Ill 

47 

29 

48000 

213.333 

48 

29 

46000 

204.444 

49 

2^ 

45500 

202.222 

50 

27 

46000 

204.444 

51 

27 

46500 

206.667 


Mean = 205,621 kg/ cm 



lA 


2.16: 

split tlst 

ON C TCROTN 

CYLINDERS POR TENSILE 


STRSJG-TE OP 

GONCEETE 


8ize 

of cylinder 

= 15 cm dia 

and 50 cm long 

SI. No. 

No .of 
days of 
c ur ing 

Failure 
load P 
in kg. 

Tensile 2P 

strength = — — 
in kg/cm^ ixdl 

1 

28 

18000 

25.465 

2 

28 

19500 

27.587 

3 

28 

17500 

24.757 

4 

28 

18500 

26.172 

5 

28 

17000 

24.050 

6 

28 

18500. 

26.172 

7 

28 

17500 

24.757 

8 

28 

18000 

25.465 

9 

28 

19000 

26.880 

10 

28 

19000 

26.880 

11 

28 

17000 

24.050 

12 

28 

17500 

24.757 


Mean 


25.582 
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TABLE 2.17^ TENSILE STRENcrTH TEST RESULTS OF REIl®'ORCEMEIT 


SI. 

No . 

Nomi- 
nal 
dia- 
meter 
of bar 
in mm 

Actual 
mean 
dia of 
bar in 
mm 

Actual 

area 

0 Z bar 
in mm^ 

Yield 

load 

in 

kg. 

Yield 

stress 

in 2 

kg/mm'^ 

Ulti- 

mate 

load 

in 

kg. 

Ulti- 

mate 

tensi- 

le 

stre- 
ngth. 2 

in kg/mm 

Percen- 

tage 

Elon- 

gation 

1 

6 

5*94 

27.712 

800 

28.87 

1200 

43.300 

36.7 

2 

6 

5.91 

27.432 

840 

30.62 

1210 

44.1 10 

35.4 

3 

6 

6.00 

28.274 

830 

29.35 

1210 

42.795 

.22.9 

Me an 





29.615 


43.40 

31 .66 

4 

8 

8.10 

51 .‘530 

1550 

50.080 

2415 

46.860 

19.5 

5 

8 

8.08 

51.276 

1650 

32.179 

2480 

48.366 

24.6 

6 

8 

8.08 

51.276 

1550 

50.229 

2470 

48.170 

22.3 

Mean 





30.829 


48.465 

22.13 

7 

10 

9.26 

67.346 

2050 

30.440 

2950 

43.80 

29.4 

8 

10 

9.64 

72.987 

2150 

29.457 

3150 

43.16 

31 .3 


10 

9.60 

72.382 

2200 

30.392 

3100 

42.83 

32.4 

Mean 





30.096 


43; 262 

31 F3 

10 

TF" 

12.12 

"F'sTFr 

3400 

29.470 

5150 

44.64 

28.5 

1 1 

12 

12.14 

115.752 

3300 

28.509 

5300 

45;79 

27.3 

12 

12 

12.10 

114*990 

3150 

27.393 

5400 

46 .'96 

29.8 

Me an 





28.457 


45.796 

28.53 

13 

16 

15.84 

197.061 

5750 

29.179 

9500 

48.21 

24 .'6 

14 

16 

15.92 

199.056 

5800 

29.137 

9200 

46; 22 

25^3 

15 

16 

15.92 

199.056 

5700 

28.635 

9400 

47.22 

27.5 

Mean 





29.187 


47.22 

25.8 





Mean 

29.637 

Mean 

45.628 







Type H prism 


Type V prism 


Fig. 21 Arrangement of brick layers in brick prisms 



Fig.2*2 Stress-strai 
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Fig. 26 Idealised stress-strain curves of brickwork 
for type V prisms 
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curves for RBM beams 








Fig, 2-13 Biaxial strength envelope of concrete 
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Fig. 2’16 Stress-strain curve for steel 


Fi 92-17 Idealised stress-strain curve for steel 


Reinforcing 
steel bar 



Fig.218 Local and Global coordinates for framing 
constitutive matrix for reinforcing steel 



CH.iPi'fiR III 


fiXPDR 'Jj IRVjJSIlLfAlIOi^ 

iSarlier inYestigators have reported that no 
interaction takes place between the brick masonry supported 
on fi.G. beams if the heiaht to span ratio is below 0.5 and 
the inplane load is applied on the top of brickwork, further- 
more, no interaction is reported for such composite construction, 
if the inplane load is applied at the junction of brickwork 
and A.C, beam for any hei 2 ht to span ratio, whatsoever. 

Yifood (7 ) has suggested for the latter case the use of 
suitable tensile connectors for the interaction to come into 
play. Moreover, the general recommendation is that the 
benefit of such interaction should only be availed, if bhe 
brickwork is cast in rich mortars, viz. 1:3. ike experimental 
investigation of interaction between the brick masonry and 
R.C, beam supporting it, 'using single legged Z shaped 
connectors, is the subject matter of study in the present 
chapter. In other words, composite behaviour of such a 
construction is experimentally observed for two types of 
inplane loading on the system mentioned above and for cases 
(l) of different mortar strength by proportioning its cement 
sand ratio, (2) of varying height to span ratio of the 
structure and ( 3 ) of different size and location of openings 
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in the brickv/orlc which is reported herein. 

3.1 PREP.iIlATIOh OF iESi SPSCIlv£&;h 
3»1.1 ReinTorced Concrete Beams 

The reinforced concrete beams have been cast in 
wooden moulds prepared for the purpose. Pigs. 3.1(a) and 
3.1(b) show respectively the details of wooden moulds used 
for casting the beams for (i) compressive loading (11) 
tensile loading. Phe size of R.G. beams is kept as 325 cm 
long, 24 cm wide and 8 cm thick in most cases. Changes in 
thickness of the beam, if any, are discussed for particular 
cases. Bending reinforcement in beams consisted of 3 
plain round mild steel bars. Diameter of bars varies of 
course based on height to span ratio. Single legged Z shaped 
vertical connectors, thirteen in number, were embedded in 
all concrete beams along their length at a uniform spacing 
of 25 cm centre to centre according to the modular size of 
the bricks used in the present work . Limited experimentation 
has been carried out by varying the size and spacing of 
the vertical connectors. Beams used for tensile loading 
had mild steel flats embedded in it all along the length 
of beam, at a uniform spacing of 25 cm centre to centre, to 
facilita.te the loading arrangement. The mild steel flats 
were 35 cm long 8 cm wide and 1 cm thick. They were 
embedded gust at the top of concrete. Each mild steel flat 



92 


had tw'" 25 mm dia holes at their ends, 30 cm centre to 
centre, '.Thus 12 such flats were embedded in each bean, 
figures 3»2(ajand 3*2(b) show respectively the beams used 
for compressive and tensile loading. The beams were cured 
by putting the I'loist empty gunny cement bags over it for 
a period of 4 to 6 days before removing them from moulds. 
However, the sides of moulds were removed after 24 hours 
only. :ifter removing the beams from mould, they were kept 
supported at four points and curing was continued till 
the construction of brickwork over them. 

3.1.2 Brick Masonry Y/alls 

Brick masonry walls, with height to span ratio less 
than 0.4, were cast is casting yard. Those, with height to 
span ratio of .greater than 0.4, were cast after placing the 
precast concrete beams over s uitable supports (two end 
supports on rollers and two temporary supports at one 
third span) right on the lest floor. The temporary supports 
were necessitated due to the fact that the strength of 
beams after one week of curing was just sufficient to 
support the self weight. Any appreciable over burden would 
have caused excessive deflections and consequent failure 
of be^ma at the time of specimen preparation. The brickwork 
was done by two experienced masons. In all cases, brickwork 
was completed in one day. The bricks were soaked well in 
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water telore using. I'he r^bio of cement and sand in the 
mortar was token ’:iy weight. The control tests on mortar 
and bricks Yrere carried out and have already been reported 
in Chapter II. Gai^e was exercised to maintain the constant 
mortar thickness. It turned out to be 1 cm for horizontal 
mortar oeds and 1.5 cm for vertical joints. Brick masonry 
was cured for about 28 days by sprinkling water 2 to 3 times 
a day over the walls, 

3.1.3 loading 

A uniformly dis'cributed load was simulated on the 
test specimen which was ensured to be horizontal, by applying 
the load at twelve points. The centre to c entre distance 
between the load points was 25 cm. Twin box type of test 
floor available in Structural Engineering laboratory at 
IIT Kanpur was ideally suited, for the purpose. Hydraulic 
jacks were fitted on the channels or box sections which 
were held dovm to the tes'c floor with the help of mild steel 
rods, passing through the fifty cm centre to centre holes 
of the test floor. These mild steel rods in turn transferred 
the load to the walls at desired level. In all, nine jacks 
were used, 3 on each channel or box section. The details 
of loading amngement for applying the compressive and 
tensile type of loads are shown in Eigs. 3.3 nnd 3.4 
respectively. 
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Ilie jacks and the associated hydraulic system 
were calibrated r 3ing a proving ring of 20,000 kg. capacity. 
The specimens were landed and unloaded several times to 
about 25 to JO percent of first crack load to check the 
performance of loading system and measuring devices. The 
loading was incremented at prescribed interval so as to 
have adequate readings in the elastic and post cracking 
stage. Electrically operated hydraulic loading was with 
remote control and therefore loads could, be recorded upto 
failure accurately. However, near the failure load, measuring 
devices had to be removed in order to protect them from 
damage due to possible collapse, 

3,1.4 Instrumentation 

In order to measure strains in brickwork and overall 
deflection of the composite construction, strain measuring 
devices were suitably mounted on the test specimen, 

3. 1.4.1 Dial gauges 

The overall deflection of composite construction at 
various load levels have been recorded with the belp of 
dial gauges of 0.01 mm ^eastcount mounLed at the bottom of 
E-.O. beams. In all, five dial gauges were used for the 
purpose. The dial gauges were placed at midspan and at 
50 and 100 cm from midspan on either side. 
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3. 1.4* 2 Strain measuring device 

Metalstuds were initially tried to measure the 
Ic^teral and vertical strains on the surface of composite 
construction, but had to be abondoned because they gave 
erratic readings. After this, attempt was made to measure 
strains in the horizontal direction by 15 cm single wiere 
■120 ohm electrical resistance strain gauges and 1.5 cm 
flat grid 120 ohm electrical resistance strain gauge in 
vertical direction with no success, Finally strains were 
recorded by fixing the dial gauges between two points 50 cm 
centre to centre at each horizontal leve"^!. For fixing the 
dial gauges 5.5 mm dia holes were drilled in the masonry 
with the help of an electric drill using masonry drill bits. 
One 80 mm long roofing bolt was fixed in each hole with 
help of a plugging compound. Between these two roofing 
bolts, the dial gauge was held. The arrangement can be 
clearly seen in the photographs for cracking pattern. In 
most of the cases, these gave satisf dctory readings, 

3.2 PiiRiU'IuTRiC STUBIES 

Experimental study of the interaction between R.G. 
beam and brick masonry wall has been conducted for both 
compressive and tensile loading for the different parameters 
such as (l) ratio of cement and sand used in mortar for 
brickwork (2) height to span ratio of composite system 
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(3) symmetric and unsymmetric openin^^s in the masonry wall. 

3.2.1 Yariation in Cement Sand Ratio in Mortar 

In order to stady the efCect of mortar strength 
on interaction of composite construction, five cement sand 
ratios of 1:8 and ranging from 1:6 to 1:3 have been experimented 
upon, keeping a constant height to span ratio of 0,33. Ihe 
size of ten specimens numbered 1 to 10, five each tested for 
compressive and tensile loading is shovm in Pig, 3.5* Ihe 
bending reinforcement is provided for a total working load 
of 9.0 tonxies and assuming the lever arm to be 0,9 times in 
depth of composite system, considering it to be an under- 
reinforced brick beam with zero tensile strength of brick 
work and concrete, 

3.2.2 Variation in Hei '•ht to Span Ratio 

The effect of variation in the height to span ratio 
on the ultimate load carrying capacity of the composite 
system has been studied on ten specimens, five each under 
compressive and tensile loading for height to span ratios 
of 0.25, 0,33, 0.4, 0.5 and 0.8. The bending reinf or cement 
had to be varied for different height to span ratios as 
described in preceding subsection.- Hovrever, the same bending 
reinforcement was provided for height to span ratio of 0.5 
and above as obtained from calculations for 'B./L of 0,5. 

The table on next page gives the details of 
specimens tested for the purpose. 
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3.2,3 Va'ciation in Size and -location of Openings 

The walls c-f such composite construction shall have 
door or window openings at suitable locations. How do these 
openings influence the interaction of composite construction 
IS the subject of experimental study. In the present work, 
a door opening of 1m x 2m and a window opening 1m x 1m has 
been considered. While the door opening starts right above 
the R.C, beam, the sill of window opening has a height of 
1m above the R.C, beam. Both of these openings have been 
placed symmetrically about the vertical centre line and 
unsymmetrically at a centre of 95 cm from one end of the 
specimen for purposes of experimentation, A precast concrete 
member, 8 cm thick and carrying 3 plain mild steel bars of 
1 0 mm dia formed the lintel over the openings. The details 
of nine specimens numbered 19 to 27 are shown in Pig, 3*6, 
Specimens number 19 to 23 were cast in mortar 1:6 and have been 
tested under compressive loading. The rest of the specimens 
were cast in mortar 1 .*3 and tested under tensile loading. 
Specimens number 23 and 27 are similar to specimens number 22 
and 25 respectively except for the fact that thickness of 
R.C, beam in specimens number 23 and 27 is increased to 16 cm. 
This was necessitated because specimen number 22 failed under 
its self -weight while making preparations for testing. All 
specimens had a height to span ratio of 0.8, 
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3.3 RESULTS 

file resuj.ts of various tests conducted are reported 
and discussed in this section according to parameters studied. 

3 . 3.1 Effect of Variation in Gement Sand Mortar Eor Briclc Work 

The oDserved test data on deflections at the bottom 
of R.G, beams at various points along the span, are given 
in Tables 3.1 to 3.5 for compressive loading and 'Tables 3.6 

a 

to 3.10 for tensile loading respectively, Figs, 3.7 and 3*8 
show the plot of load versus deflection for various specimens 
tested under compressive and tensile loading respectively. 

Load deflection curves are almost linear in the initial 
stages and become nonlinear after the appearance of first 
crack. However, in case of load deflection cux’ves for 
specimens tested under tensile loading, there is slight 
deviation from linearity even in the initial stages before 
separation sets u-. 

Tables 3.11 and 3.12 reports the observed test 
data on longitudinal strains for specimens cast in mortars 
of 1:3 and 1:6 respectively tested under compressive loading. 
The test data on longitudinal strains for specimens cast in 
mortars of 1:4, 1:5 and 1:6, and tested under tensile loading 
IS given in Tables 3.13, 3.14- and 3.15 respectively. The 
strains at a cross section of composite construction, at 
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mid span, have been plotted in Eigs. 3.9(n) and 3.9(t') for 
specimens tested under compressive and tensile loading 
respectively. Yariation of strains is linear at this 
cross section, Erom the linear load deflection curves and 
the linear variation of longitudinal strains before cracking 
starts, it can be concluded that composite material behaves 
linearly in the elastic range. Eigs, 3*l0(a) and 3 , 10(b) 
show the cracking pattern for specimens number 1 to 5 and 
6 to 10 and tested under compressive and tensile loading 
repectively. Specimens number 1 and 2 cast respectively 
in 1:3 and 1:4- mortar and tested under compressive loading, 
cracked near the mid span, indicating thereby a bending 
failure. Specimen number 3 cast in 1:5 mortar and tested 
under compressive loading has a mixed type of crack pattern 
i.e, failure is caused jointly by bending and diagonal tension. 
Specimens number 4 and 5 cast in 1:6 and 1:8 mortar failed 
due to diagonal tension as the cracks in these specimens 
developed near the supports only. As is seen from Eig, 3.10(t))y 
all the five specimens numbered 6 to 10 tested under tensile 
loading failed by horizontal separation at different heights 
and subsequent development of cracks due to diagonal tension, 

Tables 3.16 and 3.17 give the first crack load, failure 
load, and the load factor for specimens tested under compressive 
and tensile loading respectively. It is observed from these 
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tables that the first crack load for specimens tested under 
compressive loading varies from 12.26 to 22.85 tonnes (including 
self weight) , while for specimens tested under tensile 
loading it varies from 5»6 tonnes to 10,7 tonnes (including 
self weight). The load carrying capacity in elastic range 
of the thin concrete beam supporting the brick masonry is 
just its self weight which is around 150 kg. Therefore the 
composite construction is taking a load in the elastic 
range v^fhich is appro ximateT-y 37 to 152 times more than the 
load carrying capacity of the beam it self depending upon 
the mortar used and the loading condiuion. In fact, the 
contention of present work is to demonstrate that for the 
composite construction under investigation, thinnest possible 
R.G, beam can be used which is sufficient enough to provide 
suitable bond and cover to the bending reinforcement and hold 
the single legged vertical connectors in ;^ssition. hot with 
standing these high ratios, it is clearly observed that the 
composite construction behaves as a single composite structural 
element in the elastic range, Furthermore, richer the 
mortar mix, greater is the interaction in the elastic range 
between the masonry and h.C, beam connected through single 
legged Z -shaped stirrups. 

The total failure load for specimens tested under 
compressive and tensile loading varies from 12.26 to 35.24 
tonnes and 9.28 to 18.7 tonnes respectively, while the 
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ultimate load carrying capacity of R.C, beam itself is 
only 800 kg. Therefore, the interaction observed in the 
elastic range extends even beyond. In other words, richer 
the mortar mix greater is the interaction upto failure of 
the composite structural element. 

ilv4 

It IS seen from Table 3.17(the load factor for 
specimens tested under tensile loading varies from 2.08 to I. 03 . 
All specimens, except the one cast in 1:5 mortar, have a load 
factor of less than 2,0. I.S. code specifies a load factor 
of 1,5 for dead load and 2.2 for live load and therefore 
mortar leaner than 1:3 are not suitable for use when loading 
i IS tensile. 

For specimens tested under compressive loading the 
load factor is more than 3 for all specimens cast in mortar 
1:6 or richer but it suddenly drops down to 1,36 for specimen 
cast in 1:8 mortar. It is also observed from Table 3*16 that 
the first crack load and failure load for specimen in 1,8 
mortar are exactly the same meaning thereby a sudden failure 
of composite system in 1 ;8 mortar. The reason for this is 
the bond failure between brickwork and the vertical connectors. 
Therefore, it can be concluded that the use of 1:8 mortar 
for brickv^ork will not provide the desired interaction. 

Further more, keeping in view that use of 1:6 mortar achieves 
a load factor of about 3> it can be safely recommended to 
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cast such composite structural element in 1:6 mortar for cases 
where loading is compressire. 

B.-ised on the above observations, subsequent parametric 
studies have been conducted on specimens cist in mortars of 
1:6 and 1:3 for compressive and tensile loading respectively. 

3.3«2 Effect of 'J'aria'cion in Height to Span Ratio 

the observed test data on deflections at the 
bottom of R.O. beam at various points along the span are 
given in Tables 3.4 and 3.18 to 3*21 for specimens tested 
under compressive loading while the same for specimens tested 
under tensile loading are given in Tables 3.6 and 3.22 to 
3.25. Rigs, 3*11 and 3.12 show the load deflection curves 
for compressive and tensile loading respectively, Ror the 
load deflection curves, it is observed that deflections go 
on reducing as the height to span ratio increases, meaning 
thereby that the element is becoming stiffen with the increased 
height to span ratio. 

Tables 3.12 and 3.26 to 3.29 report the observed 
data on longitudinal strains on the specimens under compressive 
loading and Tables 3.30 to 3.32 under tensile loading. The 
strains at a cross section at mid span have been plotted 
in Rigs. 3, 13(a) and 3.13(b) for specimens tested under 
compressive and tensile loading respectively, The variation 
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of strains is linear npto a height to span ratio of 0 , 4 , 

It IS a bilinear curve for a height to span ratio of 0,5 for 
specimens under compressive loading and is linear under 
tensile loading, For a height to span ratio of 0.8 the 
variation of strains is not clearly obtained. The maximum 

strains at the top and bottom of specimens under compressive 

—5 

loading in the latter case are of the order of 5 x 10 ~ cm/ cm. 
The least count of measuring device leads to a strain of 

[T 

2x10” cm/cm, A couple of dial gauges upto quarter 
height above the R.C. beam recorded constant values leading 
to zero strains meaning thereby that the strain is less than 
2 X 10”^ cm/cm. In the top portion of this specimen the 
strain variation is clearly nonlinear. For the specimen 
number 18 of height to span ratio of 0,8 tested under tensile 
loading meaningful data could not be obtained. 

Figures 3.14(a) and 3.14(b) show the cracking pattern 
of specimens nuiiber 11 to 14 and 15 to 18 for compressive and 
tensile loading respectively. Specimens number 11 and 12 
having height to span ratio of 0.25 and 0,4 tested under 
compressive loading developed cracks near supports meaning 
thereby a failure in diagonal tension similar to specmen 
number 4 having height to span ratio of 0,33 cast in 1 :6 mortar 
as discussed before. Specimens numbered 13 and 14 having 
a height to span ratio of 0,5 and 0,8 failed due to diagonal 
tension and crushing of concrete and brickwork at the supports. 



105 


Piio-tographs showing the failure pattern for these specimens are 
affixed as Pig. 3.15. All the four specimens numbered 
15 to 16, tested under tensile loading, failed by horizontal 
separation at different heights and subsequent development 
of cracks due to dia,gonal tension. I'aoles 3.33 and 3.34 
give the first crack load, failure load and the load factor 
for specimens tested under compressive and tensile loading 
respectively. Prom these tables it is observed that higher 
is the height bo span ratio higher is the load factor. Under 
compressive loading, the load factor for the specimen of 
height to span ratio of 0,25 is 1,8 while for others it 
varies from 3 to 6, The reason for a low value of failure 
load for specimen of height to span ratio of 0,25 is the 
possible bond failure between the vertical connectors and 
the brickwork i.e. a bond length of about 70 cm in 1 :6 mortar 
is sufficient to transfer loads corresponding to the load 
factor of 1,8 , It is, therefore, concluded that in order 
to achieve a load factor of 2 or above a rich mortar should 
be used for height to span ratio of 0,25. 

Por specimens tested under tensile loading the 
load factor varies from 1,61 to 2.98. The load factor is 
below 2.0 for a' height to span ratio of 0.25 while, for 
others, it is more than 2.0. Por this height to span 
ratio, one additional specimen Uo, 28 was tested. The 
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specinieii was identical to specimen number 15 except that 
ib had 6 mm dia 24 connectors provided at a uniform spacing 
of 13 cm centre to centre instead of only 13 connectors. 

The observed data on the deflections for this specimen is 
given in Table 3.35. 'the load deflection curve and the crack- 
ing pattern is shown in Pigs, 3.16 and 3.17. In this case, the 
deflections are lower as compared to those for specimen 
I!'o,15 having 13 vertical connectors. Por example, the 
deflection at mid span at a load of 6.77 tonnes is 1.47 mm 
for specimen number 15 "vdiile for specimen number 28 it is 0.81 
for the same load, furthermore, the horizontal separation 
starts at a much later stage of loading in specimen number 28 
i.e. at a load of 13.5 tonnes instead of 6.77 tonnes in 
specimen number 15, showing thereby increased resistance to 
bond slip due to use of more connectors. The load factor 
for this specimen is 2.08. It is, therefore, concluded 
that in case of tensile loading and height to span ratio 
of 0.25, more number of vertical connectors have to be used. 

Motivated by the above results, two additional 
specimens numbered 29 and 30 were cast for a height to span 
ratio of 0.33 in order to study the effect of the variation 
of size and nuimber of vertical connectors on the load 
factor. Specimen number 29 had 13 number 8 mm dia and 
specimen number 30 had 20 number 6 mm dia plain mild steel 
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vertical connectors respectively. Thus specimen number 29 
had same number of increased diameter connectors while 
number 30 had increased number of same diameter connectors 
as specimen number 6, The observed data on deflections 
for these specimens is given in Tables 3.36 and 3,37 
respectively. Load deflection curves for these specimens 
are shown in Fig, 3.18, Table 3.38 gives the first crack 
load, failure load and load factor for specimens numbered 6, 

29 and 30 for the purposes of comparison. It is observed from 
Table 3.38 that in specimen number 29 separation started at 
a load of 6,768 tonnes and it failed at 12.69 tonnes while 
in case of specimen number 6 separation started at a load of 

8.5 tonnes and failure took place at 16,5 tonnes, indicating 
thereby that increase in the sizd of connectors has reduced 
the load carrying capacity. This can only be explained by 
the fact that 8 mm dia connector has not been able to develop 
sufficient bond with brickwork through mortar joint about 

1.5 cm thick normally used. In specimen number 30, separation 
started at a load of 13,536 tonnes and it failed at a load 

of 19,0 tonnes indicating higher load carrying capacity 
with increased mmiber of connectors. From this limited 
study, conducted on specimens numbered 28, 29 and 30 , it 
can be concluded that, 6 mm dia connectors provided at 
spacing closer than 25 cms yields better interaction in 



108 


situations where loading is tensile. However, in order 
to achieve spacing of connectors lower than 25 cms, breaking 
of the bricks is inevitable, 

3. 3*3 nffect of Variation in Size and location of Openings in 
Brickvfork 

tables 3.39 to 3.42 give the ooserved test data 
on deflections for specimens number 19, 20, 21 and 23 and 
tested under compressive loading. Tables 3*43 bo 3*46 
report the observed data on deflections for speciiiiens 
number 24 to 27 and tested under tensile loading. Big. 3»19 
and Fig. 3.20 show the plot of load versus deflections for 
specimens with openings tested under compressive and tensile 
loading respectively. In addition these figures also show 
the load deflection curves for solid walls o-j- height to 
span ratio of 0.8 tested under respective loading for purposes 
of comparison. Meaningful data on longitudinal strains 
could not be obtained because of lew accuracy of measuring 
devices. Figs. 3.21 and 3.22 show the cracking pattern 
of specimens number 19 to 23 and 24 to 2? tested under 
compressive and tensile loading respectively. In case of 
compressive loading, specimen Ho. 19 with a symmetric window 
opening failed in diagonal tension and crushing of brickwork 
and concrete at supports. Specimen no. 20 with a symmetric 
door opening failed in diagonal tension only. Specimen 
no. 21 having an unsymmetrj^ window opening also failed in 
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diagonal tension accompanied by shear as is evident from 
the crack pattern. Specimen no, 22 and 23 failed due to 
shear only. Photographs affixed as Pig. 3.23 show very 
clearly the failure patterns discussed. 

All specimens tested under tensile loading failed by 
horizontal separation and subsequent development of cracks 
due to diagonal tension, 

Tables 3.47 tuid 3.48 give the first crack load, 
failure load and load factor for specimens tested under 
compressive and tensile loading respectively. Prom fable 
3.47 and Pig. 3.19, it is observed that S37mmetric openings 
do not affect the load carrying capacity of the specimen 
under compressive loading. However, unsymmetric openings 
affect the load carrying capacity signif icantly, The load 
factor for an unsymmetric window opening is 4.1 against 5.9 
for symmetric window opening. The unsymmetric door opening 
affects the load carrying capacity adversely. Specimen 
number 22 failed under self -weight while making prepcirations 
for the test. The specimen number 23, which had a 16 cm thick 
R.C, beam failed at a load of 10,04 tonnes thus giving a 
load factor of only 1.61. 

The load carrying capacity of specimens under tensile 
loading is reduced due to the presence of symmetric or 
unsymmetric door and window openings. The load factor for 
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specimen witli symoietric and unsyiometric window openings 
in 2 * 4 . 4 - and 2,18 respectively. It may not be out of place 
to recall th.at specimen number 6 having height to span ratio 
of ,35 and cast in mortar 1:3 has a load factor of 1,98* fhe 
sill of the window openings as described earlier was hept 
at 1 m height above the R,C, beam meaning thereby tho.t sill 
height to span ratio is *33. It is, therefore, concluded 
that specimens with window openings fail at a load only 
slightly greater than those whose height is eq.ual to the 
sill height. 'Ihe load factor for specimen with a symmetric 
door opening is only 1.57 as against a load factor of 3 .O 
for corresponding specimen without opening. The load factor 
for specimen number 27 which has an increased thickness of 
R.C, beam viz. 1 6 cm and a symmetric door opening is 2.17. 
Rrom this it is concluded that inorder to have a load factor 
of 2.0, such composite construction carrying S 3 nnmetric door 
openings be carried out with 16 cm thick R.C, beams. 
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^ — 
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TABLE 3.11: STRAINS 


H/L MORTAR LOADING 
^35 ^TTS “cOMPiESSITi 


Pressure on jacks 

! 000 

200 

400 

600 

800 

1000 

in psi 



! 

1 






Total load 

in tonnes 

0.0 



4,13 

8.26 

12.39 

16.52 

20.65 



Reading 
of dial 
guage 

1 

11,59 

11.61 

11.635 

11»66 

11.685 

11.71 



00 










Strain 

0.0 

-4.0 

-9.0 

- 14.0 . 

-19.0 

- 24.0 

s i 


xIO*"^ : 








*" n 







O 

0 -P 

bjO 

S s 


Reading 
of dial 

9.34 

9.355 

9.37 

9.39 

9.41 

9.43 

t3 o 

Ci3 ^ 

OJ 

• 

gauge 








cn 








a 


Strain 

0.0 

-3.0 

-6.0 

-10,0 

-14.0 

-18.0 

■H O 


xl0"^ 















o Ai 

® ^l 


Reading 







^ g 
o P ^ 


of dial 

7 , 65 

7.66 

7.67 

7.68 

7.69 

7.705 

H 

-P -P q 

ClJ CQ H 


gauge 

1 






• 








o H p 
o p pq 

DA 

Strain ; 

io.o 

-2.0 

-4.0 

-6.0 

-8.0 

-11.0 



X10"5 i 

I 

i 

1 








Rea,ding ' 

1 








Of dial 12.16 

2.165 

2.17 

2.175 

2.18 

2.185 


O 

• 

gauge 

1 

( 







O 

LH 

Strain ! 0, 0 

j x10“^ I 

-1.0 

-2,0 

-3.0 

-4.0 

-5.0 


Table contd...on page 1 22 



122 


TABLE 3.11 contd, . . 


Pressure on ^acks j 000 200 400 600 800 1000 

in psi 

Total load in tonnes | 0,0 4*13 8,26 12.39 16,52 20,65 

Heading ! 

[ of dial |1.12 1.12 1.12 1.12 1,12 1.12 

i gauge i 

'+-1*1 I 

° ,5^: strain 1 0.0 0.0 0.0 0.0 0.0 0.0 


05 S 
H O 

'H 

o M 
<D U 
ri g o 

H 05 15 
-p -p o 
o5 m H 
O H i4 
o n pq 

hi'—' 


b 

Heading ! 
of dial t 

p 

I 

gauge I 

I LA 

I 

Strain i 

x1 0”^ 



Reading 
of dial 

c- 

I * 

gauge 

O 

CO 

Strain 

I 

xl0“^ 

! 

n 

I 

Reading | 
of dial j 

T- ! 

gauge 1 

I ^ ' 

1 1 

UD j 

cr\ 

1 Strain i 

r 

xlO"^ 1 


1.12 1.12 1.12 


Remark 


0.0 

0.0 

0.0 

0.0 

6.035 

6.03 

6.02 

6.0i: 

0.0 

1.0 

3.0 

4.0 

0.17 

0.16 

0.145 

0.13 

0.0 

2,0 

5.0 

8.0 

Separation 

starts 



6.015 6.005 5.995 

4.0 6,0 8.0 


0.13 0.115 0.10 

8.0 11. Q 14.0 
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T iiBLE j . 1 2 : 


STRAIiTS 


E/L MOR'fAR LOiDIRG 

0.33 n?" coi'ffifiEisrra 


Pressure on oacks 
in psi 

! 

000 

100 

200 

300 

400 500 

Total 

load in tonnes | 

0.00 

2.065 

4.13 

6.195 

8.26 10.325 

Vi 

0 ! 

04 1 

; 

Reading 
of dial 1 
gauge 

7.68 

7.71 

7.74 

7.77 

7.80 

7.83 

0 1 

4 ^ 

a 

0 

CO 

• 

Strain i 
xlO“5 

0.00 

- 6.0 - 

- 12.0 - 

-18.0 - 

- 24.0 . 

- 30.0 

u \ 

0 

CNJ 

• 

(T\ 

Reading 
of dial 
gauge 1 

10.36 

10.38 

10.4® 10.43 

10.455 

10.48 

H 

0 

0 

Strain 
xlO”^ j 

i 0.00 

I 

-4.0 . 

- 9.0 - 14.0 - 

-19.0 . 

- 24.0 

•H 

3 

CQ 

\ 

\ 

\ 

! 

1 

1 

VD i 

Reading I 
of dial j 
gauge I 

6.62 

6.635 

6 .65 

6,67 

6.685 

6.705 

0 

fciD 

1 

Cb 1 

• 

to 

I Strain { 

i X 10 "'' j 

0.0 

-3.0 

-6.0 

-10.0 

- 13.0 

- 17.0 

H ^ 

fl ?4 
0 

VI Pi 

0 

• 

I Reading 
of dial 
gauge 

14.62 

14.63 

14.64 

14.65 14.66 

14.675 

0 

0 

H 

0 ?-I 
rH PP 
-P 

0 

LO 

Strain 

I 2;10”5 

j 0.0 

1 .. _ 

-2,0 

-4.0 

-6.0 

-8.0 

-11.0 

05 

0 

0 

VI 

1 f 

Reading 

1 of dial 

1 gauge 

« „ .... 

16.98 

16.985 

16.99 

16.995 17.00 17.005 


] LA 

1 ^ 

J 

Strain 

i X 10"^ 

I 0.0 
! 

-1.0 

-2.0 

-3.0 

-4.0 

-5.0 
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Table 5.12 contd.. 




in psi 


Ph 
U1 o 

CD -P 
tiO 

3 S 
cd o 

U 

I — I 
05 

H O . 

p ^ N 
Im a 
'H H p 
O ^ 

_ p u 
flop 
o ri H 
H 05 fH 
-P4= P 
05 m 
O r^ CH 

op o 
p'--' 


on 

gacks 

j 

I 

00 

I 

100 

200 

300 

400 

500 

ad in tonnes 

o 

o 

• 

o 

2.065 

4.15 

6.195 

8.26 

10.325 

tr- 

• 

O 

OD 

Reading I 
of dial I 

r 

j 7.15 

7.15 

7.15 

7.15 

7.15 

7.15 

gauge j 

~ • I 

Strain ! 

o 

o 

« 

o 

0.00 

0.00 

0.00 

0.00 

0.00 

i 

T— I 

xl0"5 

Reading 
of dial 

5.97 

5.965 

5.955 

5.95 

5.945 

5.935 

• 

^ I 
^ i 

i 

gauge 

1 Strain 

0,00 

1.0 

3.0 

4.0 

5.0 

7.0 

! 

I 

1 X10”^ j 

I 









T:IBLE 3.13- STRAINS 


Pressure on ^aclcs 
in psi 

Total load in tonne 


o 

Pa 

o 

-p 

a 

o 

u 

a 

o 

•H 

0) 

O 


P 
CQ 

8 
0) 

N) 

I 

a 

•H 

O ?H 
O ^ 


o 

p 



heading j 
of dial 
gauge 


Strain 

1 

xlO-5 

j 

Reading 
of dial 
gauge 

CM 


• 

CTv 

Strain 


xl0~5 


Reading 


» of dial 

1 

S VD 

j g tage 1 

K\ 

* .strain ' 


1 xlO-5 1 

1 

.) 

1 

I Reading 
of dial 


O 

1 * 

gauge 

1 * 

> o 

1 LTN 

\ 

] 

I 

Strain 

X10-3 


H/L 

MORfiR 

lOiiDING 


0.53 1:4 


TjSRSIIE 



OOO 

200 

400 

600 

800 

0.0 

1.69 

3.38 

5.08 

6.77 

10.41 

10. 44 

10.47 

10.50 

10,53 

0.0 

-6.0 

-12.0 

-18,0 - 

•24.0 

2.34 

2.36 

2.38 

2.405 

2.43 

0.0 

-4.0 

-8.0 

-13.0 

-18.0 

0.47 

0.485 

0.50 

0.515 

0.53! 

0.0 

-3.0 

-6.0 

-9.0 

-13.0 

1.23 

1 .24 

1.25 

1.26 

1.27 

0.0 

-2,0 

-4.0 

-6*0 

-8.0 


900 

7.61 

10.545 

- 27.0 

2.445 

- 21.0 

0.545 

-15.0 

1 .275 
-9.0 
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Table 3*15 coirtd.... 





Location of Dial G-auge 
(Distance in cm from top of Brickwork; 


i‘^i 


1A.BLE 3.14 


H/L ' MOBTiR LOADING 

STRAINS “ — — — — 

0.33 1:5 TBNSILE 
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Table 5.14 contd..,. 


Pressure on ^acks 
in psi 


000 200 400 600 700 800 


Total load in tonnes 0.0 1.692 3.4 5.07 5.92 6*77 


Reading 

Of dial ! 9.28 9.285 9.29 9.295 9.295 9.30 


O 

ft 

o 

CD -P 

CtJ o 
0:5 U 
Ch 
I — I 

n3 S 
•H O 

” a 

S'":? 

^ <D M 

^ 9 S 
o S ^ 
H S 
+3 -p o 

05 CQ r-| 

O H H 
o nfd 


^ j gauge 

m . 

'n I Strain 
xIO*"^ 


Reading 


o 

CO 


! 


xiO* 


0.0 - 1.0 - 2.0 - 3.0 - 3.0 - 4.0 


of dial 
gauge 

4.86 

4.86 

4.855 

4.855 

4.855 

4.85 

Strain 

0.0 

0.0 

1.0 

1.0 

1.0 

2,0 




J 




I 

Reading 
of dial 

1 .27 

1.26 

1.25 

1.235 

1.23 

1 .225 

gauge i 







Strain 

0.0 

2.0 

4.0 

7.0 

8.0 

9.0 
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TABIE 3.15: STRA.MS 


li/T— MORTiiH " LOiDIlTG 
0.33 1:6 TEIISILE 


Pre 

ssare on oacks 1 

000 

200 

o 

o 

500 

600 

700 

in 

psi 








Total 

load in tonnes 



0.0 

1.692 

3.384 

4.23 

5.07 

5.92 

cm) 

! 

I 

i 

Reading 
of dial 
gauge 

1.36 

1.40 

1.445 

1 .47 

1 .495 

1.52 


CO 








H 


• 

Strain 

0.0 

- 8.0 . 

-17.0 . 

- 22.0 

- 27.0 . 

- 32.0 

u 



xl0“^ 







O 










1 



Reading 







o 

H 



of dial 

2.76 

2.79 

2.83 

2.85 

2.87 

2.89 

U 

PQ 


CVI 

• 

gauge , 







0) 


(T\ 








bJDCH 

O 


-r- 

Strain 

0.0 

-6.0 . 

-44.0 . 

-18.0 

- 22.0 

-26.0 

05 

cb Pd 
o 



XI 0-5 1 







H -P 










CU 



Reading 

I 






U 


! 

of dial 

I 7.31 

7.335 

7.365 

7.38 

7.395 

7.41 

o 

1 

• 

gauge 

1 






CD 










O P 



Strain 

1 0.0 

-5.0 . 

- 11.0 ^ 

- 14.0 

- 17.0 

- 20.0 

rH 05 
-P P^ 


1 

x 10"5 

\ 






05 CD 


L 








O *H 

33 



Reading 
of dial 

2.005 

2.02 

2.04 

2.05 

2.06 

2.07 



o 

gauge 

! 








o 


lo.o 








! ir\ 

Strain 

-3.0 

-7.0 

-9.0 

- 11.0 

- 13.0 


X10-5 I 
s 
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Table 3»15 contd.... 


in psi 


ij 

^ I • 

I m 


H a o 

p p CO 

u 

CH 

O 

CD 

^ ^ I 

o rig [ 

H 05 O I 

-P -P 
03 CO Pi 

O ‘H H • 

Op I, ^ 

P'^ y CJ^ 


L jacks ' 

000 

200 

400 

500 

600 

700 

in tonnes 

j 0.0 

1.692 

3.384 

4.23 

5.07 

5.92 

Reading 
of dial 
gauge 

12.54 

12.545 

12.555 

12.56 

12.565 

12.57 

Strain 

0.0 

-1.0 

-3.0 

-4.0 

-5.0 

-6,0 

X10 ^ 







Reading 
of dial 

1 .76 

1.76 

1 .76 

1 .76 

1.76 

1.755 

gauge 







Strain 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0 


Reading 


of dial 
gauge 

6.21 

6.20 

6.19 

6.185 

6.18 

6.17 

Strain 

0.0 

2.0 

4.0 

5.0 

6.0 

8,0 


I Separation starts 


Remark 
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--c\ 

EH 


CQ 

H 

^ S 




CQ 

U 

05 

S 

CD 

Ph 


hq 

u • 

O Pq j 
t5 4^ * 

05 O EH j 

o o5 

(i 


H 

o5 


H 
O 

o5 

U 

O nz^ 
CD 05 

u 

rH 05 


CQ 

-P 

u 

o 

CO 




CD 

c! 


qn 


H 

PH 




isD 1 

o5 

0 

0 

M 

bD Pi O 

P4 


P 

p! 


•H r-ci 


P 

P 

H 

•p 

1 


H 

o 



o 

0 

fi 

Pi 

o 

Pi 

o5 

Xi 

Pi 

03 

0 

Pi 

*H 

P* 


P 

pp 

o 


o5 

0 

o 


rH 

M ^ 
U 05 
o o 
^ H 


-F 

WbDu} 
ri H CD 

0 H 0 • 

fj rg ^ S l^ 

H 5=^ O • 

o5 H o:^ H ^H -P Cq 

-P H o5 o rH • 

O 05 O Pi 0 EH 

lEH ChH H 0«P ll 


0 



0 

0 


V- 

tP- 

CO 

KD 

H P 

p 

O 

cn 

O 

IP- 

O 

H o5 

Ph 

• 

« 

« 

• 

« 

o5 o 

Pi o 

m 

CO 

IP- 


o 

Gq H 

H P 

tn 

CM 

CM 

CM 



•H 

0 -P 
4-*^ q-l 0 
05 O Ci5 

n O 


I 


CO 

r- 

CM 


VX) 

}• 

CM 


o o5 
-P 


Pi 

0 


^c^ 


C- 

CO 


in 


0 e g Pi 
PH H 3 0 
CQ o Is; 


05 

^ 0 

o5 1 — 1 
0 H H 
Pi o5 

bO 0 
Pi -P Pi 
H Pi o 

O H 
O Ph 0 

ri 

P P) 0 
O 04^ 


in 

in 

in 

o 

CD 

cn 


CM 

o 

m 

• 

• 

• 

« 

* 

tn 

tA 

tn 

tn 

T- 



v— 


00 

MD 

CM 

V- 

CM 

cn 

CM 

• 

• 

• 

• 

• 

tn 


cn 


CM 

m 

tn 

CM 

CM 

'T- 


P M 

0 

0 

in 

CM 



VD 

0 o p 

pj 

VD 

in 

in 

tn 

o 

P 05 03 

pI 

• 

• 

• 

• 

« 

•H ri o 

p! o 

o 

VD 

in 

CM 

o 

pH op 

H P 

CM 

't- 

V- 





00 

IP- 

tp" 

IP- 

00 

p! 


P- 

tP- 

!P- 

tP- 

p- 

H 


« 

• 

• 

• 

• 

0 P 


m 

IP- 

IP- 

CJ\ 

tn 

p p 0 


T • 

• 

• 

• 

• 

03 O 0 

( 

1 

cn 

CM 

o 

in 

CO 

n EH 

! 

V- 


CM 




IP- 

IP- 

CO 


MD 


CO 

c- 

» 

CM 

00 


00 


CM tn 


in 
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3 y 

rH 

H 

rH S I 

05 O CH I 

-pH O 

-P 

O 03 o i 

H ‘H 
H n5 CO 
U 

000 
W -P 


I I 

0 o 

1 \ 



PH 

1 


H 

H 

h 

S i 

3 

0 

1 

P 

^ i 


1 

4 P 

i 

\ ; 


M j 

0 1 

it 


O l-^ I H 

^ -p « 5 ^ 

05 o pq ; p 05 

O o5 *500 
iP p E“i ; H 

, P 

H 

0 H 0 

P rc5 ^ • 

H :::3 ^ 

o5 H n:3 rH P ♦ 

-p H 05 O H |3q 

o o5 o f 0 | 0 • 

P «H H 0 EH 


IT-- VID 
tOs 


0 

H" 


Lr\ 

00 

t-” 


CM 


CM 

• 

• 

» 

• 

• 

00 

ir\ 

K>i 

CM 

(T\ 


U 

0 

H- 

0 

ir\ 

CO 

pi 

in 

CM 

0 

V- 

0 

H iiJ 

• 

• 

• 

« 

• 

H o5 

VJD 

ca 

r— 

0 

!>• 

05 0 

\ — 

v* 

T- 

*v — 


pq H 1 






Pi 






*H 0 






0 






0 OZJ ?H 

VO 


ir- 

CM 

0 

H-3 at 0} s 


VO 

D- 

cn 

H" 

0 P 0 0 

• 

• 

« 

• 

* 

H 0 H -P 

CD 


VO 

La 

ba 

SP 







cr- 

CO 

00 

CD 

D 

0 'H 


t- 

l^- 


la- 

P P P 

• 

• 

« 

« 

• 

0 0 


CM 

CM 

ta 


P 0 

« 

« 

« 


♦ 

P 

C- 

0 


VD 

ba 



CM 

CM 


V— 

SP 1 

C— 

00 

00 

00 

CO 


tr- 


t>- 

tc- 

lo- 

H 

• 

• 


« 

« 

0 P i 

CM 

r- 

r- 

CM 

ba 

-P CH 0 1 

« 

• 

p 

• 

• 

05 0 o5 

V- 



La 

r- 

p 0 1 

1 CM 

CM 

CNJ 



1 1 

i 

CO 

H" 

tr\ 

CO 

CO 

0 05 

• « 

«« 

« a 

• « 

0 a 

gP i 




x*-* 

T— 


I 0 1 

0 a g p 

PM H P 0 

ca O ^ 


CO CTn 



c\5 

s 

CD 

Ph 



0 

P 

p 

P 

0 M 

1 — i 

P 0 

H 

P o5 

03 

ft g 

P-i 
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o 

h=i 


Si 

EH 

O 


K 


H 

CO 


O 

O 


v£) 


in 

CNJ 

o 
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H 
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O 




00 
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i CD 
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rH O 

P H S 
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I 0 

oij P tjjD 
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0 pl CH P 
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I 
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o5 
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0 

H w 

P H 05 
H O re} tlO 


0 ri 

H O 

^ ^ ^ 

0 p a 

PI o H S 


0 

H bO 
05 

Vh -H 05 
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0 

! — I ttO 
- „ cc 5 ^ 
0 p H c\5 
pej rH O re} tlD 


rd 

o5 

o 

p 


I 

0 

p 

p 


CQ 


P o 

p p 


0 CQ 

P ^ 

0 p o5 ri 0 

0 O *r-o*H p 


H o 

CO fef 


0 

ir\ 


VD 

1 

KO 

CM 




0 

T“ 

A 

p- 

VsD 

A 
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• 

m 

1 

1 

1 I 
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P 




CO 
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CM 
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• 

1 

1 

1 
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P- 

A 
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y- 




0 

(M 

P- 
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CM 

0 




• 

• 

• 

• 

• 

» 

1 

i 

1 

0 

0 

0 

0 

CM 

A 




CA 

A 

00 

A 

KD 

0 




0 

A 

A 

00 

A 

y^ 

1 

1 

1 

• 

W 

# 


• 
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v- 
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T— 
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V- 
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P 
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0 
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Table 5^26 contd.r 
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Table 3.30 contd.... 
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TABLE 3.31: STRAINS 


E/L 

MORTAR 

LOADING 
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1:3 
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Table 3.31 contd... 
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TiBIE 3.32 
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Table 3.32 contd... 
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T.4BLE 3.34: PIRJST CRACK AMD FAILURE LOAD 
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Note: 24 connectors of 6 mm diameter has been used 
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PLAN 



Fig-3‘1 (a) Mould for casting R.C. beams (Compressive loadihjl 
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Fig.3'1(b) Mould for casting R.C. beams *Tens,ie 
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Fig. 3-2 a Details of beam (Compressive loading) 
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Fig. 3*2 b Details of beam (Tensile loading) 
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Rg.3*3 Loading arrangement (Compressive loading) 







Fig. 3*4 Loading arrangement (Tensile loading) 















Specimen no. 2 2 


Specimen no. 21 and 26 


Fig. 3*6 Details of specimen with openings 
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Fig.3*8 Load versus vertical deflection curves 
(Tensile loading) 
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Fig.310 Crack pattern at failure (H/L-0*33)' 












H/L=0'8 


o Denotes ftrst 


Mortar 1'.6 


L -325cm 


Deflection in mm at mid span 


.H/L- 0-8 


H/L= 0-8 


Deflection in mm at 50 cm from Deflection in mm Qt 100cm 
mid span from mid span 

Fig.3H Load versus vertical deflection curves 
(Compressive loading) 
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Fig.3-14 Crack 



pattern at fai(ur 0 




Specimen no-U H/L^O‘8 , Mortar 1:6 
Loading - Compressive 

Fig.3-15 Photographs showing failure pattern 
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Rg.3’17 Crack pattern at failure 
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H/L=^0'33 
Mortar 1:3 


•L- 325 cm 


A - 20 connectors 6mnn dia 
B - 13 connectors 6mm clia 
C - 13 connectors 8mm dia 
o - Denotes first crack 


> 4 6 8 10 
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Deflection in mm at 50cm from Deflection in mm at 100 cm from 
mid spon rnid span 

Fig.3-18 Load versus vertical deflection , 
curves (Tensile loading) 
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FiQ.3-20 Load vs vertical deflection curves 





Fig.3-21 Crack pattern at failure 
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Specimen no. 21 Unsymmetnc window opening 
H/L = 0*8 , Mortar 1'.6 , Loading - Compressive 



Specimen no. 22 Unsym.metnc door opening 
H/L = 0-8 , Mortar 1 6, Loodm g -Compressive 


Fig.3-23 a Photographs showing failure patterns 
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Specimen no 25 symmetric door opeming R.C. 
beam 8 cm thick 

H/L = 0-8, Mortar 1:3, Loading - Tensile 



Specimen no 27 symmetric door opening R.C. 
beam 16cm thick 

H/L= 0‘8 , Mortar 1'3 / Loading -Tensile 


Fig.3-23 b Photographs showing failure patterns 
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ANilYTICJili IFVBSTIG-ATlOlf 


4.1 lilTRODUCTIOi^ 

Brickwork suppori:ed on R.C. beams is a composite 
structure wkicli has been idealized as orthotropic material 
in Chapter II, In order to carry out a rational analytical 
study of this class of structure for the purposes of studying 
crack propagation and obtain load deformation characteristics 
of the composite system, finite element method seems to be 
the most usefiil tool. Hence a computer programme based on 
finite element method incorporating nonlinearity in the load 
deformation response arising due to cracking, yielding and 
* crushing of concrete or brickwork and yielding of steel 
reinforcement has been developed in the present work, The 
structure has been modelled as a plane stress problem since 
the composite system considered in the present work is a 
two dimensional planar structure. 

The finite element method is the representation of 
a continuum by an assemblage of subdivisions called finite 
elements. The characteristics and properties of all elements 
are combined to analyse the structure. The simplest element 
used for two dimensional finite element analysis is a constant 



strain triangular element, For tv 70 dimensional bending 
problems, constant strain triangular and quadrilateral 
elements are not able to represent certain simple stress 
gradients and that causes inaccuracy in the solution (54)* 
Therefore the^e is need to use higher order elements. Wilson 
et al. ( 54 ) improved the accuracy of quadrilateral element 
by adding incompatkble displacement modes in the displacement 
field, The:^e extra displacement modes violate inter element 
compatability. The other higher order and simple element in 
use is a linearly varying strain triangular element, which 
has been used in the present work. 

The finite element method employed in this work is 
based on the stiffness (displacement) approach. Since the 
method is very well documented in standard text books (55>56, 

57 ) only a brief description of the steps involved in 
obtaining the 'structure stiffness matrix* is presented in 
this chapter, 

4,2 DISPIaCEMEMC EQRlflULATION 

In this formulation the displacements u at any 
point within the element are approximated by a set of functions 
and can be expressed in terms of nodal displacement 6 as 

H = f-M'} 


(4.1) 
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where , [h ] is called the interpolation or shape function. 

Having selected the displacement function, the state of strain 
within the element is given in terms of nodal displacements as 

.■ey = [B] Isi (4-. 2) 

where 6 is the vector of strain components at any point 
within the element and [b] is strain-displacement transformation 
matrix obtained by differentiating [H] . 

The stress within the element can now be expressed 
in terms of strains and material properties as 

[a] = [ 0 ] [el (4.3) 

where is the stress vector at any point within the 
element, and [c] is the constitutive (stress-strain trans- 
formation) matrix. 

The potential energy of a deformed solid is given by 
the sum of the internal strain energy and the potential of 
the body forces and surface tractipns. Thus, for a linearly 
elastic body, 

Potential energy = strain energy -work done by external forces 

X = I / as, (4.4) 

where | X | is the body force vector; 

l^pj- is the surface load vector prescribed over ; 
and V is the volume of the element. 
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Substituting the values of a , ^u|.from equations 

( 4 «l)j ( 4 . 2 ) and (4-3) into equation ( 4 . 4 ), potential energy 
of the body can be written as 

u = ■! ^ [Cj[h]^6'^ dV [6]'^[Srf \XjdV 

-f [df (pj dS^ (4.5) 

Applying the principle of minimum potential energy for the 
element we 0 brain 

\bf (j[Bf [0][3] Sl^l^ dY - |ii dY - ^ [wf .^pj dS^)=0 

( 4 . 6 ) 

Since the variation of nodal displacements jo 6 l are 
arbitrary, the expression in parenthesis must vanish. This 
gives the equilibrium equation for the element 

[K] U,] = |PL (+.7) 

where [K] is defined as the element stiffness matrix and P 
is the nodal load vector given by 

[k] = / CbF [G] [b] dY ( 4 . 8 ) 

Y 

and {P] = ) x} d/ + f [df ipt dS. (4.9) 

Y ^ - -- ^ 

Once the element stiffness matrix is found, computing the 

overall structure stiffness matrix, imposition of boundary 

conditions and solution of equations follows the standard 
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direct stiffness method of analysis ( 55 j 57 ). 

4.3 THh] IlhfiARLY VffiYING STRAIh TRIjU 43-U1AEI(1.S.T , ) jELEIvIUiTT 

A linearly varying strain triangular element which 
has complete quadratic displacement model and satisfies all 
the convergence requirements has been used in the present work. 


4.3.1 Element Stiffness Matrix 


A linearly varying strain triangular element (56, 58) ha 
six nodes out of which 5 are primary nodes at the three 
corners and 3 secondary nodes at the mid points of three sides. 
This element has 12 degrees of freedom for plane stress problems 
ae shown in Eig. 4.1. In natural coordinate system, three 
coordinates 1 -| , ^ 2 * ^3 used to define the location of a 
point but only two of these are independent. Their relation 
to cartesian coordinates is given by 



(4.10) 


where (X^,Y^), (X 2 >Y 2 ^ (^> 1 ^ 3 ) are the cartesian 

coordinates of primary nodal points and 1 ^, l 2 > ^3 
natural coordinates of a point inside the element which 
are defined as 

I. 


'1 


= A ’ ^2 


A. 


A 


A, 

^3 ^ A 


9 


(4.11) 
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Where is the area of the triangle formed by the Tertices j, 
K and point P and A is the total area of the triangle. 

The quadratic displacement field for this element is 
given by 

2 2 

u_ = a^ + 0tpX + a_y + a .xy + a^-x + a^y 

' ^ ^ ^ ^ (4.12) 

2 2 

V = + ttgX + a^y + a^Qxy + a^^x + 


where are generalized coordinates. 

Using natural coordinates displacement model can be 
written as 



(4.13) 


^2^3’ 4li^i*^] 


[6a} = [a, Ug] 

[64'' = [-^1 ’Tg v, Yg Tg] 


and Uj_’s and V^'s are nodal displacements. 

Differentiating the expression for u and v we get the strain 
displacement relation as 


) 

>4" 

[of 


[6u} 

i “ 

[o]' 

N^J 


(ot) 


= CB][6} 


(4.14) 
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where •[sjr'*’ = [ 0^, 6^,, Yxy^ strain vector at any 

point within the element. 


C ' 'I 

l®ii= 2i 


( 41 ^- 1 ) b^' 


‘^(41^-1) a.| 

(4X2-1) b 2 
(4L3-I) b^ 

n4i < 

/ ^ 

( } 3^2 

( 4 ^^*" 1 } 

l 4 (L 2 b^+I^b 2 ) 

I 

^2*^^ 2^3^ 


4 ( X ^ 3 ' 2 'f'X 2^^ 

4(l^b^+l^b^) 


_4(X^a^H-X^a^) 


(4.15) 


where 



*>1 = ^ 2-^5 

a2 = Xi - X^ , 

bg = Yj-T^ 

5 

b, = Y^-Yg 

,Y^), (X2,Y2), 

(X^,Y^) are 


(4.16) 


primary nodal points. 


The strain displacement matrix is not constant. Since 
the strain components ame linear function of natural coordinates, 
they can be written by linear interpolation using strains at 
three primary nodes. Therefore the nodal strain vector will 
be 



T 



®X2 ®Y1 


®T2 ^13 XY1 XY2 XY^ 


(4.17) 


where ®Y YxY- strains at primary 
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no de s . 


'Ihe strains at the primary nodes can be written in 
terms of nodal displacement by evaluating the matrix [b] of 
equation (4,14) at primary nodes. This gives 


("4 



[of 

[of 



[Bnl] 


IM 

N 


= (4.18) 


where 


m 1 

"5^ 

■ ^2 

~b^ 

4b2 

0 


3b2 



4b, 



-h2 

3b5 

0 

4b 


( 4 -. 19 ) 


3a.j 

^3*2 


43^2 

0 

4a^ 


-ai 



4a.| 


0 

(4.20) 


"**^2 

39'^ 

0 

4*3.2 

43.] 



-k 


and [ B„] is called nodal strain displacement matrix. 


Assuming that the material properties are constant with 
in the element, stresses will also -^ary linearly within the 
element. The nodal stress vector is defined as 

■i?n}^ » [ ^Z2 ^Y1 °Y2 ^13 '^rfl '^XI2 ^XY3 1 
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Since both stresses and strains vary linearly within 
the element, their interpolation models will be id.entical which 
are given by 


= tStaJicTn] (+-22) 

te) = [% J[®n} (4.. 23) 


where J errand ^6^ are stress and strain vector at any point 
within the element 


[ife]= = 


and 



[of 



I of 

W'" 

loj’^ 

(4.24) 

[of 

[of 




= 1.2 


(4.25) 


The stresses and strains are related through material 
constitutive matrix [ 0] as 

(erj- =[C]^'e}- ( 4 . 26 ) 


The nodal stresses and strains, using equation (4.26) can be 
written as 



c^aP] 


] G^j[I] 

O22P] Cajfi] 
*^25 P ] P ] 


{%\ = PnJijn] 


(4.27) 


9 



204 


where O-j-^'s are element of [ 0] xn equation (4.26) 

X J 


[I] 


IS a 3n3 identify matrix and 


To ] IS called the nodal stress strain transformation 
matrix. 

The strain energy U of element is given by 


1 


U = ~ {ol 6 .Y 


(4.28) 


V 


Substituting the values of and|aj from equation (4.22) 

and (4.23) in above equation we get 
gi T 




(4.29) 


Putting the value of | from equations (4* 18) and 

(4.27) in above equation strain energy of element can be 
written as 

/n '' TT ^ 


Therefore stiffness matrix [Kjcan be defined as 


T 


[K] = [B^] [D] [C„] [B^] 


( 4 . 31 ) 


where 


[D] = f [hpf [N^] d7 


V 


(4.32) 


Por 

an element 

of constant 

thickness. 

D is 


= t / 

A 



[0 ] 

[0 ] ■ 

[P]^ 

[N^] dA = 

[0] 

[I),] 

[0 ] : 




[0 ] 

[0] 

C ] 


(4.33) 
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where = t 


/ 

A 




dA 


Ah. 

12 


"2 1 

1 2 
" 1 1 



(4.34) 


A and t are the area and thickness of the element respectively. 


Ihus element stiffness matrix of linear strain triangular 
element is given by the expression (4.51) that is by direct 
multiplication of matrices [bJ , [D ] and [ Oj , 


The matrix [bJ is based on the displacement vector 
| 6 j in which all degrees of freedom corresponding to x directior 
are grouped first and then all degrees of freedom corresponding 
to Y direction. To simplify the assembly of complete stiffness 
matrix and solution coding, the vector of nodal displacements 
IS rearranged so as to have to components (u,v) at each node 
together. The modified displacement vector is as follows 



u, 


U2 ^2 ^3 


V 


3 


u 


4 


V 

4 


u 


5 


V 

5 


^6 



(4.35) 


The nodal displacement vector d]- and | 6 'i are related through 
a transformation matrix [T'] as 

[6 J = [ T 'J (4.36) 


where , 
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[T. ] = 


1 00000 0 00000 > 
001 000000000 
00001 0000000 
0000001 00000 
000000001000 
00000000001 0 
01 0000000000 


( 4 . 37 ) 


0001 00000000 
000001000000 
00000001 0000 
0000000001 00 
000000000001 


4 . 3.2 Consistent load Vector 

Consistent load vector for l.S.T, element (56y58) is 
calculated from equation (4-. 9). Eor uniform body forces, the 
contribution to load vector is 

[OOOOOOiyxyiy ] (4.38a) 

i j 5 

and for a linearly distributed load vector, the contribution 


to load vector is 
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+ 

(2) 

loP 

3^X1 



n ^3) 


n (2) 


+ 


l,p(3) 

i- 

l-iP^^ ^ 
1-^X2 

3^X2 



+ 

1 72 

3^72 






+ 

1 p^^^ 


2^73 


+ 

pi’b 


+ 

p(5)) 

(pig' 


4lh 

2^1 (Pi4' 


pC’) 

''h (pij' 


p(2b 

% 


+ 

p(2) \ 

7-] ^ 


(4.38b) 


where 1 




"bln 

is the length of i side of a triangle , that is 
the side opposite to node i, 

IS the load intensity per unit length at nodd i in 
X direction of a linearly varying load that has been 
applied to side j of the element and 
is the load intensity per unit length at node i in 
Y direction of a linearly varying load that has been 
applied to the side j of the element. 
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4.4 THB mmARLY YARYISG STE.AM BJffi BLB1IGBT(L.S,B. ) 

The steel reinforcement is idealized as an assembly 
of one -dimensional bar elements in the present work. Bond 
between reinforcement and concrete or between reinforcement 
and brickwork is assumed to be perfect upto failure. Since 
linear variation of strains have been used for brickwork or 
concrete elements the bar element should also have linearly 
varying strains to achieve strain and displacement compatibility. 

Fig. 4.2a shows the representation of one dimensional 
bar element with degrees of freedom in both local and global 
coordinate systems. The linearly varying strain bar element (56) 
has three nodes. It has a total of six degrees of freedom in 
global coordinate system and only three degrees of freedom in 
local coordinate system as shown in Fig, 4,2b, The axis 
of element may be at some orientation 9 relative to gldbal 
coordinate system. The local natural coordinate system for 
bar element is shown in Fig, 4.2c, The relationship between 
local cartesian coordinate and natural coordinate I is 
given by 

1 1 

X = 2 (1-I‘)x^ + ~ (1+1) X2 ( 4 . 39 ) 

which gives 

X- (x^+X2)/2 X - x^ 


I = 


( x 2 - x ^)/2 
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2 

or L = ~ (x-x^) ( 4 -. 40 ) 

where , X 2 an-d x^ axe the nodal cartesian coordinates in 
local caxtesian coordinate system and 1 is the length of the 
bar elaaent. The quadratic displacement function for three 
noded bar element is given by 

2 

u = a^ +a2X + a^x 

where a^'s are generalized coordinates. 


Using natural coordinates, the displacement model 
interms of interpolation functions can be expressed as 



(4.41) 


where 

= 1 1 (1-1) 

Ug = ^ 1 (1+1) (4.42) 

= (1-L^) 

and ■]6| = ['a-] u .2 u^] is the nodal displacement vector in 

local coordinate system. 


The strain displacement relation is obtained by 

differentiating the expression for u and can be written as 
du ^ 1 2 fU- 1 ' 

ax |^2 f 

u, 


(4.43) 



210 


or 


iej = [B] [6]^ 


The stress strain relationship is 


<j — E 6 
X s X 


ioj= [C]l3j- 


(4.45) 

(4.46) 


where is the modulus of elasticity of bar material that 
is reinforcement, Therefore, stiffness matrix of a bar element 
in local coordinate system can be written as 

iT 


[k] = / [b]-^ [g] [b] 

e Y 


dV 


Isfs p 
l2 1 


(2L-1)(2L~1) 
(2L-1)(2L+1) 
- 4L(2L-1) 


(4.47) 

(2L_1)(2L+1) -4L(21-1) 

(2I+1)(21+1) -4i(21+l) Idl 


- 41 ( 21 + 1 ) 


161^ 


(4.48) 

where is the cross sectional area of the element. Integration 
of equation ( 4 . 48 ) gives the stiffness matrix for linear strain 
bar element as 


[id 


A B 

__.s_ _s 

31 ^ 


7 

1 

-8 


1 

7 

-8 


-8 

-8 

16 


( 4 . 49 ) 


The element stiffness matrix from local coordinate 
system is to be transferred to global coordinate system. Using 
the basic displacement transformation|^ 6 j 2 _ == fe-j] bo relate 
the displacements in two coordinate system, the desired stiffness 
matrix transformation is obtained as 
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[K^ = [K]g ( 4 - 50 ) 

wiiere [K ] is the ej.emeiit stiffness matrix in global coordinate 
system and 

[1-, ] is the transformation matrix which transforms the 
local displacement vector to global displacement 
vec'cor. 


]?or linear strain bar element, the transformation 
matrix is given by 



GosG 

SinO 

0 

0 

0 

0 ~ 


0 

0 

Cos9 

Sin9 

0 

0 

(4.51) 

0 

0 

0 

0 

Oos0 

Sin9 



where 9 is the inclination of local coordinate system to global 
coordinate system. 


4.5 PSEUBO LOADS DUE TO GRAdnNG OR YIElDIh'G 


The total pc jential energy of an elastic uncracked 
concrete or brick element is given by expression 

U = J + f (4.52) 

where U is the strain energy of concrete or brickwork 
W IS the potential energy due to external loads. 

When this element cracks or yields, material constitutive 
matrices will change and as a result of this the strain 
energy will be released. If is the strain energy released 
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due to crocking and due yielding of concrete or 
brickwork, tke expression for potential', energy can be 
modified CIO (od, 33) 

* = - °oo - V 

For cinite element displacement formulation, the potential 
energy er^uation ( 4 . 53 ) can be written as follows 

+/ [B]^[C^y]p 3 ] dy^y)j 6 }-j 6 }^|P} (4.54) 

^cy 

wherefO^ ] and [C„_.] are the change in material constitutive 

0 (2 Oj 

matrices due to cracking and yielding of concrete or brickwork 
respectively^ 

[ K ] is the stiffness matrix of concrete or brickwork element 

O 

[P ] IS the externally applied load vector which includes 
body force, surface traction and concentrated nodal loads. 


Minimizing the potential energy, we obtain. 


where and l^cyl called pseudo loads and are defined as 

( f ~ pseudo load vector due to cracking 

^ ^QQ ' of concrete or brickwork 

{ P !=(/ [B]'^[doJC®] '^^Gv “ pseudo load vector due to yielding 


of concrete or brickwork. 



213 


Similar 1;' ,we can derive the expression for pseudo loads 
due to yielding of steel reinforcement. The expression will 
be as follows 

ihyl = ( + -56) 

sy 

where [G_J IS the change in material constitutive matrix due to 
oy 

yielding of steel reinforcement. 

The pseudo loads in equation (4.55) are due to the 
stresses released in the element because of change in material 
properties. The pseudo loads are computed as follows. 

Considering the pseiado load due to cracking 

1^4 = 4B]ho,,][B]dV,48} 

^CC 

If [c j IS the original constitutive matrix of uncracked and 
[OcJ is that of cracked element, then ^cc ” ^ ” ^cr 

{ ^oi= ,{ [3]'^ [0] - [OJ [Bjj6} 

c c 

= [ 0 ] - 4 /^^ 

ho 

= ■^ho < + 

^CC 

where I^To} stress released due to cracking. 

Similarly pseudo loads due to yieldin^g of concrete or brickwork 
can be computed in terms of stresses released. This only 
requires modified material constitutive matrix. 
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4*5.1 Ooaputation of Pseudo Loads on Goncrete/Brickwork 
Element 


In bhe computational procedure, each element is 
checked for cracking and yielding at various stages of 
loading. One of the major advantages of using I.S.T, element 
is that quite accurate results for linearly elastic analysis 
of flexural members can be obtained with relatively coarser 
mesh as shown by Jain, A.K. (41 )• However the computation 
of pseudo load vector becomes slightly domplicated due to the 
fact that stresses are not constant over the element and 
coarser mesh may lead to less accurate results for nonlinear 
analysis if cracking, yielding etc. are checked at 'one point 
(e,g. centroid) only. To overcome this difficulty, the 
element is further ^vided into 4 subregions (programme 

'Y’\ " 

developed accomodates even 16 subregions) as shown in Pig. 4. 3 
for the purpose of pseudo load computations. The strains and 
stresses are computad at the centroid of each subregion to 
check for yielding and cracking. The pseudo load vector as 
given by equation ( 4 , 57 ) for the subregion may be written as 


where|A?gJ- is the pseudo load vector for the subregion, 

IS the volume of subregion that has cracked or 
yielded and 

is 'the vector of stress released in the subregion. 


[a do} 
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i’Tow for ii.S.f . elemeni;, it ca.n be written as 

' ’'r 

= S [Hef iA'^o) 

^r 

On integration this gives 


[ABsj= ^ '‘r [®J’' 


f\^X 

ho 


ho 


ho 

dPj 

ho 


ho 


ho 


ho 

l^XJ 

ho 

m 

^30 


(4.58) 


wliere is the area of subregion, and are 

components of stress released vector^^Ql , and and 

are natural coordinates at the centroid of subregion, 

4 , 5,2 Computation of Pseudo load Vector on Steel Reinforcement 

For the computation of pseudo loads, the steel 
reinforcement is divided into one or four subregions as shown 
in Fig, 4.4 and yielding is checked at the centroid of each 
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subregion. The pseudo load vector on the bar element 
subregion would be 

[zvhj = .[ “ 

where [B ] is the strain-displacement matrix for the bar element, 
and stress released in the subregion due to 

yielding. 

If bar is divided in one subregion as shown in 
Big, 4«4-a the pseudo load vector after integration is 

1 - 1 - 1 ' °} (+- 59 ) 


If bar is divided into four subregions as shown in 
Big. 4.4b, the pseudo load vector will be as follows. 

Pseudo load vector for yielding of subregion 1,2,3 and 4 are 




T 




■/T 




= 1- 


s/Vo 


5 
1 T 


L 2 

J -^s^o 
i. 2 

^ 2 


L 

4 

3 

4 

1 

4 

1 

4 


1 

4 

1 

4 

3 

4 

5 
4 


3 

2 

1 

2 


} 


3 

2 


(4.60) 


These pseudo load vectors are in local coordinate system. 
These are transformed to global coordinate system using the 
transf ormafcion matrix [T^Jof section 4 . 4 . Using transformation 



matrix, tLa pseudo load vector in global coordinate sys-tem 
wo uld be 



4.6 iTOxTLIJUAR ijfiiliYSIS 

'-Chcre are two sources of nonlineariby in structural 
problems, first is the material nonlinearity which results 
from nonlinear constitutive laws. Second type of nonlinearity 
arises from the violation of linear strain displacement 
relationship and is called the geometric nonlinearity. The 
geometric nonlinearity has been neglected in the present work 
as the displacements are observed to be small. Besides material 
nonlinearity, the maximum contribution to the nonlinear behaviour 
of brick masonry supported on reinforced concrete beams is due 
to cracking of concrete and brickwork. The way to treat such a 
situation in the finite element method is to physically define 
individual cracks as they occur by successively modifying the 
elements and connectivities of nodal points. But for non- 
linear analysis the modification of finite element mesh 
involving renumbering of nodes as the crack propagates has been 
found to be q.uite difficult and computationally uneconomical. 
Therefore, the cracking has been incorporated for the nonlinear 
analysis, in the present work, by changing the material 
properties of elements and keeping the topology unaltered. 
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21iere are three mam approa-ches to carry out non- 
linear analysis of structures, 

1 . Incremental methods 

2. Iterative methods 

3. Incremental iterative methods. 

All the three methods are well documented in 

standard text boohs (55, 56) and therefore are not discussed 
here. Incremental iterative methods are useful and accurate 
when complete load response of the structure is desired and 
therefore, this method is used in present work. For 
completeness of presentation it is briefly described in the 
subsequent section, 

4,6,1 Incremental-Iterative Methods 

Ihe incremental -iterative methods uses a combination 
of incremental and iterative schemes. In this procedure load 
is applied incrementally, but after each increment successive 
iterations are performed until the solution converges. The 
two incremental-iterative methods that are most commonly used 
in literature are: (i) incremental tangent stiffness method 
which updates stiffnesses, strains and stresses each time, 
and (ii) 'initial stress' method which uses the original 
elastic structure stiffness matrix in all load steps and 
updates only strains and stresses in each iteration. In 
both the methods a corrective force is applied at each step 
to restore the structural equilibrium. 
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4. 6. 1.1 Incremental tangent stiffness method 

'Ihe essence of this method is that at any stage 
a nodal force system equivalent to the total stress level is 
evaluated and compared with the prevailing applied loading 
systems. The difference between the two results is a set 
of residual forces that can be interpreted as a measure of 
any lack of equilibrium. To restore equilibrium the residuals 
are then applied to the structure and problem is resolved 
after updating the structure stiffness matrix. This process 
is repeated until residuals are sufficiently small before 
applying the further load increments. A typical step of the 
method is described with the aid of Fig. 4-. 5 which is a 
single degree of freedom problem. The point B corresponds 
to the start of the increment i, that is point B is at the 
end of the increment (i-l). The solution by using the 
tangent stiffness corresponding to point B leads to the 
point B, V 7 hereas the tame solution is at point 0, The 
total deflection corresponding to the point D in multi- 
degrees of freedom system is given by 

= (4.62) 

The total load that is equilibrated for j Q indicated by 
point D' is 

= [h-i} 


(4.63) 
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where D h' represents denoting the unbalanced load 

for the iteration 1 of load increment i. In general, for 

th 

3 iteration of i increment, the total equilibrated 


load will be 


Pi) = -Ia^o4 


(4«64-) 


Referring to I*ig. 4»5b, the total strains in the element 
corresponding lo structure displacement is given by 

[®]l> = i-^^i} (4*65) 


Use the tangent stiffness corresponding to the point B gives 
the stress value at B (Big. 4«5b) which is higher than the 
actual stress at D' equilibrated by the strains Sj) • Therefore, 
the total stress that is equilibrated for 6 indicated by 
point B* is 




}b' = 


( 4 , 66 ) 


where|^/\0' ^ j- is the correction stress for iteration 1 

th 


of load increment i. In general for the j 
total equilibrated stress will be 


increment the 



(4.67) 


It can be seen from equations (4.64) and (4.67) 

(Bigs, 4.5a and 4.5b) that the values converge to true solution 
when correction stress£ikV 7 )jand correction loadl^P^^ tend to 


zero 
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'fhe total stress increment far j iteration of 
"blo. 

the 1 load increment is expressed as 

= [C] (+-63) 

where [C ] is the constitutive matrix corresponding to 
the current state of stress. But the determination of 
current state of stress Is not possible without knowing the 
corresponding constitutive relation. To overcome this 
difficulty, the solution is found by approximating the 
material properties corresponding to the stress state of 
the previous cycle. Now writing the correct stress 
increment as 

Hi] = tO'] {Ae|(- - (4.69) 

where tc] is the material constitutive matrix corresponding 
to the previous cycle, the potential energy increment can, be 
written as 


A” = i {AQ} - {A^}* 

(4.70) 

or the equilibrium equation is given by 

[Kj] |a 6} ={_APj +[APoi- 

"tlx "tlx 

or for j iteration of i increment, equilibrium equation 
is as follows 

+ [APy 


(4.71) 
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where [2^ ] = [ ^ D > Ij 

= 3 = 1, 

" 1^1 3 > 1 


(4.72) 

( 4 . 73 ) 



eLments''^^ L^°ol^l 


d > 1 
D = 1 


( 4 . 74 ) 


The incremental solution of equation (4,18) is 
obtained by computing the correction load[A^Qj[_^ for each 
iteration, updating the structure stiffness matriz [Krji] 
due to change in material property and iterating the process 
to convergence. 


' The incremental -iterative method with variable 

stiffness requires the computation of the structure 
stiffness matriz at each step of solution procedure. This 
suffers from an economic disadvantage because a complete 
reformulation of the stiffness matrix and a now solution 
of governing equations are required at each iteration. 

This difficulty can be overcome by adopting a method v^hich 
uses the same stiffness matrix repeatedly in the solution 
procedure as described below. 


4*6, 1,2 Initial Stress Method 

As correction loads in equation (4.71 ) serve 

to satisfy the equilibrium, it is not essential that the 
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tangent stiffness matrix be used. Instead, the tangent stiff- 
ness matrix in equation (4.71) can be replaced by the initial 
stiffness matrix [Kj corresponding to the original elastic 
stiffness of the structure. This form of the solution has 
been called 'initial stress' method (26,27,55)» Numerically 
the 'initial stress' method has a unique advantage of using 
the same stiffness matrix at every stage of iteration. 

After the first iteration, each subsequent iteration can 
be performed in a small fraction of time needed for the 
first solution since the inverted stiffness matrix is 
available. The use of the 'initial stress’ method where 
the structure stiffness matrix is computed and inverted 
only once in the solution procedure has been found computa- 
tionally advantageous (27). Therefore this method is 
adopted in the present work. 

In each iteration of this method, the difference 
between the true stress level corresponding to the appropriate 
strains and that corresponding to an elastic solution is 
determined. This stress difference is redistributed elastically 
to restore equilibrium. 

Fig, 4.6 shows the essential steps of the 'initial 
stress' method. The explanation of Fig. 4*6 follows exactly 
that of Fig. 4.5 except the difference that the line ED 
remains parallel to the initial tangent. When the original 
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material properties are used, then potential energy increment 
of equation (4,70) gets modified to 

Alt K„{AQ} h{A?P (+. 75 ) 

or the equilibrium equation is given by 

(+•■ 76 ) 

where [K^] is the initial tangent stiffness matrix which 
is given by 

[K ] = / [sf [0 ] [b] dY (4.77) 

V V 

where [C^] being the initial material constitutive matrix. 

The values of nodal load vectors|/\P? ^ and^fvP^^ are given 
by equations ( 4 . 73 ) and (4,74). 

4.7 OOFTERG-MGE ORIThlllA 

In the incremental-iterative procedure based jon 
’initial stress’ method used in the present work, the 
equilibrium equations are solved for nodal displacements. 

The two most obvious criteria of measuring convergence at 
the end of an iteration are the magnitude of forces by which 
equilibrium is violated or the accuracy of the nodal displace- 
ments* The violation of equilibrium is measured by the 
magnitude of the residual unbalanced nodal forces. The 
accuracy of the nodal displacements can be measured by the 
magnitude of the additional increments of displacements. 
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In the present investigation increments of nodal 
displacements and residual unbalanced nodal forces both 
are considered to check the convergence, fhe maximum vector 
norm is used to measure the error for each component of 
residual displacement and unbalanced force. The ercors are 

1 1 Sjyl 1 = max 1 I Q^l I for displacement vector; and 

I 1 Epl 1 = max j 1 i unbalanced force vector, 

n 

When both the errors become smaller than their 
convergence criteria, the program will stop the iteration 
and go on to the next load increment. The residual imbalan- 
ced nodal forces will be carried over and added to the 
next load increment. 

The tolerance for convergence could be decided on 
the basis of the required accuracy of the results which is 
often dictated by the accuracy of measurement. In the 
present v/ork the tolerence on unbalanced force has been 

taken as 1.0 kg. while the same on displacement has been 

—6 

taken as 1x10” cm, 

4.8 EESUITS OF EINITE ELEi'lEET iUIkLYSIS 

The computer programme developed in Fortran lY 
in the present work based on foregoing description has been 
used to analyse the specimens investigated and parametrically 
studied experimentally. The analytical results are obtained 
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on DEG system 1090 at III' Kanpur and these are compared 
with the esperimental results reported in Chapter III. 

Ihe composite structures and load thereon is simulated to 
represent identical system as investigated experimentally. 
Half symmetry of the system has been exploited all through 
except where unsymmetric openings existed. 

Experimental investigation has revealed that the 
composite system behaves linearly in the initial stages 
until the load approaches the first crack load. This 
has been verified analytically as is shown in Pig. 4«15 
and Pig. 4»15. Hence, inorder to carry out nonlinear 
finite element study of the composite system, starting 
load does not have to be zero load. On the contrary the 
starting load can be oust below the first crack load for 
each specimen. This of course is not apriori known. 

Therefore, a numerical experimentation was carried out 
in each case to decide upon the initial load which was 
almost 3 tonnes below the first crack load for each specimen. 
This saved the computational effort as well which would have 
been otherwise required if the load on the system would have 
been incremented starting from no load condition. The load 
was incremented in each case by 3 tonnes in the present 
analytical study. Por each load increment an upper limit 
of 300 on the number of iterations was imposed to allow 
the pseudo-loads and deflections to converge. While at 
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tile failtore load, pseudo-loads converged in all cases within 
this I^-imit on the rimber of iterations the corresponding 
deflections failed to converge in 3OO iterations. Increasing 
the upper limit on the number of iterations to ensure 
convergence of deflections at higher loads near the failure 
would have resulted in prohibitive use of computer time. 
Therefore, this was not done in the present work. The 
termination of the programme at 300 iterations or the 
crushing of concrete whichever occurred earlier was taken 
to be the on set of failure of the composite system in the 
present analytical study. 

4.8.1 Effect of Variation in Cement Sand Mortar for Brickwork 

Specimen number 1 to 10 (described in the previous 
chapter) are S3mimetric about a vertical centre line passing 
through mid span. Therefore, half of the specimen size was 
modelled by 42 I.S.T, elements for brickwork and concrete and 
21 L.S.B, elements for bending reinforcement and vertical 
connectors as shown in fig, 4.7. 

figs, 4.8(a) to 4.12(a) respectively show the plot 
of longitudinal stress, cr , at various cross sections for the 
specimens number 1 to 5 (cast in mortars of 1:3, 1 !4) 1^5, 

1:6 and 1:8 respectively) at the initial load, from these 
plots it is seen that the supporting R,C. beam is always 
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under te.-sion except very ne^r the supports. The compressive 
stresses are borne by the brickwork. The vertical stresses, 

^ and the shear stresses, have been plotted in Pigs* 

4.8(b) to 4.12(b) and 4.8(c) to 4.12(c) respectively for 
the five specimens analysed at the initial load. Prom these 
figures it is observed that vertical stress and shear stress 
concentrate in the R.C. beam and brickwork 3ust above it 
very near the supports. Pigs. 4.8(d) to 4.12(d) show the 
stress along the bending reinforcement at various loads 
starting from initial load upto failure. The bending 
reinforcement in four of the five specimens (i.e, all except 
one cast in mortar ratio of 1 ;8) reaches thd yield stress 
at failure. The stress in bending reinforcement increases 
rapidly as the concrete cracks. Pigs. 4.8(e) to 4, 12(e) 
show the vertical deflection at the bottom of specimen along 
the span at various loads starting from initial load upto 
failure. Prom these plots it is seen that the deflection 
shape has point of inflexion near the supports meaning 
thereby that the specimens deflect as simply supported beams 
with small over hang on either side. This is in accordance 
with the plots shown in Pigs, 4.8(a) to 4.12(a) where the 
nature of longitudinal stresses change from near the supports 
to the mid span. Pig. 4,13 shows the plot of load versus 
deflection at mid span obtained analytically. The corresponding 
experimental deflections obtained are also shown on this 
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figure for the sake of comparison. Prom these plots it can 
be concluded that experimental and analytical results of 
load versus deflection compare very well all through upto 
failure load. Pig. 4.14 gives the plot of longitudinal 
strains for various specimens at the mid span cross section 
under the initial load. The variation of longitudinal 
strains is linear as observed experimentally. Table 4«1 
gives the first crack load and the failure load for specimen 
1 to 5 obtained analytically. The first crack load obtained 
analytically for specimen number 1 cast in 1 :5 mortar is 
about 25 percent lower than the corresponding experimental 
value. However, for specimens number 2 to 5 cast in mortar 
1:4 to 1:6 and 1:8 respectively these results are approximately 
10 percent lower than corresponding values obtained experimentally 
as given in Table 5*16. The general reason for this discrepancy 
is probably the fact that specimens investigated experimentally 
actually cracked earlier than observed. However the large 
difference in case of specimen number 1 is only attributable 
to variation in material properties and workmanship of the 
experimental specimen. The analytical values of failure load 
IS lower by 10 to 20 percent in case of specimens number 1 to 4 
cast in mortar 1:3 to 1:6, while for specimen number 5, cast 
in mortar 1:8 it is higher by 50 percent as compared to 
corresponding experimental value of failure load. 
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The failure load anaj.ytically obtained corresponds 
to either crashing of concrete or termination of the programme 
for each load increment at 300 iterations, whichever is 
earlier. Por the specimens number 1 to 5, the failure load 
obtained did not result into the crushing of the concrete. 

In other Virords, it is corresponding to the premature termina- 
tion of the programrae only due to the limit imposed on the 
number of iterations without attaining the convergence on 
deflections. This, however, was unavoidable in the present 
work because of prohibitively high computational time which 
would have been otherwise required. It is this premature 
termination of the programme which has resulted in the 
lower failure load obtained analytically for specimens 
number 1 to 4 as compared to the failure load recorded 
experimentally. However, it would be recalled that during 
experimental investigation, the dial gauges were removed to 
ensure safety of tne devices and the observers, at a load 
where the deflections had significantly increased than 
those in the previous load increment and the failure was 
expected to be Imminent. After the removal of the dial 
gauges the experimental specimens were further loaded to 
observe the actual failure load. Significantly, the analytical 
failure loads obtained are comparing very well with those at 
which the dial gauges were removed in the experimental 
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investigation as shown in Table 4»2. Prom this it can be 
safely concluded that the analytical method used in the 
present work is correctly predicting the imminent failure 
load. Furthermore , the analytical result so obtained is 
conservative and hence can be safely used for design 
purposes. Specimen number 5 cis.t in 1:8 mortar give a 
higher failure load analytically as compared to the correspon- 
ding experimental value. The reason for this, as discussed 
in the previous chapter, is the load failure between the 
connectors and the brickwork in 1 :8 mortar when investigated 
experimentally while in analytical investigation a perfect 
bond has been assumed between the connectors and the brick- 
work, Naturally, therefore, the latter value turns out to 
be higher. 

Pigs, 4.8(f) to 4.12(f} show the crack pattern 
obtained at failure for the five specimens. These crack 
patterns have been plotted from the computer print out 
taking into account the information regarding the subregion 
of the element cracked and the direction of this cracking at 
the centroid of the subregion. Anolyt^ical plot of crack 
patterns show^ more number of cracks in a specimen as 
compared to those obtained experimentally. The reason for 
this is that in the anilytical work the finite element 
subregion in which the tensile stress exceeds the .allowable 
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limit is assumed to have cracked and a crack is shown to 
initiate from that particular element or subregion perpen- 
'cular to the direction of the principal tensile stress. 
Experimentally, however, number of cracks do not develoo 
but the crack already formed initially goes on v/idening. 

In any case, the nature of crack obtained both experimentally 
and analytically is principally the same. 

The above analytical study firmly reinforces the 
conclusions drawn after the experimental investigation in 
this regard. 

Table 4.3 gives the initial load applied, number of 
load increments and computer time taken for the analytidal 
study of specimens number 1 to 5. 

Pig. 4.15 shows the analytical load versus deflection 
curve under tensile loading for the specimens number 6 to 10 
(cast in mortar 1 : 3 , 1 : 4 , 1 : 5 , 1:6 and 1:8 respectively) at 
the mid span. These results are under the assumption that 
perfect bond between the connectors and brickwork exists. 

The corresponding experimental curves are also shown on the 
same figure for purposes of comparison. Prom these curves 
it is observed that deflections obtained experimentally are 
much higher than those obtained analytically. It is also 
observed that cracks obtained analytically are linear upto 
the first crack load while experimental curves show non- 
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linearity right from zero load. The discrepancy in the 
analytical and experimental results is due to the fact 
while experimentally bond slip between the vertical connectors 
and the brickwork is established, analytically perfect bond 
between the two has been assumed . Therefore, to obtain 
a correct analytical load deflection response, introduction 
of bond link elements is necessary. This however, has not 
been attempted in the present work and any further analytical 
study for specimens under tensile loading has been abondoned^^ 

4.8.2 Iffect of Variation in Height to Span Ratio 

Specimens number 11,4,12,13 and 14 having height to 
span ratios of 0.25, 0.33, 0.4, 0.5 and 0,8 are analysed 
under compressive loading by the computer programme developed 
in the present work and the results so obtained are compared 
with the corresponding experimental values, ill these 
specimens being symmetrical about the vertical line passing 
through the mid span, half S3mimetry has been exploited in 
the analytical study. Specimens number 11, 12 and 13 were 
modelled by 42 L.S.T, elements for brickwork and concrete and 
21 L.S.B. elements for bending reinforcement and connectors. 
Specimen number 14 having a height to span ratio of 0.8 was, 
however, modelled by 56 L.S.T, elements for brickwork and 
concrete and 28 l.S.B. elements for bending reinforcement 
and vertical connectors. 
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Fxgs. 4.16(a) to 4.19(a) respectively slaow th.e 
plot of longitudinal stress, cr , at the initial load at 
various cross sections for the specimen number 11 to 14« 
iProm these plots and from Pig. 4.8(a) which is fo'r the 
longitudinal stress , cr^, for a height to span ratio of 
0,33 it is seen that the supporting thin R.O. beam is alv^/’ays 
in tension except very near the supports as was the case 
for specimens number 1 to 5. It is further observed that 
the variation of longitudinal stress, m brickwork is 
almost linear upto a height to span ratio of 0.5. Por 
specimen number 14j having a height to span ratio of 0.8 
the longitudinal stress at various cross sections are very 
small and these do not vary linearly. Infact the brickwork 
is relatively free from the longitudinal stress in all 
specimens investigated. It is further observed from these 
plots that the neutral axis at initial load (i.e. for 
uncracked section) is at a depth almost equal to 2/3rd the 
height of ccmposite system. Vertical stresses, a » and the 
shear stresses, have been plotted in Pigs. 4.16(b) to 

4.19(b) and 4.16(c) to 4.19(c) respectively for the four 
specimens number 11 to 1 4 , It is observed from these plots 
that concentration of vertical stress and shear stress takes 

place near the supports in the R.C. beam and the brickwork 

¥ 

just above it. Pigs. 4.16(d) to 4.19(d) gives the stress in 


steel along the bending i^einf orcement at various loads starting 
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from initial load upto failure, it failure, the suress 
in bending reinforcement in all specimens reaches the 
yield stress. Pigs. 4.16(e) to 4.19(e) show the vertical 
deflection at the bottom of composite system along the 
span at various loads starting from initial load upto 
failure. Prom these curves it is seen that hxgher is the 
height to span ratio lower are the deflections for the 
same load. Pig'. 4.20 shoves the plot of load versus deflection 
at mid span obtained analytically. The corresponding 
deflections obtained experiiaentally at various loads are 
also shown on this figure for purposes of comparison. Prom 
these curves it can be concluded that the experimental and 
analytical results of load versus deflection compare very 
well upto the failure load except for specimen number 11, 
which has a height to span ratio 0.25. In case of specimen 
niMber 1 1 , the experimental and analytical load deflection 
curves upto first crack load overlap each other but beyond 
first crack load experimental deflections are much higher 
than obtained analytically. The reason for this difference 
is that, the length of connectors is not sufficient , to 
transfer the load and therefore bond slip starts after the 
first crack load and this bond slip causes higher deflections 
and subsequent early failure of the specimen. However, 
since perfect bond has been assumed in the analytical work, 
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analytical results show lower deflections and higher failure 
loads, fable 4.4 gives the first crack load and the failure 
load for specimens number 4 and 11 to I 4 having various height 
to span ratios. The first crack load for all specimens are 
slightly lower than observed experimentally for the correspond- 
ing specimens as given in Table 3.53 except in the case of 
specimen number 12 having a height to span ratio of O .4 
where the analytical value is around 20 percent lower than 
the experimental value. This difference in t^o values for 
specimen number 12 is only attributable to variation in 
material properties and worknanship of the experimental 
specimen. The failure load obtained anaijrfcically is about 

10 percent lower than the corresponding failure load obtained 
experimentally except in case of specimen number 11 where 
experimental failure load is lower than the analytical value 
for reasons explained earlier. The reason for the lower 
analytical failure load is again the premature termination 
of computer programme because of the limit imposed on the 
number of iterations for a load increment as discussed in the 
preceeding section. Figs. 4.1 6(f) to 4.19(f) show the 
crack patterns obtained for the four specimens. The nature 
of cracks is basically the same as obtained experimentally. 
However, analytical plot of crack pattern show more number 
of cracks in a specimen as compared to th.ose obtained 
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experimentally for reasons explained in proceeding section. 

This study reinforces the conclusions drawn from 
the experimental work that higher is the height to span 
ratio, higher is its load carrying capacity, Furthermore, 
the analytical results for the first crack load as well 
as the failure load are slightly conserrative from designers 
point of view. 

Table 4.5 gives the initial load applied, number 
of load increments and computer time taken for the analytical 
study of specimen number 11 to 14. 

4.8,3 Effect of Variation in Size and location of Openings 

Specimens number 19 to 23 having openings as detailed 
in the table on next page , have been analytically analysed by the 
computer programme developed in the present work. Half 
symmetry has been exploited wherever possible. The table 
also gives the details of finite element model for each 
specimen. 

Figs. 4.21(a) to 4, 24 (a) gives respectively the 
plot of longitudinal stress, o , at the initial load at 
various cross sections for specimen number 19 to 22. From 
these plots it is seen that behaviour of the composite system 
with symmetric door and window opening as well the unsymmetric 
window opening is more or less same as that of composite wall 
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DliTAlLS OF FIIIT3 FLBMBMI’ .iOFEL CF SPECniElirS 


Y/ITH OPENINGS 


Speci- 

men 

number 

Type 

Of 

Location 

No .of 
nodes 

Number of 
elements 

L.S.T. 

Number 

elements 

l.S.T 

opening 



In 

concre- 

te 

In 

brick- 

work 

In bend- 
ing 

re inf or- 
cement 

In 

conne 

ctors 

19 

Window 

Sjnrmetric 

161 

20 

46 

10 

23 

20 

Door 

Symmetric 

157 

20 

42 

10 

21 

21 

Window 

Unsymme- 

tric 

31 2 

40 

92 

20 

47 

22 

Door 

Unsymme- 

tric 

8 cm thick 
R.C.beam 

298 

40 

84 

20 

45 

23 

Door 

TJnsymme- 

tric 

556 

68 

84 

20 

45 


1 6 cm thick 
R.C.beam 
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without opening except at the lintel level. Refering to 

Rigs, 4.21(a) and 4.22(a) we observe that the lintel is 

under longitudinal compression and the brickwork above it 

is almost free of longitudinal stress. On the contrary we 

see in Rig. 4.23(a) that the lintel is under longitudinal 

tension and the brickwork over it is in, .longitudinal compression. 

The reason for this is the fact that the load through the 

brickwork above the unsymmetric window opening is directly 

transferred on to the lintel while for symmetric openings the 

/ 

load gets transferred through arching action. Conmfting to 
Rig, 4.24(a) we observe that the supporting R.C. beaa 
experiences longitudinal tension at the centre o' the opening. 
This is easily explained from Rig. 4.24(e) where the corres- 
ponding deflected shape shows hogging moment at this cross 
section. The behaviour at the lintel level and above in this 
case is identical to the one discussed for the case of 
uns3rmmetric window opening, Vertical stress, o . is plotted 
along the span at various heights of the wall for specimens 
number 19 to 22 and is shown in Rigs. 4.21(b) to 4.24(b) 
respectively. Concentration of vertical stress takes place 
near che supports and brickwork above it as observed in 
corresponding solid composite system and also at the ends 
of the R,G, lintel provided at the top of the openings. In 
case of symmetric door and window openings, the magnitude 
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of vertical stress at the supports and brickwork above it 
IS same at various heights under identical load. Moreover, 
the values also compare well at the corresponding heights 
of the composite system without openings under identical 
load. This is as anticipated since load transference in 
either case is through arching action meaning thereby that 
symmetric openings about the vertical centre line of the 
composite system do not affect the over all behaviour of the 
composite structure, kigs. 4 . 29 (c) to 4 . 24 (c) show respectively 
the plot of shear stress, at various cross sections for 

specimens number 19 to 22, The concentration of shear 
stress in all cases is incidentally at locations where the 
concentration of vertical stress occurs viz. at the supports 
in R.O, beam and at the two ends of R.C. lintel above the 
openings. Pigs, 4*21 (d) to 4 . 24 (d) show^ the plot of stress 
in steel along bending reinforcement at various loads starting 
from intial load to failure load. In specimens number 19 
and 20 having symmetric windovf and door openings bending 
reinforcement reaches the yield stress at the failure load 
but for the remaining two specimens number 21 and 22 having 
unsymmetric openings, the stress in bending reinforcement is 
below yield value even at the failure load. Pigs. 4.21(e) to 
4 , 24 (e) show the vertical deflection at the bottom of specimens 
all along the span at various loads starting from initial to 
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failure load. Plots of Pigs, 4.21(e) and 4.22(e) are almost 
identical to the corresponding plot of composite system 
without opening. The slight difference in the values of 
deflections at various points under the identical load is 
qualitatively accounted for from the fact that the stiffness 
of composite system due to openings reduces. Prom Pig. 4 . 24 (e) 
it is observed that there is a marked increace in the deflections 
at various points along the span from the initial to the 
failure load. 

Pig. 4*25 shows the plot of load versus deflection 
at mid span of all the specimens obtained analytically. The 
corresponding deflections obtained experimentally are also 
shown on this figure for purposes of comparison. Prom 
these curves it is concluded that the analytical and experi- 
mental results of load versus deflections compare well upto 
failure for all specimens except number 22 for which experi- 
mental load deflection curve is not available. Since it 
failed under its self weight. Table 4.6 gives the first crack 
load and the failure load for specimens number 19 to 23 . 
Analytical value of first crack load are slightly lower than 
those obtained experimentally (given in Table 3.47) in case 
of specimens number 19, 20, 21 and 23. Analytical value of 
failure loads is lower by about 10 percent to 15 percent for 
all specimens except number 22 for which no experimental data 
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IS at hand. The reason for the conservative analytical 
failure load has already been discussed. Fig. 4.21(f) to 
4.24(f) show^ the crack pattern obtained from analytical 
results. The crack pattern is ba,sically the same as obiainted 
experimentally. 

From this study it is concluded that symmetric openings 
do not appreciably affect the load carrying capacity of the 
composite system. However the unsymmetric window opening 
affects the load carrying capacity significantly while 
unsymmetric door opening affects the load carrying capacity 
adversly. It would be recalled that same conclusions were 
drawn from experimental investigations. 
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TABIE 4.1 


EIRST CilACK MD FAILURE LOAD 


H/L LOADING 
^.T3 GOMPEESSITB 


Specimen 

Number 

Mortar 

1st Ora- 
cle load 
in tonnes 

Failure 

load 

in tonnes 

Remarks 

1 

1:3 

15.0 

27.0 

Bending reinforcement at 
mid span has yielded. The 
first crack was at mid span 

2 

1 :4 

15.0 

24.0 

-do- 

4 

1 :5 

12.0 

24.0 

Bending reinforcement at 
mid span and brickwork at 
top has yielded. The first 
crack was at mid span. 

5 

1 :6 

1 2.0 

21 .0 

-do- 

6 

1 :8 

9.0 

15.0 

The first crack in concrete 
was in the mid span region 
while that in brickwork was 
near the support. The ultims 
failure took place because 
of yielding of brickwork 


HA LOiUING 

TABLE 4.2 ; HffilFEiJI FAILURE LOADS 

0.33 COMPRESSIYE 


Specimen Ilortar Analytical Experimental- 

number failure load imminent 

in tonnes failure load 

in tonnes 


1 

1 :3 

27.0 

28.91 

2 

1 ;4 

24.0 

24.78 

3 

1 :5 

24.0 

23.67 

4 

1 :6 

21.0 

22.71 

5 

1 :8 

15.0 

8.26 
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TABLS 4.3 


OOLiPUTER TIliU 


E/h 

104DIIG 

0.33 

GOlimBSSlW 


Specimen Mortar Initial Humber of Computer time 
number load in subsequent 

tonnes load 

increment 


1 1 :3 

2 1 :4 

3 1 :5 

4 1 :6 

1 :8 


3.0 8 

3.0 7 

9.0 5 

9.0 4 


58 min 30 secs 
52 min 24 secs. 
50 min 27 secs 
41 min 23 secs 
34 min 53 secs 


5 


6.0 


3 
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MORTAR lOADING 

TABLE 4.4^ FIRST GRaGK AlTD FAILURE LOAB 

1 :6 COBiPRBSSira 


Specimen 

number 

hA 

First 
crack 
load in 
concrete 
in tonnes 

First 
crack 
load in 
brick- 
work in 
tonnes 

Failure 
load in 
tonnes 

Remarks 

11 

0.25 

9.0 

9.0 

18.0 

Brickwork at top and 
bending reinforcement 
reaches the yield 
stress 

4 

0.33 

12.0 

12.0 

21 .0 

Bending reinforcement 
at mid span has 
yielded 

12 

0.40 

15.0 

15.0 

24.0 

Brickwork at top and 
bendin,g reinforcement 
mid span has yielded 

13 

0.50 

21.0 

21.0 

30.0 

-do- 

14 

0.80 

24.0 

24.0 

45.0 

Concrete at the supports 
has crushed and bending 
reinforcement has 
yielded 
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i/IORTJ^ LOiiDIi^G 

TABIi) 4.5: COi'tPUTER TIME 

1 :6 COt/lPRESSITE 


Specimen 

number 

H/L 

Initial 
load in 
tonnes 

Number of 

subsequent 

load 

inbrements 

Computer Time 

11 

0.25 

6.0 


35 min 40 secs 

12 

0.40 

12.0 

4 

56 min 57 secs 

13 

0.50 

18.0 

4 

58 min 54 secs 

14 

0,80 

21.0 

8 

96 min 41 secs 
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Table 4.6: qrack Am failure load i.OAi>iRa 

Q«Q 1 >6 QOBaPRE.qsnm 


Specimen 

number 


Type of 
opening 


First crack 
load in 
tonnes 


Failure 
load in 
tonnes 


Remarks 


19 

Symmetric 

window 

24.0 

45.0 

20 

Symmetric 

door 

24.0 

45.0 

21 

Unsymmetric 

window 

18.0 

24.0 

22 

Unsymmetric 

door 

3.0 

9.0 

23 

Rusymmetric 

door 

16 cm R.C. 
beam 

3.0 

9.0 


Concrete at 
supports has crus 
and bending rein- 
forcement has 
yielded 

-do- 


Shear failure 


-do- 


-do- 
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Fig.A-lOa Longitudinal stress distribution at various 
cross sections 
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Fig. 4-10 d Variation of stress in bending reinforcement 
along spon at various loads 



Distance in cm from (eft end of woU 

r- ^ 100 1 C 




,H/l 1 Manor kVn.hnJ-nrt,, 

Cnntprpr, ,ivp qt j 


i i 


Fig.4'10e Vertical deflection of wall along span 
at various loads 
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Fig.4-12e Vertical deflection of wall along span 
at various loads 
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Fig. 4-13 Load vs deflection curves 
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Fig.4'16a Longitudinal stress distribution at 
various cross sections 
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Fig. 4-16 c Shear stres'i distribution at various 
cross sections 
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Fig. 4-18 d Variation of stress in bottom reinforce 
ment along span at various loads 
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Fig.421a Longitudinal stress distribution at 
various cross section 
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Fig.4-21 b Vertical stress distribution along 
span at various heights 
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Fig.422 b Vertical stress o-y distribution along span 
at various heights 




Fig.4-22 d Variation of stress in bending reinforcement 
along span at various loads 
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Fig. 4*23 a Longitudinal stress distribution at various 
cross section 
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dt various heights 


21T (Symmetric door 



Fig.4’23 c Shear stress distribution at various cross sections 
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Fig. 4*23 e Vertical deflection of wall along span 
at various loads 
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Fig. 424 b Vertical stress redistribution along span at 
various heights 
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SFi'i'VLlRY, OOITCLUSIOFS iHD Rii)COM,iEi©iill]IOiTS 
5 . 1 SUMtiARI 

Earlier investia^ators have established that brick 
masonry on reinforced concrete beams is not an overburden 
but behaves in a composite manner providing additional 
strength to the composite system through arching action. 

However benefit of such an interaction has been recoimaended 
only for composite construction of height to span ratio 
greater than 0,5 subject to compressive loading. This is 
probably because sufficient arching action does not come 
into play for height to span ratio below 0,5. Moreover, 
in this range the brickwork tries t'o slip over the reinforced 
concrete beam, Purther more, it has been reported that such 
composite construction subject to tensile loading enjoys 
no benefit of composite action. It is obvious that under 
tensile loading, not only the adV'intage of arching action 
is absent but also the supporting reinforced concrete beam 
separates from the supported brickwork uruTer tensile loadi:ig. 

The present effort demonstrates 'chat such composite construction 
by use of single legged 2-shaped vertical connectors at 
equal spacing all along the length of the composite system 
behaves as a single composite structural element with much 
increased interaction under both compressive and tensile 
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loading. This has been established, in the present work, 
both by experimen al studies and analytical verification. 

The resul-tsfrom the present study lead to certain 
conclusions which are grouped in section 5.2, A design 
procedure is recommended for the design of such composite 
structural element in section 5.3. A set of recommendations 
to extend the scope of present study is presented in section 
5.4-. In any case the present work brings a consciousness 

that inorder to be able to make rational recommendations 
on the composite action of such elements extensive experimental 
and analytical work is warranted. 

5.2 COITGIUSIONS 

1 . Single legged 2-shaped vertical connectors of 
plain reinforcing steel available in minimum size 
i.e. 6 mm diameter provided at a uniform spacing of 
one brick length make the brickwork and supporting 
R.G. beam to act as a single composite structural 
element , 

2. Substantial interaction between the brickwork and 
the supporting R.G. beam in presence of single 
legged Z-shaped vertical connectors as described 
in (l) above exists. In other words, the composite 
structural element is much stronger and stiffen than 
the one without the vertical connectors. Howerer, 



the extent of this interaction depends upon various 
paraxaeter e.g. the mortar strength, the height to 
span ratio and the type of loading. 

Richer is the mortar mix used f >r the brickwork in 
the composite element, more is the interaction i.e. 
higher is the load carrying capacity and lower are 
the deflections. 

Higher is the height to span ratio of the composite 
element more is the interaction i.e. higher is the 
load carrying capacity and lower are the deflections, 
for composite element under compressive loadin/^ 

(load transferred at the top of brickwork) e,g. 
verandah beams, grade beams etc. use of 1:6 cement 
sand mortar mix gives a load factor of more than 
2.0 and hence is recommended for height to span 
ratio of 0.33 or more. However, for height to 
span ratio below 0.3'^ use of richer mortar mix is 
recommended to achieve a lo.ad factor of 2.0 or 
above . 

for composite elexnent under tensile (load transferred 
directly to the supporting reinforced concrete) e.g. 
the multistorey framed buildings, use of 1:3 cement 
sand mortar mix is necessary irrespective of height 
to span ratio inorder to achieve a load factor of 
2.0 or above. 



ise of 6 mm dxameter plain reinforcing steel vertical 
connectorr is sufficient to achieve the desired 
interaction (load factor of 2.0 or above). Increasing 
the di,em-3ter of connectors does not help in increasing 
the load carrying capacity of the composite system 
under either compressive or tensile loading, 
ii'or composite element under compressive loading, 
the uniform spacing of the connectors at one brick 
length i.e. around 25 cms centre to centre is 
sufficient to achieve the desired interaction 
irrespective of mortar strength and height to span 
ratio. 

For composite element under tensile loading, the 
uniform spacing of the connectors at one brick 
length is sufficient to achieve the desired interaction 
for height to span ratio of 0.35 more. However, 

for height to span ratio below 0.33 which may not be 
frequently occurring in practice, closer spacing 
shall be necessary. 

Normal size door and window openings considered in 
the present work symmetric about the vertical centre 

If 

line of the composite element do not affect its load 
carrying capacity under compressive loading. 



-‘he load carrying capacity of a composite element 
having a symmetric door opening and subject to 
tensile loading is adversply affected and use of 
thicker supporting reinforced concrete beam is 
necessary. 

The load carrying capacity of a composite element 
having a symmetric window opening and subject to 
tensile loading is marginally affected. However, 
desired interaction is still available if the sill 
level IS at a height eo,ual to a third of span or 
more. 

Unsymmetric window openings do affect the load 
carrying capacity of the composite element 
signif icantly under both types of loading. However, 
desired interaction is s till available if the sill 
level IS at a height epual to a third of the span 
or more for both types of loading, lor purposes of 
design, the height upto the sill of window opening 
should only be c-onsidered for computing the height 
to span ratio of the composite system. 

Unsymmetric door openings afCect the load carrying 
capacity adversly. Experimental work show that 
thicker supporting reinforced concrete beams are 
necessary if unsymmetric door openings are to be 
used in composite elements under compressive loading. 
The same obviously follows for the tensile loading. 
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15. "file computer pro-iraoime developed in the present 
work is suitable to predict the load response 
characteristics of composite element under 
compressive loading. 

16. First crack load and the failure load analytically 
predicted by the use of finite element method for 
composite elements under compressive loadiaag are 
conservative and hence can be safely used for 
design purposes. 

5,3 ISlSIGiJ iiEThODOLOGY FOR COMPOSITE ELEICEOT 

Oarlier investigators have recommended a suitable 
reduction factor in the design moments to design the supporting 
E.O. beam inorder to account for the benefit of interaction 
in the brickwork supported on R.C, beams. However, these 
recommendations vary markedly. The reason for these variations 
is the variation in mechanical properties of brickwork and 
the dimensions of supporting beam used by these investigators 
in their experimental work. In the present v/ork, it has been 
shown that the thinnest possible supporting beam (thickness 
being just sufficient to provide suitable cover to bending 
reinforcement and hold the single legged Z-shaped connectors 
in position) interacting with the brick masonry over it in 
presence of vertical connectors behaves as a single composite 
structural element. Therefore, for the pvtrposes of design 
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the composite system can he assimied as an under-reinforced 
A.B. (reinforced brickwork) beam provided the depth of 
brickwork is more than 20 percent of the span of the composite 
system. This limit gets fixed from the fact that a balanced 
R.B. section has a height to span ratio of approximately 0.2 
for normal loading as considered in the present work. The 
depth of the composite system corresponds to the height of 
the composite element. In that event the parameters 
required to be fixed while designing are (l) bending 
reinforcement (2) spacing of 6 mm, dia single legged vertical 
connectors and (3) the mortar mix to be used for brickwork. 

The latter two are to be fixed from the conclusions summarised 
in the preceding section. The method for the determination 
of bending reinforcement follows the well known procedure 
of designing an under-reinf orced section. However, the 
maxim-um value of the lever am taken in the present v^ork 
has been limited to 0.4-5 B where 1 is the span of the 
composite system. 

5.4 PROPOSiiii POd PURTHBR E-STBiSIOH OP .PRBSB^fT WORK 

Based on the knowledge gained during the present 
work, the following are the recommexidations for further 
investigation. 

1. Effect of the variation in the height of connectors 

should be studied experimentally as well as analytically 
for both type of loading conditions. 



The bond slip between brick walls and supporting 
R.O, beams and between brick layers themselves has 
bean noticed under tensile loading during the 
experimental work. Therefore, bond slip or linkage 
elements should be included in the finite -element 
model so as to predict correctly the response of 
composite system under tensile loading. 

In the present work the load is either compressive 
or tensile. However, in real life e,g, for multi- 
storey framed construction, a combination of two 
types of independent loadings considered in the 
present work occur. Therefore, study should be 
undertaken by simultaneously applying in suitable 
combinations the load at the top of brickwork and 
at the level of R.C, beam at the first stage. 
Subsequently, even the self weight of the brickvrork 
can be analyxically modeled as a body force varying 
with height. 

In real life situations composite elements shall not 
be simply supported. On the contrary it will be a 
composite infilled frame and should be studied an 
such. In the latter case the effect of unsyrametric 
door and window openings is not expected to be so 
adverse as observed in the present work. 



l3.teral loading to simulate the earthquake and wind 
load effects should be studied through experimental 
and analytical investigations on composite infilled 
frames. 

ihe present and proposed studies should be repeated 
for other mortars e.g. lime surkhi mortar, lime sand 
cement mortar, poztolana cement sind mortar etc. in 
varying thicknesses of brickv/ork. 
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Ntl6)=L5*2 

NCll J»L6»2-1 

Na2)sL6*2 

IEMZ=12 

G0T0326 

IEMZa6 

pg325Jal,IEMZ 

D0325Kal,IEMZ 

KKsNCK)-JJ+I 

IF(KK,LE,0)GOT0325 

CONTINUE 

IFdTER.NE.UGO TO 390 

jl?NSKNH?lMf3r“""*''”™BPLIED''coKci»TRMiD"LOADS TO DP 

DO330Xai , n6ld 

KsNQDLODCl) 

DP(2»K-l)aXLDIN(IdNCB)+0P(2»K»l) 

iSrioS ™ EOUIVALENTS 

ISINODCL) 

JaJNODCL) 

KjKNJpiu 

JJa24a 

KK*2»K 

DXaCXtKD-CXCI) 

DyaCYCKI^CyCI) 

AL«SQRT(OX*»2+Dy**2) 

PXT=SURTRX(L,l)*Ar» 

PXK=SURTRX(L,2)*AL 
PYtaSURTRYCLj) I^AL 
PYKSvSURTRYa,2)»AL 
DP(II-n=DP(Il"U + PXI/6,0 


pp(JJ-n=DP(JJ-l) + (PXl + PXK)/3,0 

DPfKK-l)sDP(KK-l)+PXK/ 6,0 

DP(TI)=UP(TI)-fPyi/6.0 


DP(TI)=UP(H)-»-Pyi/6,0 

DP(JJ)sDP(JJ)+(PYI+pyK)/i,0 

DP(KKJ=DP(KK)+PYK/6,0 

continue 
CALL GEUMBC 
RETURN 
END 


SUBROUTINE INSTMD C IFAIL/ ICONVG) 

COMMON /AREA1/NNP,NEL£N01.D,N0NBC,NMAT,NB0DY,NSL,C 
COMMON/AREA2/Ej.ANU,WTPUV ««» 

COMMON/AREA3/KQDEL,LOC „ ^ 

C0MM0N/AREA6/yiELDS,YXBL,STENBl,YYBT,STEN82,ySHR,STENSN,ECU,FCP, 
lYSTCf DSTRN 

C0MM0N/AREA7/L1C,L2C/L3C,LBC 
CQMM0N/AREA8/NDF,IBAND,1NCR,ITER 
COMMON/AREA9/QtD5,DP 
CDMM0N/AREA12/ISUBRC,NSR,NSRB 
C0MM0N/AREA13 /IjAMDA,B.AREA,BT 
C0MM0N/AREA14/ALtTM,TAT,N 
CUMMON/AREAIS/TSIG.TSTR 
C0MM0N/AREA16/TH,ERPSUD,ERDISP 
C0MM0N/AREA17/C1B,CIC,CN, ANGLE 

common/areai9/esx,esy;esxy,esz 


30 

30 

)0 

)0 

)0 

)0 

)0 

)0 

)0 

10 

10 

)0 

)0 

to 

lO 
• 0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

D 

i 

3 

3 

) 

3 

3 

3 

3 

3 


500 


505 

506 


510 

§12 


515 




- ,^„ aMDA ( 12 . 123 » 

r 


REWIND ITAPEl 
IFAIL=0^ 

IC0NVG»0 

ihdNCR . EQ . D.AND 
GOT0506 
DQ505J«lrNEL 
D0505K*! ,NSR^ ^ 
TSIGCJ#Kfl)»0,0 
ta,K,2)a0,0 
a, K, 3)30,0 

j ' kIDsOIo 

J,K, 2)30,0 
| ajK , 3i . o.o 


(ITER,EQ,1))GOT0500 


TSIG 

TSIG 

TSTR 

TSTR 

TSTR 

CONT 


D0595IsljNEL 
MC»K0DELU) 

LlaLOC(Id) 

L2 » LgC ( l4 ) 

L3sLOCa , 3 ) 

L43L0C(I,4) 

L5sl,OC(I,5) 

N(4)aL2*2 , 

N(5)3L3*2*1 
N(6)sL3*2 

If!hC,EQ. 3)GOTO510 
N(7)»L4f2-l 
nU )« u 4*2 
N(9)«L5»2-1 
NhO ) aL5*2 
N(1U»L6*2-1 
N(l22»G6»2 
JfL3l2 
G0T0512 

C^NDENCE ELEMENT NODALDISPLACEMENT 
D05l5a»lf JL 
JjsNCJ) 

SQCJa )= t * WCJJ ) 

CS.'/ilNUE ron 

COWipEi LiE^T -^.TPH-tuS AT THE TH<EE\CORNER 



FROM GLOBAL NODAL DISPLACEMENT 


NODES OF TRIANGLE 


REGION OF ELEMENT 


DSTRN(li {)sESTRNU »1)^LIC(J 
DSTRN<2, 1)=ESTRN(4,1) 

DSTRN b » n rESTRN ( 7^1 K 


♦LlCjJ 


__CCJ)+ESTRM 
*,>fL2C(J)+ESTRN 
i5♦LiC^ J)+ESTRN(8ii)^L2CCJ)+ESTRN 
CENTROID OF ska ELEMENT 


i+estrn{2,1)»l; 
|♦ ESTRN ( b4 ) 




516 
C 

517 


522 

C 


COMPUTE' STRESSES 

CALL^M|?AuL(CIB,DSTRN,DSlGf i»3» If 3,3r3»lr3fl) 

§! iuKriE “ T ^ liL ' Adu ' STl < il * i ' A *' inf^^ENlROID Of SUBRf CION 

TSTR(l»jlK)=T5TRClf JfK)|OSTHN(K,U 

TSIG(lfiJfK)=T5IG(lf Jf K)+DSIGtR, 13 

CALcSlaIe PRINCIPAL STRESSES AT THE CENTROID OF SUBREGION 


j-rv-u^rs-iiauDHClX iiJJ 

SIGX=T5IG(I,J,n 

SIGXaXSlGaIj;2) 
SIGXY=TSIGCi ,5,3) 
CALL PSTRES 
1F(ISUBRC(I,J).EQ, 
IFCMC.EO.nGOT65« 


CHECK FOR xIELDlNG OF CQNTRP^Tf 

riiL'ipM(?KEfS^Jj2-SIGXWIGyt3.0»(SIGXY»*2) 

IFCFYEL,GE,0.0)GC)T0532 
IFCISUbAc ( i ,5) ,NE,0)GOT0548 

MATRIX CN CORRESPONDING TO CURRENT 

CALL CNE«(I,J) 

A»GtiJX,G)-i»G 

CALL CYXELDCI.J) 

ISUiRC(I,J)=l 

angle(i,5)«ang 


0))GOTU530 

o*(siGxy=»*2) 


igf O’' 

cpcv «ff*e ^ ^ A < 


ESI »isiRa,j;2i 

1F(FCRUS,GT,5,0)GOT0534 ^ 

GO TO 535 
ISUBRC(I,0)a4 
IFAIL=i 

ESS(lTll94l‘^CR.ITER,I,J.ANG 

2»^^Efo^3) * SUBREGION n 6,«» , I2,3X, *HAS CRUSHED* , li, 'A nGLE 

CHECK FOR CRACKING OF CONCRETE 

ICRACI®0 

iCRAC2aO 

1 F C PSX Gl f LE « V « 0 ) §0 TO 545 
IFCPSiGl ,Ge*SXEi^5W)iCRACl»l 

iPRD ® I 

Ft? ClCHACl ) PPlAf TlS^Jir Ir IMCRjfITER. ANG 

F* URmAF ( 1 K 7 % T2 r ^6 f HAS CRACKED IN^fl2> 

At fxOAD INCR^AMENT fio I St^ITiRATION mj,h 

^fF?fsi|iLE.6!J)§gJo546 

IFCPSIGS ,GE,STEWSN) 1CRAC2«1 
I PRDs2 

ISOpRC V I f J J *2 
GQTOSif ^ ^ 
lSU0RC(l,J)a3' 

G0T0548 

IF(ICRACl.E9,UISUBRC(I,0)*3 

5^8!fRE8§^%l 6«ni0536 
DSTRNCl,l)sT5TR(l,J,l) 

DSTRN(2,l)»TSTRa,0,2) 

IF((iPCDSR.EQ,2),AND,(ISUBRC(I , J) ,EQ, 3 ) )G0T0537 
GnT0535 

INSIG(l,n=TSlG(X,J,n 
INSlG(2,l)aTSlGtI, J,2) 

INSIG(3,l)sTSTG(I, J,3) 

GOTU540 

CALCULATE INITIAL STRESSES DUETO YIELDING OR CRUSHING OR DUE TO 
NONLINEARITY OR DUE TO CRACKING OF CONCRETE 
D0533Kal ,3 
D0533LS1 ,3 

CN(K,L)=CIC(K,L)-CN(K,L) 

continue 

CALL MATMUL(CN,DSTBN,INSIG,3,3,l,3,3#3rl»3,l) 

GOT0540 

CHECK FOR YIELDING AND CRACKING OF BRICKWORK 
TF(ISUBRC(I,J) ,EQ,UG0TU6b2 


,J),EQ,3))G0T0537 


DUE TO 


V * ' W * ~ V 

ICRAC2=0 


|^^|^|JS'^«Gt.&TtNB2)ICRAC2sl 


tp?fTrRA?}^^.^‘';5FsOJGOT064b' V. 

IFCaCRAci F§*li AK!n*?TrDAri»&^R*9^^^°'^°550 

G0T0648 

IP'(IPCDSR,f4L0)GOTO654 

BilSSB;} if 

DSTg«l|,n=rsTRh:S:iS 

SL»TSIGa,J,l)/YXBL 


Sti®(fSIG(I-J.l)p*DSTG (1-1)1 /YXBt. 
§§c.i f 1) ^5§IQXI # ii) /¥SHR 


D0915K=1,3 
D0915G=l,3 
CN(KiG)=ESZ»CH(K,Ii3 
lSGBRCCl,J)=i ' 

Pg|gT705aNCR,ITER,I,J,ANG 

CALL C¥IEGD(I, JJ 

D0656K=1 4 
D06S6L=1 »3 

CN(K,L)aCI0CK-L3-CN(K,L) 

CALL M^TMUL(CM,DSTRN£lNSIG/3»3, l,3,3,3,l ,3,1 ) 

So 54 ?kL 3 cihMb'ut iitH'mkaioii - 

TSIGCl7j!K)»ISIGCI,J,K)-INSlG(lt,U 

CALCULATION OF PSEUDO LOAD VECTOR DUE TO NONLINEARITT OF MATERIAL 
ARlAOF^THfsUBlillOi!*^^^'^'^^^® MATERIAL AgrCRUSHING OF MAtMia(. 

AN3SR=NSR 

AREAS=ARRA/ANSR 

Cl'=rh(n»AREAS 

DP5(1,1)=CT»INSIG(1,1)»L1C(J) 

DPS(2,U=CT*XNS1G(1,1)»L2CCJ) 

UPSC3,l}=CT»lNSIG{l,l)»L3Cta5 

DPSCArDsCT^INSIGCZfD^LICCJ) 

DPS(5,n=CT»lNSIG(2,l3»L2C(U) 

DPS(6,l)=CT4INSlG(2,l)»L3CtJ) 

DPS(7,U=CT»INS1G(3,U»L1C(J) 

DPS(8,1)=CT=»1NSIG(3,1 )*L2C(a) 

DPS(9,l)=CT + INSIGt3 U»L3C(a5 

D05b2K=l,12 

D05b2L=l ,9 

BT(K,)U)=B(L,K) 

CALL MATMUL(rtT,0PS,PLV,l2,9,l,12£l2,9,l,12,l) 

ASSEMBLE GLOBAL PSEUDO LOAD viWR »*»**'*' 

D0555K=1 ,12 
L=N(K) 

DP(L)sDP(L)+PLVCK,1) 

CONTINUE 
CONTINUE 
CONTINUE 
GO TO 596 
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oo 
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?0 

50 

50 

)0 
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50 
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)0 

)0 

50 

50 

>0 

)0 

)0 

>0 

>0 

10 

>0 

10 

10 

10 

to 

to 

lO 

to 

lO 

50 

to 

lO 

'0 

0 

lO 

0 

0 

0 

0 

0 

0 
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0 
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0 

0 
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561 

c 


562 

C 


569 

701 


570 

565 

C 


574 

575 
C 


576 


577 


578 


579 


580 

581 


583 

582 


585 

595 

C 


C 

C 


RK An f T Ta PI.’ 1 

Kt* Au V i 1 Ar L I ^ ^ ) AL # C C TM r R . f, ") R***! 1 t 

CALE^MATMUMTM^S0?SOlI°?^fi V LOCAL COORDINATE SYSTEM 


ALL matmul(tm7so;sql,3.6 1 


MC=KODEL(n 


SUB-BAR 




0+2,0*LBC(K))»SQL(2,U 


CONTINUE 
COMPUTE TOTAL 
D0562K=1 ,NSRB 

TSTr f I ' K ' lI-l«Tr r T ' ' 1 3 +ESTRN (K, 1 ) 

CHECK EACH SUB- 
Da570K«l,NSRB 

GOT0570 
IX(K)sl 
IFIISUBRCC; 

FORMAT (2X,ui 
i'BAR ELEMENT 
ISUBRCCl, 

CONTINUE 
D0565J=1.3 
PLVLCJ, 15=0,0 

D0575K=PnSRb'^^^^ stresses DUE TO YIELDING OF STEEL 
INSIGM(k5=0.0 
IFCiy(K5|Eg.OJGpTD575 


STRAINS AND STRESS AT THE CENTROID OF EACH SUB-BAR 


•REGION FOR YIELDING 


.M|NT^NO,= %I3,3X,*SUBAEGforNO:=Ml27lx%AS’nfeLDfep5 


IFCfSXG i7K 1 lEfrO 0)GQT0574 

|NfI?HCK)fe-Y||LO|?JsiG(I,KjU) 

JgJC|^^K,l)*TSIG(I,K,l)-INSIGM(K) 


|NfIGMCK)7TSIGafK,l)tYIELDS 
^oifrN&r ^ ^ 5"TSIGC If K(r 1 j - INSIGM(K) 

calculate PSEUDO LOAD 
iFCNSRB,EQ.UGOT0579 
IFClYdJ ,EQ,0)GOTQ576 
CT=TH(i5^INSIGM(1)^0,5 
PLVLCl fl5=PiiVLQ f 15-1 ,25*CT 
PLVL(2,l5=PLVLC2ll5-0 25*CT 

i,5»CT 


VECTOR DUE TO YIELDING OP STEEL 


pLVLafi5=PLVLCi7i; 
PLVL(2,l)=r ' ' 


PLVLO.i 
IFClYCaJ 

ct=th(i5 

PLVLdfl 

PLVL(2fl 


, sPLVLdfl 
)=PLVL(3fl 
,EQ.0)GO TC 
fiNSIGMO: 
JsPLVLC 
WPLVLI 


cr!:?Hd5fiN§iGM(4)fo, 


W,3 

-0,75»CT 

:yEig 

,50»CT 


PLVLdfnsPLVLI. 
PLVL(2fi)«PLVL(; 
PLVL|3gl)=PLVL(; 


GOTO- . . 
IFClYil 
CT=THa 
PLVLCl,- 
PLVLC2,1 


rl)-l,! 


;S*CT 

1S*CT 

iO^CT 



,f03G0TO580 

pCvEj3;iie§;o*'^* 

D058iK*lfNSR8 

IFCIY(k5.NE.03G0TO583 

continue 

G0T0595 

S*SIEH TO 

D0582K=1 ,6 
DQ5t}2L = l ,3 
TMTCK, L5=’rM(L,K) 

CALL MATMUL(TMTfPLVL,PLV»6f3,lf6,3,3flfl2.l3 

add THIS PSEUDO LOAD VECTOR TO GLOBAL PSEUDO LOAD VECTOR 

DD5d5tV“l ,6 

L=N(K) 

DP(L)=DP(L)+PLV(K,1) 

CONTINUE 

CONTINUE 

CHECK FOR CONVERGENCE 

TCOWV=0 

ICOND=0 

CHECK whEATHER PSEUDO LOAD VECTOR ,8 ELEMENTS ARE LESS THAN SOME 
SPECIFIED LIMIT 


351 


(ICOND, 


,NDF))GOTO600 


1210 


IK UBSCDPCJ ) ) .LE.ERPSUD) lCONVslCONV+1 
CHECK^CONVERGENCE OF DISPLACEMENTS 

^’l^^t^^^I^ISPUCONDsICUNDtl 

continue 

IFClTER,EO,l)GOTOb90 

*AND,(1COND,EO.NDF))GOT0600 

CyOTOo & 

XFCIC0NV,EQ,NDF)G0 TO 600 

GO TO 635 

XCQNVG=1 

CALL PRINTR 

RETURN 

END 

PSTRESS ^ ^ » 

SUBROUTINE PSTRES 

C0MMgN/AREA20/SX,SY,TX¥,SIGl,SIG2,ANG 

XPY=(SX+SY)I0.5 

XNY=(SX-Sy)IO,5 

CT=SQRI(XNVII2+TXYI=I2) 

sigi=xpy+ct 

S1G2=XPY-CT 

IFCCT.Eg.O,0)GOTQ300 

ANG=( 180. 0/6. 2831853 )IACOS(XNY/CT) 

JFaXYjLT,0,0)ANG=180,0*ANG 

ANG=0,0 

RETURN 

END 

♦ ]!(♦«♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦:);♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
^^♦♦♦♦♦♦♦♦♦♦♦♦♦^♦♦f ♦♦♦♦♦♦♦♦♦♦♦!(♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

SUBROUTINE MATMUL (ML,MR,PM ,N,M ,K,N1 ,Ml ,M2 ,K2 , N3 , K3 ) 

REAL ML(N1,M1),MRCM2,K2) ,PMU3,K35 

0012101=1, N 

DD1210J=1,K 

PM(I,g)aO,0 

D01210L=1,H 

PMCItU)=PM(I,J)+MLCI,L)*MR(L,U) 

CONTINUE 

RETURN 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*»l♦♦♦♦♦*l♦♦♦ 

♦♦♦♦♦♦♦♦l♦♦♦♦S^I♦♦l♦♦♦♦♦♦♦♦♦♦♦♦f♦♦l^♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»♦♦♦♦♦♦♦♦ 

C0MM0N/ARiAl/NNplfe^N6LDlN0NBC,NHAT,NB0DY,NSLC 

E8BKrRiJl?^^i°KG „ ^ 

CDMHON/AsiAli/llS.ERfsODjERilSP ryfinnl lilin »(12 131 

5?tl!icxa2, 

A2»CXCLi5-CX(L35 

A3»CX(L2)-CX(L15 

11:85 81 :8nk 

AREA=ot5*(Al*B3-A3*B12 

ASSEMBLY OF NODAL 0ISPLACE.MEN1 MATRIX B 
DOiOOJsl ,3 
DOIOOK*! »6 
B(J,K)-b6(J,K) 

B(a K+b)s0,0 
BCJ+3,K)=0,0 , , 

B(J+3,K+61=riVCJ.K) 

R(J+6,K)=BV(J,K5 

H(J+6,K+6)=BU(J,K) 

CONTINUE 

TRANSFORMATION OF B MATRIX , 

CALL MATMULCil,LAMDA,B’r,9,l2,12#9, 12»i2rl2f I2f 12) 

D0102J=1 f9 
DU102!S.=1 ,12 ^ 

B( J,K)=BTIJ,K) 

TRANSPOSE OF H MATRIX 
DQ105J=1 ,9 
D0105X=l»12 ^ 

BT(K,J)=HtJ,K) 

CONTINUE 



J. ** * 4 . 1 ^ V 9 jUUO#lU w UUni/ 4 .^#% A C» 


553 


SYSTEM 

CALL DMATklfAREA,!) 

CALL CMATHXa) 

Daiioj=i’^i2 
D0110K=ll9 
BTCJ,K)=KGCJ,K) 

D0115K=ij9 

BTCJ,K)=KG(J,K) 

RETURN ^ ^ ^ fl2f9|12fl2fl2f9fl2rl2fl2) 

END 

♦ ♦♦ + 4 : ♦ilf****^!^ »sfi 

* SUBROUTINE ESMAT2 4 

44444444444444»*44 4 * 4 ^^Hc 4 :<(:*! 4 !)j(*t *444 444444 » 44 » 44 *»)|‘ 4444 » 4 * 4444 » 44»4 
SUBROUTINE ESMAT2 ( I.Ll , L2 ,L3 ) 

common/areau/ka,kg 

rMtA di f*! At y A n I » -4 il / A V ^ rmt 


COMMON /AREA14/AL,TM,TMT 
CQMM0N/AREA16/TH,ERPSUD 


r 



,Cy(400>,TM(3,6) ,TMTC6,3), 


ERDISP 
' “,NC 12 ) 

MC»KODELCX 2 

Dx»cxa 2 )^cxcLn 

DY=CYa2)-CYCH) 

AL=SQRTCDX**2+DY*»2) 

CALCULATION OF TRANSFORMATION MATRIX 

COSAs'DX/AL 

5INA=DY/AL 

D0400U=1 ,3 

DQ400K®! .6 

TM(J,R)=6,0 

TM(1 ,1)*C0SA 

TM(1,2)=SINA 

TMC2,i)=C0SA 

TM(2,4)=SINA 

TMC3,&)sCOSA 

TMC3i6)=SlNA 

TRASPOSE OF TRANSFORMATION MATRIX 

004050*1,3 

D0405R*i ,6 

TMTCK K) 

ELEMENT STIFFNESS MATRIX OF BAR ELEMENT IN LOCAL COORDINATE SYSTEM 
CT= CTH ( iT^E (MC , 1 ) ) / C AL»3 , 0 } 

KL(l,l)a7,0kT 
KLUf2)sCT, 

- ^ 2 , 1 ) 5 KLC 1 , 2 ) 



kIiEM Jlel^otcT 

TRANSFORMATION OF STIFFNESS MATRIX^IN GLOBAL COORDINATE SYSTEM 
CALL MATMULCTMT,KL,KG,6,3,3,6,3,3,3,12,12) 

D 0410 J= 1,6 

DD4iOKsl,3 

CALL AAiM(JL{TSTiTM,KG,6,3,6,6,3,3,6,l2,12) 

RETURN 

FMD 

♦ !J!it:t !4 + * 4 !* 444 »» 4 »**t 44444 ** 44 » 44 » 44 * 4 * 44 »» 4 * 44 » 4444 *********** 4**44 

SUBROUTINE BMATRI C AREA,Al ,A2 , A3 ,HM) 

DIMENSION BM(3,6) 

CT=0.b/AREA 
BM(1,1)*3,04A1*CT 
HM(i ,2 )s:-a 24CT 
BM(1,3)=-A3»CT 
BM(i,4)=4.04A2*CT 
BM(l,b)=0,0 
BMCl ,6)a4,04A3*CT 
BM(2,l)=-Al*CT^ , 

BMC2,2)=3,0*A2*CT 

BMC2,i)=-A3»CT 

BM(2,4)=4,04A1*CT 

BM(2 blsA.O^AI^CT 

BM(2,6)=0,0 

BM(3,ns-Al^CT 

BM(3,2)=-A2»CT 


ALPHABsl.O 
AljPriAC = 0,4 
D0500K=1,3 
005000=1*3 
CN(K,0)=0,0 

IFCISUBRC(1,J3 ,EQ.3)G0T0599 
ANG=ANGOE(ltJ) 

THEATA=(ANGV3, 141 59265/180,0 

COSA=CDS(THEATA) 

SINA=SINCTHeATA) 

MC=KQDEOa) 

ir(MC,E0.2)GQT0550 

TRANSFORM THE ELASTICITY MATRIX 

CRACKED DIRECTION 

TM(l,l)=C0SA**2 

TMCl ,2)=SINA*»2 


FROM GLOBAL COORDINATE SYSTEMTO 


Jo? 

300 

)ll - 

300 C 
300 

100 505 

200 
300 


""J-WVi / J LU'd«'+ Si IK 

TM(2,2)=TMaf 1) 

TM(2,3)=-TM(i,3)^ ^ 

TMC3;a)=-2,0*TM(l ,33 
TM(3,2)=-TM(3.1) 

TMC3,33=TMC1,1)"TM(1,2) 

TRANSPOSE OF TRANSFORMATION MATRIX 

DO505Lsl,3 

D0505K=I,3 

TMT(K,t.3*THa,K) , , , , 

CALO MATMUL(TMT,CtB,PM,3,3,3,6,3,|,|,3. 
CAOO MATMULCPM,TM,CN,3,3,3,3,3,3,6,3,33 


33 


rrr- 

)0 

)0 

)0 

50 

DO 

00 

OO 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

■00 

.00 

.00 

00 

100 

»00 

>00 


550 

c 

c 

560 


540 


599 


C 

C 

C 


70S 


mAtion matrix 


700 


C 

C 

C 


CN(1,2)=0,0 
CNCi ,3)s0,0 
CN(2,1)=0,0 
CN(3^)=0.0 

CNC2,2)=ECMCf U 

CN(3,3)=AJjPHAC»(E(MC» 1)/C1 ,0+ANUCMCf 1 ) ) )*0.5 
cdordimate^syIIem^^^ matrix from crAcked system to global 

TMC1,1)=CDSA*^2'^ 

TMII ,2)sSINA**2 
TM(1»3)=C0SA*SINA 
TM(2,l)sTM(l,2) 

TM(2,2)=TM(lJn 
TM(2,3)s-TM(l,3) 

TM(3,1)=-2.0*TM(1,3) 

TM(3,2)5-TM(3,U 
TM(3,3)=TM(ia)"TM(; 

TRANSPOSE OF TRANSFC 
D0540LS1,3 
D0540K = 1 ,3 
TMT(K,L)sTM(L,K) 

CALL MATMUL(THT,CN,PM,3,3,3f6,3.3,3t3,3) 

CALL MATMUL(PM,tM,eN,S,S,SJ,5 

HEiTU HIm 

END 

SUBROUTINE CIMAT 

SUBROUTINE CIHAT 

common/areai/nnp,nel,nold,nonbc,nhat#nbody,nslc 

COMMON /AREA2/E,ANUtWTPUV 
COMMON/AREA17/CIB,CIC,CN, ANGLE 
CQMMDN/AREA18/MATCOD#NCODE 
INTEGER MATCODO) .Nc6dE(400D 

REAL|(3^2)^ANU(3,2),WTPUV(3),CIB(3,3),CIC(3,3),CN(3,3), 

^DO700Isl,^MAT 

MC*MATCaD(I) 

IFCMC,EQ,35gOTO7O0 
IF(MC.ES,l)GgT0705 ^ 

CTaE(MC.i)/a.0-ANU(MC,U**2) 

ClC(lfi5®CT 

gig 

CIC(3»2)®0,0 

CIC(3jt3)sO,5*Cl*0-ANU(MCfl))*CT 

GOTO700 

Slol SsicTlCE (McS l)+E?MC^252ijoUQRT(ANU CMC, t)*ANU(MC, 2) ♦£(«€, !)♦ 
lECMC “ 

CIB( 

CIB( 

CIB( 24 )s»ci 8 (l, 2 )„ 

CIB(2,2)5CT*ECMC,2) 

CIB(2,3jaO,0 
CIBC3,l5a0,0 
CIBC3,25«0,0 
CIB(3»3)aG 
CONTINUE 
RETURN 

5 ftt SURRniTTTNE CNE*W ^ 

SUBROUTINE CNEW(I,J) 

COMMON /AREA2/ErANU,»<TPUV 

CUMK/ARiA6/5lELD^!^?m,STENBl,YYBT,STENB2,YSHlR,STENSl*,ECU.FCP, 

1Y5TC,DSTRN 

COMMON /AREA 1 5 /TSIG,TSTR 
COMMON /AREA17/CIB,CIC,CN, ANGLE 
COMMQN/AREAn/EX.Ey.EXY^EZ , 

INTEGER K0DEL(240)iL0C(240,b, 



aCT*E(MC#l) 

aCT*E(MCfU^ANUCMC,2) 


igji:"^f 3 ^^^^wf.narrifuyni.DSTfN( 3 ,n,m 

13),CIB(3,3) ,CICC3,3),CN(3,3),ANGLE(240,4) 
MC=KODELU) 


EX=TSrR(I,J,l) 
EX=TSrP(I,J,2) 
EXyaTSrHCl.J,3) 

VsANlKMC »n 

CALCULATE EQUIVALENT STRAIN 


'9200 

'9300 

i9400 

9500 

,9600 

9700 

9800 

9900 

OOOO 

0100 

0200 

0300 

0400 

0500 

0600 

0700 

0800 

0900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

,2400 

,2500 

i2600 

12700 
12800 
L2900 
13000 
13100 
13200 
13300 
■3400 
3500 
3600 
00 
00 
00 
4000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
4800 

5000 
5100 
5200 
5300 
5400 
5500 
5600 
700 
800 
900 
000 
'iOO 
. ,'200 


EMN=FCNi 000, 0*(1 ,0-500, 0*ESAV) 

EZ=EMN/C1,0-5'**2) 

CNC1,11=EZ 
CNa,25 = V*EZ 
CNC1,3)=0,0 
CN(2,l)sCN(l,2) 

CNC2 2)=EZ 
CNC2,3)aO,0 
CN(3ri)=0,0 
CNC3, 25=0,0 

END 

cyiELD ♦ 

SUBROUTINE CYIEliDCl.J) 

C0MM0N/AREA2/E,ANU,WTPUV 

COMMON/AREA3/KODEL,liOC 

COMMON/AREA6/YIELDS,YXBL,STEHBl,YYBT,STENB2,!fSHR,STENSN,ECU,FCP, 

lYSTCfOSTRN 

COMMON/AREA12/ISUBRCj.NSR,N5RB 

CQMMON/AREAIB/TSIG.TSTR 

COMMON/AREAi7/ClB,ClC,CN, ANGLE 

COMMON/AREA19/EX.EY,EXY,EZ 

INTEGER KODELC2403j.LOC(240,61 ,ISUBRCC240,43 

REAL Bp,2),ANU(3,2) ,WTPUV(3) IdSTRNU.U ,TSIG(240,4,3 J ,TSTR(240,4, 
13),CiB(3,3),CIC(3,3),CN(3,3),ANGLE(240,4) 

D0650K=1 ,3 
D0650Lal,3 
CN(K,L)s6,0 

IF(k6dEL(I3,EQ,1)G0T0690 

AMaO.O 

IF(ISUBRCCI,J),NE,0)GOT0690 

MCaKODEL(I) 

EX*TSTR(I,J,n-DSTRNCl,l) 

EY»TSTR(I,J,2 5-DSTRN(2,n 
EXYsTSTRCl/Jraj-DSTRNcS,!) 

VaANUlMC^U 

EZ»*CV*(EX+EY))/C1,0-V) 

ES1 = (CEX-£Y)*4'2 + CEY-EZ)*»2 + (EZ-EX)*»2+1,5*CEXY»»2)) 
ESt=SQRT(ESl/2,0)/(l.0+V) 


ES=(CEX-EY)*»2+(EY-EZ)»»2 

ES=(SQRTCES/2,0))/(l,0+V) 


•EX)**2 + 1.5»=(EXY**2)) 


|6500 


m 

llioo 

||900 

17000 

w 

■fj 


54(EXY*»2)) 


,0*ESAV) 


,5=t<C1.0-V)*EZ 


+ JI:f2Ezi4=*2 + tEZ-EX)*»2 + l,54(EXY**2)) 


EXaTS|Ra,0,l) 

EYaTSTRCI/Jt^J 

EXyaTSTRCLj»3 

EZs-CV»(EX+EY) 

ES2=((EX-EY)*» 

ES2»SQRTCES2/2 


0)/(l,0+V) 


vAfV/Tf/ 

AMa(YSTC-ESl)/CES2-ESl3 

D065|K=1,3 

D0655Lai,3 

CNCK,L)=CIC(K,L)*AM 

RETURN 

COMMON/ARiAl/NNp%EL,NOLD,NONBC ,NMATrNB0DY#N5LC 
C0MM0N/AREA8/NDF,IBAND, INCH, ITER 
COMMON/AREA9/0,DQ,DP 

COMMON /ArIaI 8 AMAfCQD , NCODF, 

800),^ 

IF((INCR.E0,1) . And; (ITER, t;0,l))G6T0400 

APPLY^DISPLACEMENT BOUNDRY conditions in stiffness MATRIX 

D0405I=1 ,NNP 

IFCNCODECI) ,EQ. 030010405 

IF(NCODE(I),EQ,23GOT0410 

KA(2*i-l,l)=KA(2*l-l;l )*0,IE15 

IF(NCODE(I),EO,1)60T6405 ^ 

KA(2*I,l)aKA(2»I,l)»0.1E15 

APPLY^ofsPLACEMENT BOUNDRY CONDITIONS IN LOAD VECTOR 
D0420I=1.NNP _ 


IF(WCODEh),EQ,O3GOTO420 
IFCNCODECI) ,Eg, 2 5GOT0425 


4000 

4005 


4010 


4015 


4025 

4020 

4030 


4070 


4040 

4036 

4035 

C 


4045 

4050 

4055 


^ TiA ^ -u ,y 

XF(mCUDECI),EQ.13GOTO420 

DPC2*1)=0,0 

CQHTINUt; 

RETURN 9M1 

END W 

* SUBROUrUNE PRINTR ♦ 

SUBROUTINE PRINTR 

COMMON /AREA1/NNP,NEL,N0LD,N0NBC ,NMAT,NB0DY,NSLC 
COMMON/AREA3/KODEL.LOC ouua,i cuu 

C0MM0N/AREA8/NDFtIBAND,INCR,ITER 

COMMON /AREAR/OtDQ, DP 

COMMON /AREAl2/ISUBRCf;NSR,NSKB 

C0MM0N/AREA15/TSIG,TSTR 

COMMDN/AREAIS/MATCOD.NCODE 

COMMgN/AREA20/SIGX.SIGY,SlGXY,PSIGl,PSIG2.ANG 

INTEGER MATCOD (3 ), NCgDE(400),RODEL( 240 ), l6c (240, 61, ISUBRC (240,4) 

REAL^^(8OO),DQ(8OO),0P(BOO),TSIG(24O,4,3),T5TR(24O,4,3) 

F0RMAT(/1X,13I(1H*)) 

PRINT4005 tUCR ITER 

FORMAT C/sL'INCREAMENT* NO.aM3/5X,17ClH-)/22X, 'ITERATION NO.sMi 
1 / 

PRI«T40ig,(I,Qp»I-l) ,QC2»I) ,DP(2»I«1) ,DPC 2 »I),I® 1 ,NNP) 

FORMAT (/2X, 'NODAL DISPLACEMENT OF THE STRUCTURE AND PSEUDO LOADVE 
ICTOR IN global COORDINATE SYSTEM'/2X,87aH-)//2X. 'NODE NUMBER'. 8X, 
rU=X-DISPLACEMENT',8X,'VsY"DISPLACEMto%5X, 'PSEUDO LOAD IN X-.6lR^ 
3CTI0N',5X, 'PSEUDO LOAD IN Y-DIRECTION ' / , (5X,I4, 12X,E15 ,8 , 9X,E15,8 , 
412X,E15*8 ,loX,E15«8) ) 

PRINT4015 

F0RMATC//2X, 'STRAINS AT THE CENTROID OF EACH SUBREGION OF ATRIANGU 
ILAR element V2X,66(1H*)//5X, 'element NO,’ ,3X, ^SUBREGION NO,',8X, 
2'STRAIN EX', 12X, 'STRAIN EY' , 12X, 'STRAIN EXY' ) 

D04020lal,NEL 

IFCKODELCI),EQ,3)GOT04020 

PRINT4025,I,CJ,TSTR(I.J.l),TSTRei,J,2),TSTRCl,J,3),J=l,N5R) 

FORMAT(/8X,I4,l2X,I2,13X,El5,8,6X,Ei5,8,6X,ElS,8/(24X,I2,13X, 

PRINTtoSo'"^’® " 

F0RMAT(//2X, 'STRESSES AT THE CENTROID OF EACH SUBREGION OF A TRIAN 
IGULAR ELEMENT'/2X,67(lH-)/lX.’ELEHEMT',2X,*SUBREGinN',2X, 'SUBREGIO 
3ION' ,4X, 'STRESS SIGMAX',4X, ’STRESS SIGMAY' ,4X, 'STRESS SIGMMY/,3 
3X, 'STRESS SIGMA1',4X, 'Sites SIGMA2 ' ,6X, 'ANGLE PHI*/3X, 'NO, ' ,6X, 
4'NO, ',4X, 'C0NDITI6 n'/22X, 'CODE') 

DD4035Isl,NEL 
IFCKODELCI) .EQ,33GOT04035 
PRINT4070 
FORMATC/) 

D04036Jsl,NSR 

SIGX=TSIG(I,0,1) 

5IGY»TSIG(I,J,25 

SIGXY=TSIG(I,J,3) 

CALL) PSTRES 

PRINT4040,I,J,ISUBRC(I,J),SIGX,SIGy,SIGXY,PSlGl,PSIG2,ANG 

FORMAT(2X,I4,6X,I2,9X,Il,6X,E15,8,2X,E15,8r2X,E15,8,2X,E15,8,2X, 

1E15,8,2X,E15,8) 

CONTINUE 

CONTINUE 

STRESSES AND STRAINS AT THE CENTROID OF EACH SUB-BAR 

0040451=1,. NmaT 

IF(MArCUDCI),EQ,3)GOTO4050 

CONTINUE 

G0TU4095 . 


4065 

4060 

4095 


lOF A BAR ELEHENT'/2X,70C1H-)//5X,'E:LEMENT NO,',^X,:PU»RKUiU« NU,;, 

22X, ^SUBREGION CONDITION CODE',3X,' AXlALSTRAINEX' ,5X, 'AXIAL STRESS 
3SIGMAX') 

D04060I=1,NEL 

PRINT^ObsH J(jIlSUBRCClAJ)fTSTR(I,a,l) ,TSIGCI,J,i) , J»1 ,NSp) 
FaRMAl(/8X,|4,{0X,I2,l7i,Il,14X,El5,8,7XrEl5,8/(22X,I2,l7X,Il#14X, 
lE15,8,7X,Ei5,8)) 

CUNIINUE 

PRINTAOOO 

RETURN 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦* 

SUBROUTINE DSOLVE 

C0MM0N/AREA8/N,IBW,INCR,ITER 

C0MM0N/AREA9/0,X,B 

REAL^QctoO) J b(86o) ,X(800) ,A(800,56) , KG (12, I 2) 


rr^ 

<000 

5 

'|09 


'200 


300 


400 


500 


600 


700 


800 


900 


000 


100 


200 

US 

300 

110 

400 


500 


600 


700 

120 

BOO 

150 

?00 

100 

)00 

10 

100 


!00 


100 


too 


»oo 


iOO 


fOO 


100 


100 

165 

100 

160 

.00 

155 

•00 

[oo 


100 


ioo 


lOO 


00 


100 


100 


100 


00 


!00 

180 

iOO 

175 

00 

170 


125 


126 

00 



D01001 = 1 ,N 
IP=M«I+1 

TFC|Bi^,LT,IP)IP=lBW 
DOlSOJsiflP ' I '* 

IQsXBwj 

SUM=A(I, Jj 

IFCie.K9.0)GOTOUO - 
D0115K=1 ,IQ 
IK=I«K 
JK=J+K 

IF(J.NE.1)G0T0120 

jntS 

GDTOi&O 

CONTINUE 

D015SI=1,N 

J=I-IBrf+l 

IKsi-K+l 

SUM=!SUM-ACK,IK)*X(K) 

CONTINUE 

X(I3=SUM/A(I,1) 

CONTINUE 

001701*1, N 

TRsN+l-I 

JaIR+lBW-1 

IF(J,GT.N)JsN 

SUMaX(IR) 

IMlaIR+1 

IK=K-IR+1 

SUMsSUM-AdR^IKl^rXCK) 

CONTINUE 

X(IR)=SUM/A(IR,1) 

CONTINUE 
RETURN 

PRINT126,5UM 
F0RMAT(/2X, ^MATRIX 15 
STOP 
END 


340 


NOT POSSITIVE DEFINITE V/2X,*SUM»%E15, 8) 



